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Chapter 1 
Introduction 



The huge field of partial differential equations (PDEs) is a rapidly devel- 
oping branch of mathematics which has a countless number of applications 
in natural science, engineering, economics, etc. The formulation of physical 
theories would be entirely impossible without the language of partial differ- 
ential equations. A mathematically rigorous investigation of PDEs started 
in the 20 th century although many corresponding physical theories are much 
older. While in the old days the studies were mainly confined to the con- 
struction of explicit solutions (e.g. in terms of integrals), a paradigm shift 
took place with the introduction of functional analytic methods in the 20 th 
century. It has been realized that for most equations it is hopeless to try 
to construct explicit solutions and even if this is possible, the arising formu- 
las are completely useless due to their complexity. Thus, the new methods 
focused on studying properties of solutions, function spaces they belong to 
and estimates they satisfy without knowing them explicitly. This approach 
turned out to be very fruitful and the highly successful modern PDE theory 
is entirely based on it. 

A special class of PDEs are so-called evolution equations where, roughly 
speaking, one of the independent variables is interpreted as time. For that 
kind of equations one considers the Cauchy problem, i.e. one prescribes ini- 
tial data at a certain instance of time and studies the future development. 
Typical results are related to the existence and uniqueness of solutions as 
well as their dependence on the data. In this respect one has to distinguish 
carefully between local (for small times) and global (for all times) existence. 
It may well be possible that the Cauchy problem is locally well-posed but 
it does not admit global solutions, at least not for arbitrary data. This fact 
leads to the problem of formation of singularities, i.e. the breakdown of so- 
lutions in finite time. Once one has a local existence result one can ask the 
question: Given regular data, is it possible that the corresponding solution 
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ceases to exist after a finite time in the future and if so, how does this hap- 
pen? Closely related to this is the question under what circumstances the 
solution extends globally in time. It is worth mentioning that the possible 
breakdown of solutions (singularity formation) is a feature which is com- 
mon to many nonlinear evolution equations arising from various branches of 
physics, chemistry, biology, etc. 

This thesis is devoted to the investigation of some aspects of the Cauchy 
problem for wave maps from Minkowski space to the three-sphere S 3 . The 
wave maps system, a nonlinear generalization of the wave equation, arises 
from a geometric action principle. Nowadays, local well-posedness for the 
Cauchy problem of the wave maps system is well-understood. Furthermore, 
some major progress concerning global aspects has been made quite recently. 
Some of these results are reviewed. However, in this particular model of 
wave maps with S 3 target, singularity formation is possible. One can explic- 
itly construct self-similar solutions that "blow up" at a prescribed finite time. 
Imposing additional symmetry assumptions, this wave maps system reduces 
to a single semilinear wave equation which exhibits a rich phenomenology 
despite its simplicity. In particular, it shares some properties with Einstein's 
equations in connection with gravitational collapse. Thus, this semilinear 
wave equation is used as a toy model for the much more involved evolu- 
tion problem in general relativity since it is simple enough to be accessible 
for analytic techniques. The main goal of this thesis is to obtain a bet- 
ter understanding of singularity formation for wave maps from Minkowski 
space to the three-sphere. To this end, spectral theoretic aspects of self- 
similar solutions and their linear stability are studied. The mathematical 
machinery required for this consists of functional analysis, operator theory, 
Sturm-Liouville theory and semigroups. Since this thesis is intended to be 
essentially self-contained, all the mathematical requirements are carefully in- 
troduced although not all results are proved since this would go beyond the 
scope of this work. However, if the proof of a certain result is omitted then 
there is given an easy accessible reference or at least the idea is sketched. 
The organization of the thesis is as follows. 

In ch. [2] we define the wave maps model. First, we recall some basic 
notions from differential geometry required to state the action functional. 
Then, we derive the associated Euler-Lagrange equations which constitute 
the wave maps system. This system reduces to a single semilinear wave 
equation under the assumption of co-rotationality. Finally, we give some 
historical remarks and applications of wave maps. 

In ch. [3] we collect some basic mathematical concepts. We introduce 
Banach, Hilbert and Lebesgue spaces and note some elementary properties. 
Then, we briefly recall the basics of linear operator theory on Banach spaces. 
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Next, we turn to some aspects of spectral theory for closed operators and 
in the last part of this chapter we give a short overview of basic semigroup 
theory. 

In ch. [4] we focus on Sturm-Liouville operators. First, we state some 
general properties of self-adjoint operators. Then, we turn to the study 
of Sturm-Liouville operators. We discuss endpoint classification, calculate 
adjoints for regular and singular operators and state the Weyl alternative. 

In ch. [5] we apply semigroup techniques to abstract wave equations. By 
this we obtain a well-posedness result for the linear homogeneous case. Then, 
we turn to the inhomogeneous problem whose solution relies on a Bochner 
integral version of the variation of constants formula. In order to state this 
result we recall some basic aspects of measure theory. Finally, a nonlinear 
equation is investigated to illustrate the application of the Banach fixed point 
theorem to such problems. 

In ch. [6] we give a short overview of some known results concerning the 
Cauchy problem for wave maps. We state a local well-posedness theorem 
for general semilinear wave equations which is applicable to the wave maps 
system. Then, we mention some global results and give historical remarks. 

Ch. [7] is devoted to the study of self-similar solutions. We discuss blow 
up, criticality class and finite speed of propagation for the wave maps model. 
Then, we state a result concerning the existence of a countable set of self- 
similar solutions which provide explicit examples of blow up solutions. Fur- 
thermore, we give some numerical evidence for the universality of the blow up 
profile. We introduce a new coordinate system adapted to self-similarity and 
linearize the evolution equation. Using Sturm-Liouville theory we construct 
a self-adjoint operator which governs the linearized flow around a self-similar 
solution. Applying the previously developed operator techniques we prove 
well-posedness of the corresponding Cauchy problem. 

In ch. [8] we study the spectrum of the operator A which drives the 
linearized flow around the first self-similar solution. Using estimates for 
solutions around the singular endpoints of the eigenvalue equation and an 
oscillation argument we show that the point spectrum is empty. Furthermore, 
we identify a continuous spectrum. These results lead to the best possible 
growth estimate for solutions of the linearized equation around the first self- 
similar solution. Finally, we give an intuitive explanation why no further 
improvement of this result can be expected. 

In ch. [9] we introduce the functional calculus, a powerful method to define 
functions of self-adjoint operators. This yields an alternative approach to the 
study of abstract wave equations. We apply this method to the problem un- 
der investigation and reproduce the result previously obtained by semigroup 
techniques. 
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In ch. 10 we study the spectra of the operators A n which drive the 
linearized flow around the n-th self-similar wave map for n > 1. We show 
that A n has exactly n negative eigenvalues and derive a rough lower bound. 
Then, we numerically calculate the point spectra of A n and reveal a certain 
convergence of the eigenvalues. Finally, we construct a "limiting" operator 
Aqo. A numerical investigation strongly suggests that the spectrum of is 
the limit of the spectra of A n . 

In the last chapter we give an outlook, sketch some ideas how one could 
proceed and mention further results. 



Chapter 2 
Wave Maps 



In this chapter we define the model to be studied and note some of its basic 
properties. 

2.1 The Wave Maps Equation 

The wave maps equation is a system of semilinear wave equations where 
the involved nonlinearities have a geometric nature. We define an action 
functional for wave maps from Minkowski space to a Riemannian manifold. 

2.1.1 Basic Definitions and Notation 

We collect some basic notions from differential geometry required to define 
a wave map. The main purpose of this section is to fix notation in order 
to avoid confusion due to different conventions. For more information see 
e.g. [49], [12] or any other textbook on differential geometry. 

Smooth manifolds Let M be a set and d : M x M — > R a distance 
function on M, i.e. d(x,y) > 0, d(x,y) = <^ x = y, d(x,y) = d(y,x) and 
d(x,z) < d(x,y) + d(y,z) for all x,y,z G M. Then (M,d) is said to be a 
metric space. We define a manifold as a metric space M which is locally 
homeomorphic to IR m , i.e. for every x G M there exists a neighbourhood 
U C M of x, an integer m and a continuous bijective mapping : U — > M. m 
such that is continuous as well. If the number m is constant for all 
x G M (which we will always assume), the manifold M is said to be Tri- 
dimensional and we write dim(M) = m. An open subset U C M together 
with a homeomorphism u : U — > «([/) C M m is said to be a chart, denoted by 
([/, u). Two charts (£/, u) and {V, v) are said to be compatible if the mappings 
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uoy- 1 : v(U nV) ^ u(U n V) and v o u' 1 : u(C/ n V) -> u(t7 fl V) are 
smooth. A family ^4 = : z G /} of mutually compatible charts such 

that M = Uj e /C/j is called an at/as. Two atlases Al and A2 are said to be 
equivalent if any chart in A\ is compatible with any chart in A<i. A manifold 
together with an equivalence class of atlases is called a smooth manifold. 

Derivations, tangent space Let M, N be smooth manifolds. A continu- 
ous function / : M — > N is said to be smooth if for every x G M there exists 
a chart (U,u) on M with x G £7 and a chart (V, i>) on N with /(x) G V 
such that the mapping v o f o u -1 : u(C/ fl /^(V)) — ► u(V) is smooth. 
The set of smooth functions from M to N is denoted by C°°(M, N). A 
derivation £ x at x 6 M is a linear mapping £3. : C°° (M, R) — » R such that 
Ufg) = Uf)g(x) + f(x)Ug) for all f,g G C°°(M,R). It turns out that 
the set of derivations at x G M form an dim(M)-dimensional vector space. 
This vector space is referred to as the tangent space T X M of M at x. 

Tangent map, tangent bundle Let M, N be smooth manifolds, dim(M) = 
m and / G C°°(M, N). For x G M we define the tangent map at x T x f : 
T X M -> T /w iV by & i-> TJ ■ £, and T x f ■ £ x {g) := £,(<? ° /) for £ x G T X M 
and (7 G C°°(N, R). The tangent bundle TM is defined as the disjoint union 
of all tangent spaces, i.e. TM := {J xeM {x} x T S M. There is a canonical 
projection p : TM — > M given by (x,^) 1— > x. One defines the tangent map 
Tf : TM — > TiV as (x, £3) 1— > (f(x),T x f-£ x ). Given an atlas {(U^Ui) : i G /} 
on M there exists a canonical atlas on TM, namely {(p _1 (£7j), T?/j) : % G J}. 
This leads to the observation that TM is a 2m-dimensional smooth manifold. 

Vector fields, metric, pullback Let M be a smooth manifold. A vector 
field £ on M is a smooth mapping £ : M — > TM such that po( = idM- The 
set of smooth vector fields on M is denoted by X(M). A metric g on M is a 
mapping which associates to every x G M a scalar production p _1 ({x}) (or, 
equivalently, on T X M) such that the mapping x 1— > g(x)(^(x), C(%)) '■ M — > R 
is smooth for all £, £ G X(M). To simplify notation one usually writes g(£, () 
instead of <?(•)(£(•), C("))- ^ the scalar product is positive definite, g is said to 
be a Riemannian metric. If it is indefinite, g is called a pseudo-Riemannian 
metric. (M,g) is called a (pseudo-) Riemannian manifold. 

Given a (pseudo-) Riemannian manifold (N, g) and a smooth mapping 
/ G C°°(M, N) one can "pull back" the metric g on V to M. The pullback 

■"■By a scalar product on a real vector space X we mean a symmetric bilinear mapping 
b : X x X — ► M which is non-degenerate, i.e. b(x,y) — for all y G X is equivalent to 
x = 0. 
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metric f*g is defined by /*<?(*) CW) := g(f(x))(Tf ■ £(x),Tf ■ C(x)) 
for i 6 M and £, C G X(M). 



Cotangent bundle Given the tangent space T X M at x G M of a smooth 
manifold M with dim(M) = to, one denotes the dual space or cotangent space 
(in the sense of linear algebra) of T X M by T*M. T*M consists of all linear 
functional on T X M and is again an m-dimensional vector space. Using a 
similar construction as in the case of tangent spaces, the disjoint union of all 
cotangent spaces can be viewed as a 2m-dimensional manifold, the cotangent 
bundle T*M. 

Similarly to vector fields, one defines a one-form u as a smooth mapping 
u : M ■ T*M such that p o u = id M where p : T*M — > M denotes the 
canonical projection. 

Identification of TM with T*M Let (M,g) be a (pseudo-) Riemannian 
manifold. Then there is a natural way of identifying elements of T X M with 
elements of T*M: At each x G M the metric <? induces a scalar product ^ x on 
T X M. Let ^ G T X M and consider the mapping £ x h-> 5 r a .(^, •) : T X M -> T X *M. 
Due to the properties of a metric, this mapping is linear and bijective and 
therefore provides a natural way of identifying TM with T*M. 

2.1.2 Wave Maps on Minkowski Space 

Minkowski space Consider the set IR 4 . The identity id : IR 4 — > IR 4 yields 
a global chart on the metric space M := IR 4 and therefore M can be viewed 
as a 4-dimensional smooth manifold. Consider the i-th partial derivative for 
functions in C°°(M, IR) which induces a derivation at x G M, i.e. an element 
of T X M. We identify this element with the i-th unit vector in M. Therefore, 
the tangent space T X M at x G M is naturally isomorphic to M itself. We 
define a pseudo-Riemannian metric rj on M which induces a scalar product 
on T X M given by 



for v = (v°, v 1 , v 2 , v 3 ), w = (w°, w 1 , w 2 , w 3 ) G IR 4 . The pseudo-Riemannian 
manifold (M,rj) is called Minkowski space. 

The trace of the pullback metric Let (N, g) be a Riemannian manifold 
and (M, rj) be Minkowski space. For a smooth mapping / : M — > N consider 



3 




i=l 



14 



2.1. THE WAVE MAPS EQUATION 



the pullback metric f*g. At each i 6 M it induces a bilinear mapping 
(/*ff)s : T xM x T^M -> R. This defines a linear mapping £ x ^ (f*g) x ($, x , •) : 
T X M — > T*M = T X M, i.e. a linear operator on the 4-dimensional vector 
space T X M (we have implicitly identified T X M with T*M using 77). We denote 
the trace (in the sense of linear algebra) of this operator by tr v f*g(x). Then, 
tr v f*g is a smooth real-valued function on M. 

Variational formulation We define a functional S on the set C°°(M, N) 



for / e C°°(M, N). The integral is understood with respect to the ordinary 
Lebesgue measure on M 4 . We are interested in critical points of that func- 
tional, i.e. we intend to study compactly supported variations. However, for 
a function attaining values in a manifold it is not clear what it means to 
have compact support. Moreover, in the standard approach one considers an 
expression like S(f + etp) which makes absolutely no sense if the functions / 
and tp have values in a manifold. 

There are essentially two ways to go around this difficulty The first 
possibility is to restrict oneself to a single chart on N, i.e. one actually solves 
the problem for functions / having values in some R n . Clearly, this is a local 
notion and one does not obtain a global statement in this case. However, 
this approach is sufficient if the whole target manifold can be covered "almost 
completely" by a single chart. 

A second possibility is to embed N into some R n . Of course, in this case 
one has to ensure that all variations / + eip have values in N which restricts 
the set of admissible test functions ip. Hence, one has to make sure that the 
set of test functions is "large enough" (cf. [46]). 

Definition of wave maps For our purposes the first approach is sufficient, 
i.e. we restrict ourselves to a fixed single chart (V, v) on N . Via v _1 : v(V) — > 
V C N we pull back the metric g to v(V) C 1R" where n = dim(JV). The 
open set v(V) is a smooth manifold itself and (v ~ r )*g is a Riemannian metric 
on v(V). Hence, (v(V), (f" 1 )*^) is a Riemannian manifold. 



by 




f 



- (N, g) 




V 



(v(v),(v-rg) 
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Fix $ G C°°(M, v(V)) and consider $+ey? for e G R and smooth (p : M -> M n 
having compact support. In general the function $ + £(/? does not have values 
in v(V) but for e sufficiently close to it does. Having these considerations 
in mind we can define an action functional for wave maps. 

Definition 2.1.1. Let (M, 77) be Minkowski space, (N,g) an n-dimensional 
Riemannian manifold and (V, v) a chart on A. Define a functional S on 



for $ G C°°(M,v(V)). A function $ G C°°(M,v(V)) is said to be a critical 
point of the functional S if 



for all tp G C^°(M, W 1 ). A wave map on Minkowski space with target N is a 
mapping v _1 o $ : M — > A" such that $ is a critical point of S. 

Generalizations It is obvious that this definition can be extended to the 
case where M is a more general pseudo-Riemannian manifold. Such a gen- 
eralization is essential to study self-gravitating wave maps, i.e. to use wave 
maps as a matter model in general relativity. This has been done by various 
authors (e.g. [2], [50], [11]) although most of the work in this direction is 
confined to numerical studies or heuristic arguments due to the complexity 
of such systems. 

2.1.3 Local Expressions 

We intend to derive a local expression for the action functional which de- 
fines wave maps. As before, let (M, g) be Minkowski space, (A, g) an n- 
dimensional Riemannian manifold and (V, v) a chart on A. 

Lagrangian in local coordinates The /z-th partial derivative for func- 
tions defined on IR 4 yields a smooth vector field on M which is denoted 
by (/i = 0, ...,4). We set 77^ := T](d^d v ) and thus we have (77^) = 
diag(— 1, 1, 1, 1). Similarly, the A-tYi partial derivative on W 1 defines a smooth 
vector field 8a on v(V) C lR n (A — 1, . . . , n). Consider a smooth function 
$ G C°°(M,v(V)). The A-ih component of $ is denoted by $ A . By care- 
fully inserting the definitions one readily calculates T x <^-d^\ x = ^jr(x)dA\>s>(x) 
for x G M where Einstein's summation convention is assumed in the sequel. 



C°°(M,v(V)) by 




d 

de 



S($ + eip) 
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Thus, we have ^(y- 1 )* g(d^ d u )(x) = [{v^fg^x^T^-d^T^-d^) 
d$ A ( \ a® 1 

Q x v- \ X ) dx 1 - 



d * A {x)^{x)g AB {<b{x)) with g AB := ( V - l )*g(d A ,d B 



Vectors and one— forms We define the one-forms dx^ 1 by dx^(d u ) = 5^ u . 
Consider a vector field £ M <9 M . We assign to it the one-form u^dx^ where 
'■= t](j,v£, u ■ Since this is the usual multiplication of a matrix with a vector, 
the inverse operation is given by £ M = "q^Up where rf" are the components of 
the inverse of the matrix {f]^ v ). Thus, the identification of vectors with their 
duals takes place via multiplication with the matrix (rj^) and its inverse 

(iT). 

The wave maps system Taking the above considerations into account 
we arrive at 

tr^- 1 )*^) = ^(^$ a )(x)(^$ b )(x)(7 A b($(x)). 
Hence, the action functional is given by 

t v {d^ A ){d v ® B ){g AB o<S>). 



Following the standard approach (cf. [46], [19]) one interprets $> A and d^ A 
as independent variables and defines the Lagrangian £ by £($, <9$) := 
(dn$ A )(d u <& B )(g AB o $). Then, critical points (with respect to compactly 
supported variations) of S satisfy the Euler-Lagrange equations 

dC dC 

^d{d^ A ) d$ A 

A straight-forward calculation yields the equation 

□$ A + ^ v V A BC {^){d^ B ){d v ^ c ) = (2.1) 
where D$ A := r]^ u d^d u ^ A and 

r A Bc := ^9 AD {d B 9cD + d c g B D ~ d D g BC ) 
is a Christoffel symbol of the second kind associated to the metric g. Eq. 



(2.1) is known as the wave maps system. 

We remark that one can also study wave maps on Minkowski spaces 
R m+1 with different spatial dimension (m / 3). It is obvious how the above 
considerations are generalized to this case. 
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2.2 Wave Maps from Minkowski Space to the 
Three-Sphere 

Now we choose the three-sphere as a target manifold and derive an explicit 
expression for the action under certain symmetry assumptions. 

2.2.1 Basic Definitions 

The three-sphere The three-sphere S 3 C IR 4 is defined as 
S 3 := ^(x\x 2 ,x 3 ,x 4 ) e R 4 : X>*) 2 = X | " 

As a subset of the metric space IR 4 it is a metric space itself and via stere- 
ographic projections it is locally homeomorphic to IR 3 . Thus, S 3 is a 3- 
dimensional manifold. 

To define a chart (V,v) on S* 3 we give an explicit expression for v~ l . We 
define the mapping v" 1 : (0,7r) x (0, 7r) x (0,27r) — > S 3 by 

(-0,6,0) i — ^ (sin -0 sin O cos 0, sin -0 sin O sin 0, sin -0 cos O, cos ■0). 

The function v is a homeomorphism onto its image and therefore, (V, v) is a 
chart on S* 3 where V := im(t> -1 ). Similarly one can construct a whole atlas 
for S* 3 and therefore, S 3 is a smooth manifold. 

Let 5 be the standard metric on IR 4 which induces the Euclidean scalar 
product on every X^IR 4 = IR 4 . We consider the natural embedding % : S 3 — > 
IR 4 , i(x) = x for x G S 3 and pull back the metric 5 via i. The pullback metric 
g := i*5 is a Riemannian metric on S* 3 and therefore, (S 3 , g) is a Riemannian 
manifold. 

Co— rotational maps We formally define the usual spherical coordinates 
on Minkowski space (M, i]). Let u" 1 : IR x (0, oo) x (0, ir) x (0, 2tt) -> M be 
given by 

(t, r,9,<f) I— > (t, r sin ^ cos 99, r sin sin </?, r cos 6) . 

Then, -u -1 is a homeomorphism onto its image and therefore, (U, u) is a chart 
on M where U := im(-u _1 ). 

Consider the coordinate representation v o f o u~ 1 of a smooth mapping 
/ : M — > S* 3 . The function vo fou~ 4 assigns to a 4-tuple (t, r, </?) a 3-tuple 
(ip, O, 0). Therefore, -0, O and are functions of t, r, 9 and tp. In what follows 
we restrict ourselves to mappings / such that -0 is a function of t and r only 
and O = 9, = ip. Such maps are called co-rotational. 
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2.2.2 Explicit Local Expressions 

Coordinate representations We define $ := v o f and calculate the 
local expression h := (tr v <&* (v^ 1 )* g) o u^ 1 . However, in order to restrict 
ourselves to co-rotational maps it is desireable to work with the mapping 
$ o u^ 1 : u(U) — > v(V) rather than with $ itself. To do so we use the fact 
that the following diagram commutes. 

(M, V ) f - - (S\ g) * (IR 4 , 8) 

tr„$'>- 1 )*<? / u 

/ (u(U), (u-^—^iviV),^- 1 )*^ 




Hence, we read off the useful identity 

(tr^*^- 1 )*^) o u x = tr (tt -i).„($ o u- l Y{v- l y g . 

The metric on the three-sphere The i-th partial derivative (i = 1, 2, 3) 
for functions in C°°(v(V),M.) defines a smooth vector field on v(V) which 
will be denoted by <9j. We define g^ := (v~ 1 )*g(d i: dj), i,j = 1,2,3. A short 
calculation yields 

/ 1 

(gij&O, </>))= sin 2 ^ 

\ sin 2 i/j sin 2 





Minkowski metric in spherical coordinates Similarly, one defines the 
vector fields <9 M , \i = 0, . . . , 3 induced by the partial derivatives on u(U). 



The components r}^ v := (u 
coordinates are given by 



-i 



)*rj(dfj, : d v ) of the Minkowski metric in spherical 



( -1 








1 

















r sin' 



\ 



Calculation of the Lagrangian We make the ansatz ($ok r, 9, ip) = 
(ip(t,r),9,4>), i.e. we restrict ourselves to co-rotational maps. Calculating 
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the local expression h^ u := ($ o u l )*{v 1 )*g(d fl , d u ) for the pullback of the 
target metric we obtain 



( IPt AA 

AA A r o 

sin 2 if; 









\ 



\ sin 2 ip sin 2 © J 
Then h is given by h = ^(rj^hxv), explicitly 

2 

h(t, r) = -ipt(t, r) + ip 2 (t, r) + — sin 2 ip(t, r). 



The action functional Since the chart (U, u) covers the whole Minkowski 
space apart from a set of Lebesgue measure 0, integration over M is equiva- 
lent to integration over U . Applying the standard substitution rule we obtain 

S{$) = [ ti^*(v- 1 )*g= [ {tr ri <S>*(v- 1 )*g)ou- 1 \detJ u -i\ 

JU Ju(U) 

where J u -i denotes the Jacobian matrix of the function u~ l between the two 
open sets u(U) and U of 1R 4 . Thus, critical points of S can be obtained by 
calculating critical points of the functional 



if) i— > Att 



OO POO 



oo JO 



ijj 2 (t, r) - ip 2 (t, r) - — sin 2 ip(t, r) r 2 drdt (2.2) 



for if) : u(U) — > R and setting $ o u 1 (t,r,9,ip) := (if)(t,r),9,ip). Then, 
/ := v^ 1 o $ : M — > S* 3 is a co-rotational wave map. 



M^^S 3 



R 



Compactly supported variations of the functional (2.2) together with in- 
tegration by parts yields the equation 



2 , sin(2^) 

At ~ Ar A H 5 = 



(2.3) 



for critical points if). Thus, in the co-rotational case the wave map problem 
reduces to the single semilinear wave equation (2.3) for the function if). 
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2.3 Additional Remarks and References 

To be precise one distinguishes between wave maps on a Lorentzian manifold 
and the analogous mappings defined on a Riemannian manifold. The latter 
ones are usually called harmonic maps and were introduced by Fuller [21]. 
In particle physics wave maps first appeared in a paper by Gell-Mann and 
Levy [23] as nonlinear sigma models. For further applications in physics 
and historical remarks see the surveys [37], [38] and references therein. Mis- 
ner [37] points out that the wave maps equation contains nonlinearities in a 
very natural geometric way. Due to the lack of a vector space structure on 
the target manifold, the resulting field equation has to be nonlinear. There- 
fore, wave maps provide a rich source for nonlinear field theories. Moreover, 
they include well-known equations of mathematical physics as special cases. 
For instance, harmonic maps from Euclidean IR 3 to R satisfy the Laplace 
equation while wave maps from Minkowski space to R are solutions of the 
wave equation. The simplest nonlinear example is a harmonic map from an 
open interval (a, b) C R to a Riemannian manifold M. The resulting field 
equation is the geodesic equation on M. Therefore, the wave map functional 
provides an elegant and unified way for the derivation of many interesting 
equations (see [38] for more examples). 

For our purposes the most important feature is the geometric nature of 
the involved nonlinearities which classifies wave maps as "naturally nonlin- 
ear fields" [38]. Another example of such a "naturally nonlinear" field theory 
is general relativity and therefore, wave maps are promising candidates for 
modelling the more involved nonlinearities of Einstein's equations. Hence, 
the study of wave maps can provide insights into mechanisms which might be 
common to many field theories of that kind, in particular to general relativity 
which is the main physical motivation. From the mathematical point of view 



eq. (2.3 ) provides a simple example of a semilinear wave equation. The study 
of nonlinear wave equations is a large field which is rapidly developing and 
many difficult questions are still unanswered, even for examples as simple as 



eq. (2.3). Therefore, the wave map equation (2.3) provides an interesting 
object for a rigorous mathematical analysis. For a recent survey on open 
mathematical questions concerning general wave maps see [52]. An intro- 
duction to semilinear wave equations and wave maps as partial differential 
equations is given in [46]. 



Chapter 3 

Mathematical Preparation 



Before continuing the analysis of wave maps, we have to introduce the math- 
ematical machinery required later on. In this chapter we give an overview 
of function spaces, linear operators, spectral theory and semigroup theory. 
These tools are required for the treatment of evolution equations in an elegant 
abstract way. 

3.1 Function Spaces 

We define the standard Lebesgue and Sobolev spaces required for the analysis 
of partial differential equations. However, we restrict ourselves to the special 
case of spaces of functions defined on (open subsets of) the real line which is 
sufficient for our purposes. 

3.1.1 Banach Spaces 

Norms, completeness We start with the basic definition of a normed vec- 
tor space. Let X be a (not necessarily finite-dimensional) nontrivial vector 
space over the field K where K is either lor C. A norm on X is a mapping 
|| ■ || : X — > R that satisfies 

• > and \\x\\ = <^ x = 

• ||Ax|| = |A|||a;|| 

• + < \\x\\ + \\y\\ "triangle inequality" 

for all A G K and x,y G X. A vector space X together with a norm || • || on 
X (denoted by (X, || • ||)) is said to be a normed vector space. For brevity 
we will write X instead of (X, || ■ ||) if the norm we refer to follows from 
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the context. A normed vector space is a special case of a metric space (in 
the sense of point set topology) since d(x,y) := ||x — y\\ defines a distance 
on X. On a metric space there exists a distinguished topology defined by 
the distance function d which is called the metric topology. The open e balls 
B e (x) := {y G X : \\x — y\\ < e} form a basis of this topology. Therefore 
it follows by definition that x i— ► is a continuous mapping from X to R 
with respect to the metric topology. 

A Cauchy sequence in X is a sequence (xk) C X with the property that 
for every e > there exists a N G N such that — < e for all 
fc, j > N. If every Cauchy sequence in X has a limit in X, the normed 
vector space X is said to be complete. A complete normed vector space is 
called a Banach space. Not every normed vector space is a Banach space. 
Heuristically spoken it is possible that a Cauchy sequence "converges to a 
hole" . We give several elementary examples. 

Example 3.1.1. The vector space lR n together with the Euclidean norm ||a;|| 2 := 
J2k=i x 1 f° r x = • • • > x n) G M n is a Banach space which follows from the 
completeness of the real numbers. Similarly, every finite-dimensional normed 
vector space is complete. 

Example 3.1.2. The space C[a, b] of continuous functions from an inter- 
val [a, b] to C together with the norm ||«||c[a,6] : = su Pxe[a,6] 1^(^)1) u e 
C[a,b] is a Banach space. Similarly, the space C k [a,b] of fc-times contin- 
uously differentiable functions from [a, b] to C together with |H|c fc [a,f>] := 
2~^j=o su Pxe[a,b] \ u ^\ x )\i u ^ C k [a,b] is a Banach space. However, the space 
C 1 ^,^] equipped with the norm ||«||c[a,6]) u G C 1 ^,^] provides an example 
of a normed vector space which is not complete. 

Dense subsets, completion A subset A C X is said to be dense if for 
every x G X there exists a sequence (a^) C A such that x^ — > x in X. A 
Banach space is called separable if there exists a countable dense subset. 

Example 3.1.3. WeierstraB' approximation theorem states that polynomials 
are dense in C[a, b] and therefore, C[a, b] is separable since the set of all 
polynomials is countable. 

For every normed vector space (X, || ■ \\x) there exists a completion, i.e. 
a Banach space (X, || ■ ||) such that X is dense in X and \\x\\x = for all 
x <E X. A completion is unique up to isomorphisms, q 

1 To be more precise one requires the existence of a Banach space (Y,\\ ■ \\y) and a 
norm-preserving linear mapping i : X — » Y (an isometry) such that is dense in Y. 
Then Y is called the completion of X. In our definition we have implicitly identified X 
with i(X) which is possible since ||z(x)||y = \\x\\x for all x e X implies injectivity of i. 
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Example 3.1.4. Consider the space C^°(a,b) of infinitely often differentiable 
functions from (a, b) to C having compact support. For p > 1 we define 

/ nb \ l /P 

\\u\\Lv{a,b) ■= ( / \ll(x)\ P dx 



for u G C£°(a, b). It is a consequence of Minkowski's inequality that || • ||ip( a ,b) 
satisfies the triangle inequality and therefore, || • ||zp( a ,6) is a norm on C£°(a, b). 
The normed vector space (C^°(a, 6), || ■ ||ip( a ,6)) is n °t complete. Its completion 
is denoted by L p (a, b). 

Dual space, reflexivity A bounded linear functional on a normed vector 
space (A, || • \\x) is a linear mapping / : A — > K for which there exists a 
constant c > such that |/(x)| < c||a;||x for all x & X. We define the 
(topological) dual space X* of X as the set of all bounded linear functionals on 
X. One can define a norm on X* by := sup{|/(x)| : x & X, \\x\\x < !}■ 

It is easy to prove that (X*, \\ ■ \\ x *) is a Banach space (even if X is not 
complete). There exists a natural inclusion map % : X — > X** from X to the 
bidual space X** given by (i(x))(f) : = /(x) (x G X, f G A*). The linear 
map i is an isometry, i.e. it is norm-preserving. However, in general % is not 
surjective and therefore one cannot identify A with A**. A Banach space A 
which can be identified with its bidual space A** (i.e. the inclusion % is an 
isometric isomorphism) is said to be reflexive. 



Example 3.1.5. Consider the Banach space A := C[0, 1] (ex. 3.1.2). The 
linear mapping / : C[0, 1] — ► C given by f(u) := u(0) for u G C[0, 1] is 
bounded since ||w||c[o,i] < 1 implies |«(0)| < 1 and therefore ||/||x* < 1 
(recall the definition of || • \\x*)- Hence, / G A*. 

Example 3.1.6. Let A := C[0, 1] and define a norm on A by 

IMIx : = / \u(x)\dx 
Jo 

for m G A. Consider the mapping / : A — > K given by f(u) := tt(0). Then 
/ is linear but it is not bounded: It is easy to construct a sequence («&) of 
continuous functions with ||wfe||x = 1 for all k and itfe(O) — > oo for fc — + oo. 
Therefore, ||/||x* does not exist and / ^ A*. 

3.1.2 Hilbert Spaces 

A Hilbert space is a special case of a Banach space defined by the existence 
of an inner product. Many issues are dramatically simplified in that case and 
therefore Hilbert spaces deserve special attention. 
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Inner products Let if be a vector space over K. An inner product (-|-) is 
a mapping from H x H to IK that satisfies 

• (Ax + ny\z) = X(x\z) + fj>(y\z) 



• (Av) = (y\ x ) 

• (x\x) > if x 7^ 

for all x,y, z E H and A, \i E K. It is easily seen that \\x\\ := y/{x\x), x E H 
defines a norm on H. H is said to be a Hilbert space if (ii, || ■ ||) is a Banach 
space. 

Example 3.1.7. The vector space IR n together with the Euclidean inner prod- 
uct (x\y) := Y2=\ x kVk for x := (xi, . . . , x n ), y := (y 1 ,...,y n ) E R n is a 
(finite-dimensional) Hilbert space. 



Example 3.1.8. The space L 2 (a,b) (cf. ex. 3.1.4) equipped with the inner 
product 



(u\v) L 2 {a)b) := / u(x)v(x)dx 

J a 

is a prominent example of a Hilbert space. 



Dual space The characterization of the dual space of a Hilbert space is 
drastically simplified by the so-called Riesz representation theorem. 

Theorem 3.1.1. Let H be a Hilbert space overK. and f : H — > K a bounded 
linear functional on H. Then there exists a unique y E H such that f(x) = 
(x\y) for all x E H where (-|-) denotes the inner product on H . 

The proof is similar to the finite-dimensional case. Due to the Riesz 
representation theorem it is possible to identify the dual of a Hilbert space 
with the Hilbert space itself. In particular it follows that every Hilbert space 
is reflexive. 



3.1.3 Lebesgue— and Sobolev Spaces 

Now we switch from abstract Banach spaces to concrete function spaces. 
However, as mentioned in the beginning we will restrict ourselves to func- 
tions defined on one-dimensional real intervals. For the general definition 
of Sobolev spaces on subsets of lR n see e.g. [19]. For an extensive treatment 
of weighted Sobolev spaces as well as more general function spaces we refer 
to [55]. 
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Lebesgue spaces Let (a, b) C R be an open interval and p > 1. In ex. 



3.1.4 we have already defined the Lebesgue space L p {a^b) as the completion 



of compactly supported smooth functions with respect to the L p -norm. We 
note the important fact that in general a function u G L p (a, b) cannot be 
evaluated at a point, i.e. u(x) for x G (a, b) is not well-defined! The reason 
for this lies in the process of completion. By construction the completion 
is in fact a quotient space and therefore it consists of equivalence classes 
of functions rather than functions themselves. Two functions belong to the 
same class if they coincide almost everywhere, i.e. up to a set of Lebesgue 
measure 0. A single point x G (a, b) has Lebesgue measure and therefore, 
one can assign arbitrary values to u(x) without changing the function class. 
Thus, it makes no sense to speak of the value of a function u G L p (a, b) at a 
point x G (a, b). An important exception exists if a function class contains a 
continuous element. In this case one can define the value at a point by the 
evaluation of the unique continuous function in the class at that point. In 
the sequel we will do that implicitly. 

Additionally we remark that one often encounters statements like n Uj — > u 
in L p (a,b) implies Uj(x) — ► u(x) for almost all x G (a, b) n although, strictly 
speaking, this makes no sense. For convenience we will adopt this sloppy 
formulation having in mind what is actually meant by this, namely that 
Uj — > u in L p (a, b) implies that Vj(x) —>■ v(x) for almost all x G (a, b) and for 
all representatives Vj and v of the function classes Uj and u, respectively. 



Weighted Lebesgue spaces A weight function is a continuous, positive 
function w : (a, b) — > R. Define the weighted norm 

\\u\\ L p w (a,b) 

for u G C%°(a,b). The completion of C^°(a, b) with respect to the norm 
II " || is referred to as the weighted Lebesgue space L^(a, 6) with weight 
w. Again, the space L^(a,b) is a Hilbert space. 



\u(x)\ p w(x)dx 



Weak derivatives In order to apply functional analytic methods to partial 
differential equations it is necessary to generalize the notion of differentia- 
bility. The space L p oc (a, b) of locally Lebesgue integrable functions is defined 
by " 

Lf oc (a, b) := {u : (a, 6) -> C : u G L p (c, d) for all [c, d] C (a, b)}. 
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Inspired by the integration-by-parts-formula one says a function u G L p (a, b) 
is weakly differentiable if there exists a v G L\ oc (a,b) such that 

/ u(x)ip'(x)dx — — v(x)ip(x)dx 

J a J a 

for all ip G C^°(a, 6). In this case one sets u' :— v and u' is said to be the weak 
derivative of u. A weak derivative is unique (as a function class). As usual 
we adopt the notation for the j-th weak derivative and set := u. 

Sobolev spaces The Sobolev space W k ' p (a, b), k G N is defined to be 
the vector space of functions u G L p (a, b) such that weak derivatives 
for all j = l,...,k exist and belong to L p (a, b). The normed vector space 
(W k ' p (a,b), || 

' \\w k 'P(a,b)), where 

x i/p 

|«||w*.i>(a,6) : = ( Yl W u0) \\LP(a.h> 

is a Banach space. The special cases H k (a,b) := W k ' 2 (a, b) are of particular 
interest since they are Hilbert spaces with respect to the inner product 

k „b 



«|u)ff*(a,6) : = / U ij \x)v^)(x)dx. 

3=0 Ja 



Note the set inclusions C c °°(a,6) C IU fc ' p (a,6) C L p (a,b). It follows that 
W k ' p (a, b) is dense in L p (a, b) for all fceN. 

Weighted Sobolev spaces For a (k + l)-tuple u> = (u>o,tui, . . . , Wk) of 
weight functions we define the weighted Sobolev space W k ' p (a, b) as the vector 
space of functions u G L p w (a, b) such that weak derivatives for all j = 
1, . . . , k exist and G L^^a, b). The weighted Sobolev norm \\ ■ \[ w k,p^ ab ^ is 
defined by 

/ k \ Vp 

\i=o 

Again, ^W*' p (a, 6), || • || w fc .p( a 6)) is a Banach space. In particular, H k (a, 6) := 
H^;' p (a, b) is a Hilbert space with respect to the inner product 



u 



v) H ^a,b)-=Yl / W j {x)u ij) {x)v^)(x)dx. 

j=o Ja 
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3.2 Operator Theory 

In this section we give a brief review of well-known facts about basic operator 
theory. All of the material presented here can be found in standard textbooks, 
e.g. Kato's classic [26]. 



3.2.1 Linear Operators 

Linear operators, extensions, restrictions A linear operator A on a 
Banach space (X, \\-\\x) is a linear mapping from a subspace T>(A) C X (the 
domain) to X. We write A : T>(A) C X — > X to emphasize that topological 
properties of A are considered with respect to the topology on X and not 
to a possible different topology on T>(A). It is essential for the definition of 
a linear operator to specify the domain T>(A). Let A and B be two linear 
operators on a Banach space (X, || • \\ x )- A and B are considered to be equal 
if both their domains and their mapping rules coincide, i.e. A — B if and 
only if T>(A) = T>(B) and Ax = Bx for all x G ^(A). -B is said to be an 
extension (a restriction) of A if D(A) C T>(B) {T>(A) D T>(B)) and Ax = -Bx 
for all x G £>(A) C\V(B). 



Example 3.2.1. Let X := C[0, 1] (cf. ex. |3.1.2p , D(A) := C7 X [0,1] C X and 
Am := w', u G £>(A). Then A is a linear operator on X. Define T>(B) : = 
C 2 [0, 1], Bu := u' for u G V(B). Then 5 is a restriction of A. 



Sums and products The sum A+B of two linear operators on X is defined 
by V(A + B) := V(A) n V(B) and (A + B)x := Ax + Bx for x G £>(A + B). 
Note that A + B might be trivial, i.e. T>(A + B) = {0}. The product of A 
and B is defined by 

V(AB) := {x G 15(5) : Bx G P(A)} 

and ABx := A(5x) for x G T>(AB). Again, A5 might be trivial although A 
and B are not. 



3.2.2 Bounded and Closed Operators 

Bounded operators A linear operator A : T>(A) C X — > X is said to 
be bounded if P(A) = X and there exists a constant c > such that 
||Ax||x < clWU f° r a ll ^ G X. Since A is linear, boundedness is equiv- 
alent to continuity. The vector space of bounded linear operators on X is 
denoted by B(X). It is possible to define a norm on £>(X) by 

\\A\\b(x) ■= sup{||Ax|| x : x G X, \\x\\ x < 1}, 
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the so-called operator norm. It turns out that (B(X), || • ||s(x)) is a Banach 
space (provided that X is complete). 



Example 3.2.2. Let X = C[0, 1] (cf. ex. [3X2] ) and define V(A) : = X, 

(Au)(t) :- 



Then 

\Au(t)\ 



/ u(s)ds, u G X. 
Jo 



/ u(s)ds < / \u(s)\ds < / \u(s)\ds < \\u\\x 
Jo Jo Jo 



for all t G [0, 1] and therefore ||Au||x < \\u\\x for all «6l. Since A is linear 
this shows A G B(X). 

Closed operators Without requiring additional analytical properties, the 
study of unbounded linear operators reduces to linear algebra. An important 
property which is often present in applications is closedness which is a certain 
generalization of continuity. A linear operator A : T>(A) C X — > X on a 
Banach space (X, || • is said to be closed if for any sequence (x n ) C T>(A) 
with x n — > x and Ax n — > y it follows that x G T>(A) and Ax = y. If A is 
bounded then T>(A) = X and x n — > x already implies Ax n — > Ax. Hence, 
any bounded operator is closed. The converse, of course, is not true. 



Example 3.2.3. Consider again the Banach space X := C[0, 1] (cf. ex. 3.1.2). 
Define 

V{A) := ^ueX:t^ ^ EX 
and 

(Au)(t) :=^,i»6 2?(4 

Then, A is a linear operator on X. It is easy to show that A is closed but it 
is not bounded. 

The closed graph theorem states an important result concerning closed 
operators. 

Theorem 3.2.1. Let A : T>(A) C X — ► X be a closed linear operator on a 
Banach space X. IfV(A) = X then A is bounded. 

This theorem has many useful applications. If A : T>(A) C X —>■ X is 
bijective we can define the inverse operator A^ 1 : X —>■ X (A" 1 (Ax) = x 
for all x G T>(A) and ^(yl -1 ^) = x for all x G X). One easily shows that 
the inverse of a bijective closed operator is closed. Since ^(A^ 1 ) = X it 
follows by the closed graph theorem that A~ l is bounded. This observation 
is important for the spectral theory of closed operators. 



3.2. OPERATOR THEORY 



29 



Closeable operators, closure, graph norm A linear operator A : V(A) C 
X — > X on a Banach space (X, || • ||x) is said to be closeable if it admits a 
closed extension. If A is closeable then the smallest closed extension of A is 
called the closure of A and denoted by A. There exists an easy criterion to 
decide whether A is closeable or not. It turns out that A is closeable if and 
only if x n — > for x n G P(-A) and Ax n y imply y — 0. 

Let A be closed. A subspace F C T>(A) is said to be a core of A if there 
exists a closeable operator B : T>(B) Cl^I with T>(B) = Y and B = A. 

It is possible to define a new norm on the domain of A, the graph norm 
|| • || A given by 

||u|ft := \\Au\\ 2 x + |H|x, « e ^(^)- 

"D(A) equipped with || • ||a is a normed vector space which is complete if and 
only if A is closed. 

3.2.3 Self-Adjointness 

The adjoint Let A : V(A) C ff -> if and £ : £>(£) c if -> if be 
linear operators on a Hilbert space (H, (-\-)h)- B is said to be adjoint to A 
if (Ae|j/)h = (x\By) H for all a; G V{A) and y G £>(£). 

Example 3.2.4. Let if := L 2 (0, 1), := C^O, 1] and Au := v! for u G 

"D(A). Then A is a linear operator on the Hilbert space H. Define T>(B) : = 
{u G C^O, 1] : u(0) = u(l) = 0} and 5m := -v! for u G V{B). Then 5 is a 
linear operator on H which is adjoint to A (integration by parts). 

In general there are many operators which are adjoint to a given operator 
A. However, if A is densely defined (i.e. T>(A) is dense in H) then there 
exists a unique maximal operator A* adjoint to A (i.e. B C A* for all B 
adjoint to A). A* is called the adjoint and is constructed as follows. 

V(A*) := {y G H : 3z G H such that (Ax\y) H = (x\z) H for all x G V(A)} 

and whenever (Ae||/)# = (x\z)h holds for all x G P(-A) we write A*y = z. 
Since f^(A) is dense in H this construction yields a well-defined operator A* 
on if. It turns out that A* is always closed and therefore A** is a closed 
extension of A. 

Symmetric operators, self-adjointness A densely defined linear opera- 
tor A : T>(A) C ff — > ff on a Hilbert space (if, is said to be symmetric 
(or hermitian) if = (x|A|/)jy for all x,y £ 'D(A). It is said to be 
self-adjoint if A = A* (i.e. A is symmetric and f^(^4) = T>(A*)). Obvi- 
ously, self-adjoint operators are closed. A symmetric operator A is said to 
be essentially self-adjoint if its closure A is self-adjoint. 
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Example 3.2.5. Let H := L 2 (0, 1), V(A) := {u G C^O, 1] : «(0) = «(1) = 0}, 
Au := —ivl. Then A : T>(A) C H — > if is symmetric but it is not self- 
adjoint: Let f G C^O, 1] with t>(0) 7^ 0. Then we have —iv' G and 
(Au|u) H = (u\ - iu% for all u G This shows that v G but 

v <£ V(A) and hence, V(A) ^ V(A*). 

3.3 Spectral Theory 

We give an outline of the spectral theory of closed operators. Again, we 
refer to [26] for a detailed treatment of this subject. In what follows we 
assume A : T>(A) C X — > X to be a closed linear operator on a Banach space 
(X, || • || x ) over C. 

3.3.1 The Resolvent 

Resolvent set, resolvent For an operator on a finite-dimensional vector 
space the spectrum is defined as the set of all eigenvalues. In principle this 
definition could be carried over to the infinite-dimensional case without any 
modifications since the notions "eigenvalue" and "eigenvector" still make 
sense. However, if one intends to generalize results from finite-dimensional 
spectral theory it turns out that this definition is too restrictive. It is more 
fruitful to consider the resolvent instead. One defines the resolvent set p(A) 
of the operator A as the set of all A G C such that the inverse of the operator 
A — A exists and is bounded. The resolvent set is always open. For A G p(A) 
we define the resolvent -Ra(A) := (A — A)^ 1 . 

Properties of the resolvent Since A is closed it follows that A — A is 
closed as well. Therefore, (A — A)^ 1 : X — > X is automatically bounded (if it 
exists) because it is the inverse of a closed operator (closed graph theorem). 
This means that A G p{A) is equivalent to A— A being bijective. The resolvent 
Ra is an operator-valued function on p(A): It assigns to every complex 
number A G p(A) a bounded linear operator itU(A) : X — > X. In abstract 
notation this means that Ra is a mapping from p(A) to B(X). One can show 
that Ra ■ p{A) — > B(X) is piecewise holomorphic (piecewise because p(A) is 
not necessarily connected). Holomorphy in this context is defined as follows: 
Let / : U cC^ B(X) be a mapping from an open subset U of C to the 
Banach space of bounded operators on X. Then / is said to be complex 
differentiable at z G U if 



Um »/(*)- /(*b)l| B (*) 

z^zo Z — Zq 
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exists. / is said to be holomorphic on U if it is complex differentiable at all 
zeU. 



3.3.2 Spectra 

The spectrum The spectrum u{A) of A is defined as the complement of 
p(A) in C, i.e. cr(A) := C\p(A). Therefore, the spectrum is the set of all 
A G C such that the operator A — A fails to be bijective. In the finite- 
dimensional case the failure of A — A to be bijective is equivalent to A — A 
having a nontrivial kernel. In infinite dimensions there are more possibilities. 
Hence, there are different types of spectra. 

• The point spectrum a p (A): We say that A is in the point spectrum 
(7p(A) of A if there exists an x G V(A), x ^ such that Ax = Xx, 
i.e. A — A fails to be injective. Like in finite dimensions, x is called an 
eigenvector and elements of o~ p (A) are called eigenvalues. 

• The continuous spectrum a c (A): A G C belongs to the continuous 
spectrum o~ c (A) of A if A — A is injective but fails to be surjective and 
the range of A — A is dense in X. 

• The residual spectrum a r (A): A G C is said to be an element of the 
residual spectrum o~ r (A) of A if A— A is injective but fails to be surjective 
and the range of A — A is not dense in X. 

These three cases cover all possibilities. Therefore we have cr(A) = o~ p (A) U 
a c (A) U o~ r (A). Note that a p (A), cr c (A) and a r (A) are disjoint. In the finite- 
dimensional case we have o~(A) = a p (A). We remark that there is no overall 
standard for these definitions. Other classifications of spectra are possible 
and may be found in the literature. 



3.3.3 Structure of the Spectrum 

The spectra of closed operators In a finite-dimensional vector space 
over C the spectrum of a linear operator is bounded and nonempty due to 
the fundamental theorem of algebra. This is not true anymore when one 
considers unbounded operators in infinite-dimensional Banach spaces. We 
give two very simple examples to illustrate different spectral behaviour which 
one has to expect. 



Example 3.3.1. Let X := C[0, 1] (cf. ex. |3.1.2[ ), V(A) := C l [0, 1] and define 
Au := v! for u G T>(A). Then A : T>(A) C X — > X is an unbounded linear 
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operator on X. We show that A is closed. Consider a sequence (uk) G T>(A) 
with Uk — > u and Auk — > f ■ Note that 

lim Uk(t) = lim I itfc(O) + / u' k (s)ds 

Since Au^ = u' k —> f uniformly, we can interchange the limit with the integral 
sign and obtain 

u{t) = u(0) + / /(s)tfe. 

Hence, u G 12(^4) and Au = f which proves closedness of A. 

Now we calculate the spectrum of A. For every A G C the equation 
(A — A)u = has a nontrivial solution u(t) := e xt in P(A). This shows 
that <j(A) = <J P (A) = C and p(A) = 0. We remark that the term "point 
spectrum" might be misleading as this example shows. It is possible that the 
point spectrum consists of the whole complex plane and therefore it is not 
discrete in general. 

Example 3.3.2. Let X := C[0,1], V(A) := {u G C^O, 1] : w(0) = 0} and 
Au = u' for u G T>(A). Then A : T>(A) C X — ► X is an unbounded linear 



operator on X. A similar reasoning as in ex. 3.3.1 shows that A is closed. 
To calculate the point spectrum of A we consider the equation (A — A)u = 0. 
Formally, this equation has the solution u(t) := ce xt , c G C but u £ T>(A) 
as long as c ^ 0. Therefore, the only solution which is in T>(A) is the trivial 
one. This shows that ct p (A) = 0. Next we show surjectivity of A — A. Let 
/ G X and consider the equation (A — A)u = f. A formal solution is given 
by 

u(t) = -e xt [ e~ Xs f(s)ds 



o 



as is easily checked by direct computation. The so-defined u is continuously 
differentiable, satisfies u(0) = and hence, is an element of T>(A). Therefore, 
A — A is surjective for all A G C which, together with <7 P (A) = 0, shows 
a{A) = and p{A) = C. The resolvent R A (X) : X -»• X is given by 



(R A (X)f)(t) = -e xt f e- Xs f(s)ds. 
Jo 



These two examples show that the spectrum of a closed operator can fill 
the whole complex plane or it may be empty. Roughly speaking, everything in 
between is also possible. Therefore, the spectra of general closed operators 
can be very complicated. However, for many differential operators arising 
from mathematical physics one rediscovers the familiar situation of having a 
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pure point spectrum which is discrete. The analysis is additionally simplified 
if the problems can be formulated in a self-adjoint way. In this case one has 
a good knowledge of the structure of the spectrum. 



The Spectrum of a self-adjoint operator Like a symmetric operator 
in finite dimensions, a self-adjoint operator can only have real eigenvalues. 
This statement is already true for symmetric operators and the proof is the 
same as in finite dimensions. However, in the case of a self-adjoint operator 
a much stronger result holds. It turns out that the whole spectrum has to 
be a subset of the real axis and not only the point spectrum. Furthermore, 
self-adjoint operators do not have residual spectra as the following argument 
shows. Consider a self-adjoint operator A : T>(A) C H — > H on a Hilbert 
space H. Suppose A G a r (A). It follows that A is real and by definition 
of a r (A), A is not an eigenvalue and the range of A — A is not dense in H . 
Hence, there exists a y G H, y ^ such that ((A — A)x\y)u = for all 
x G V(A). This implies (Ax\y)n = (x\\y) H for all x G V(A) which shows 
that y G V(A*) = V(A). Therefore, we have (x\Ay) H = (x\Xy) H for all 
x G T>(A) and, since T>(A) is dense in H, we conclude that Ay = Xy which 
shows that A is an eigenvalue of A, a contradiction. 

Moreover, due to their importance for quantum mechanics, self-adjoint 
operators have been investigated systematically for a long time and there 
exists a well-developed general theory (cf. [41], [42]). For instance, there 
is a spectral theorem and it is possible to define functions of self-adjoint 
operators with the help of the so-called functional calculus. These properties 
simplify the analysis tremendously and therefore it is desireable to formulate 
the problem to be studied in a self-adjoint way whenever this is possible. 



3.4 Semigroup Theory 

A thorough introduction to the theory of semigroups of linear operators can 
be found in [17]. We also mention the classic [40]. For a treatment of semi- 
group theory with a focus on physical applications (in particular relativity 
and astrophysics) we refer to [5]. 



3.4.1 Motivation 

Let X be a Banach space, A : V(A) C X — > X a linear operator and 
u : [0, oo) — > X a function. We are interested in solving the evolution 
equation 

ju{t)=Au{t) (3.1) 
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with initial data u(0) = Uq. Formally, the solution is given by u(t) = 
exp(tA)uQ. Therefore, our aim is to give precise meaning to the exponential 
of a linear operator. There is no problem if the operator A is bounded. In 
this case the exponential series converges absolutely and therefore, exp(tA) 
can be defined as usual via a power series. This situation is completely anal- 
ogous to the case of a linear operator on a finite-dimensional vector space. 
However, as soon as the operator A is unbounded, serious problems arise. Un- 
fortunately, this is the most important case for applications. For instance, if 
we intend to write a partial differential equation in the form of eq. (3.1 ) then 
A is a differential operator and hence, it is unbounded. Semigroup theory 
provides a framework for tackling these problems and yields a method for 
the treatment of ordinary differential equations on Banach spaces. 



3.4.2 Generation of Semigroups 

Strongly continuous semigroups Let X be a Banach space and S a 
family of bounded operators on X depending on a real nonnegative param- 
eter, i.e. S(t) G B(X) for t > 0. S is called a one-parameter semigroup if 
S(0) = idx and S(t + s) — S(t)S(s) for all t, s > 0. The semigroup S is 
said to be strongly continuous if the mapping t i— > S(t)x : [0, oo) — > X is 
continuous for all i6l. We remark that this is a weaker notion of continu- 
ity than the more natural-looking requirement of S : [0, oo) — > B(X) being 
continuous. The latter is referred to as uniform continuity and will not play 
a role in our further considerations since it is too restrictive. 



Generators Let S : [0, oo) — > B(X) be a strongly continuous one-parameter 
semigroup of linear operators on a Banach space X. We define 

V(A) := jx G X : lim - % exists j . 

For x G T>(A) define Ax := lim t ^o+ j(S(t)x — x). One easily shows that 
A : T>(A) C X —>■ X is a linear operator which is called the generator of 
the semigroup S. It turns out that the generator determines the semigroup 
uniquely, i.e. different semigroups have different generators. Another useful 
property of semigroups is the fact that x G T>(A) implies S(t)x G T>(A) for 
all t > which will be important for the application to partial differential 
equations. 



Generation results An obviously important problem is to determine whether 
a given linear operator generates a semigroup or not. It turns out that gen- 
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erators can be completely characterized by spectral properties. This result 
is known as the Hille- Yosida theorem. 

Theorem 3.4.1 (Hille- Yosida). Let X be a Banach space and A : V(A) C 
X — > X a linear operator on X. Then, A generates a strongly continuous 
one-parameter semigroup S on X satisfying ^S{t)^B{x) < 1 f or a ^ t > 0, if 
and only if 

• T>(A) is dense in X 

• A is closed 

• (0,oo) cp{A) 

• \\R A (X)\\B(x)<jforallX>0. 

A semigroup S : [0, oo) — ► B(X) which satisfies < 1 for all 

t > is called a contraction semigroup. We remark that semigroups which 
satisfy a more general growth bound given by ||5'(£)||b(.x') < e wt foruEM. 
can be reduced to the Hille- Yosida case by rescaling, i.e. considering the 
semigroup t i— > e~ w *S'(t) instead of S. For even more general bounds given 
by HS^IIepO < Me w ', M > 0, there is a similar generation theorem but one 
has to take powers of the resolvent into account (cf. [17]). 



3.4.3 The Abstract Cauchy Problem 

Statement of the problem Let X be a Banach space and A : T>(A) C 
X — > X a linear operator. We consider the abstract evolution problem 

f t u(t) = Au(t) for t > (3 2) 

u(0) =u 

for a function u : [0, oo) — > X and initial data Uq G X. 

Definition 3.4.1. A function u : [0, oo) — > T>(A) C X is called a classical 



solution of eq. (3.2) if it is continuously differentiable and satisfies eq. (3.2) 



(in particular we must have u$ G V(A)). 



Well— posedness In order to obtain a reasonable Cauchy problem we re- 
quire eq. (3.2) to possess the following property: For given initial data there 
exists a unique solution which depends continuously on the data. This as- 
sumption is usually referred to as well-posedness. Now we give the precise 
definition. 
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Definition 3.4.2. The abstract evolution problem (3.2) is said to be well 
posed if for every u G T>(A) there exists a unique classical solution u 
[0, oo) -> V(A) of eq. Q 
u ^ : V(A) C X -»• 

intervals [0, T], T > 0. 



such that for each fixed t > the mapping 
X is uniformly continuous in £ on compact 



We remark that there is no overall standard and different definitions of 
well-posedness might be found in the literature, e.g. one sometimes includes 
growth bounds for the solution. In this sense our requirements are weaker. 



Well— posedness and semigroups Suppose A generates a strongly con- 
tinuous semigroup S on X satisfying [^(t) < C(t) where C is a contin- 
uous positive function. Let uq G T)(A). Then the unique classical solution u 

S(t)u (recall that S(t)u G V(A) for all t > 

\\S{t)u \\ x < C{t)\\u \\x < 
u(t) : V(A) Cl^I 
T > 0. Hence, if A generates 



of eq. (3.2) is given by u(t) 
if u 



G V(A)). Moreover, we have ||«(t;i|x 
sup tg j 0T ] C(t)||uo||x and therefore, the mapping u 
is uniformly continuous in t on each [0,T], 



a strongly continuous semigroup S satisfying HS^Hape) < C(t) then the 
abstract Cauchy problem eq. (3.2) is well-posed. 



Furthermore, the existence of a semigroup extends the notion of a solution 
of eq. (3.2). Since S(t)uo is defined for general uo 
u Q G T>(A) we can define (generalized) solutions u 



G X and not only for 

[0,oo) 



X of eq. (3.2) 



by u(t) := S(t)uo for u G X. 



3.4.4 Second Order Cauchy Problems 

Reduction to first order system We are concerned with wave equations 
and hence it is desireable to have a semigroup formulation for second order 
Cauchy problems. This can be achieved by the usual reduction of a sec- 
ond order equation to a first order system. Consider informally the Cauchy 
problem 

f ^ u (t) = Au{t) for t > 

\ u(0) = U O ,U t (0) = Ml 

This second order equation can be rewritten as the first order system 

d ( u{t) \ = ( 1 \ / u(t) \ 
dt V ut(t) J \A ) \ ut(t) J 



with initial data (u(0),u t (0)) = (uq,u\). 
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Function spaces We have avoided a precise consideration of the involved 
function spaces because this issue can be very subtle. Instead we will prove 
a generation result for the case when A satisfies certain additional conditions 
which will be enough for our purposes. However, higher order Cauchy prob- 
lems on Banach spaces have been systematically investigated and we refer 
to [20] and [58] for more information. 
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Chapter 4 

Self— Adjoint Operators 



We continue collecting mathematical prerequisites. In the first part of this 
chapter we note some important properties of self-adjoint operators. Then 
we turn to Sturm-Liouville theory which deals with self-adjoint operators 
that are generated by symmetric ordinary differential expressions of second 
order. 

4.1 Properties of Self-Adjoint Operators 

We recall some well-known properties of general self-adjoint operators. Ev- 
erything in this section can be found in standard textbooks, e.g. [26], [59]. 

4.1.1 The Square Root 

Let A : V(A) C H — > H be a self-adjoint operator on a Hilbert space H. 
Then, for u G T>(A), we have (Au\u)h = {u\Au)h = (Au\u)h and hence, 
(Au\u)h G K for all u G T>(A). A is said to be nonnegative if (Au\u)h > 
for all u G T>(A). For a nonnegative operator A there exists a square root 
A 1/2 with nice properties. 

Theorem 4.1.1. Let H be a Hilbert space and A : V(A) C H — > H a 
self-adjoint operator which satisfies (Au\u)h > for all u G T>(A). Then, 
there exists a unique self-adjoint operator A 1 / 2 such that V(A) is a core of 
A 1 ' 2 , A = (A 1 / 2 ) 2 and (A 1/2 u\u) H > for all u G V(A l/2 ). Moreover, 
BA 1/2 C A 1/2 B whenever B A C AB for B G B(H), i.e. A 1/2 commutes with 
any bounded operator that commutes with A. 



Proof. The claim follows from [26], p. 281 Theorem 3.35 together with [26], 
p. 279, Problem 3.32. □ 
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In connection with the square root we note the following triviality which 
will be useful later on. 

Lemma 4.1.1. Let A : T>(A) C H —>■ H be a nonnegative self-adjoint 
operator on a Hilbert space H. If (Au\u)h > j(u\u)h for some 7 > and 
all u E V(A) then (A 1/2 u\A 1/2 u) H > j{u\u) H for all u E V{A l/2 ). 

Proof. From (Au\u)h > j(u\u)h we conclude that (A 1 / 2 u\A 1 I 2 u)h > 7 (u\u)h 
for all u E V{A). Let u E V(A 1/2 ). Since V(A) is a core for A 1 / 2 by Theo- 
rem 4.1.1, there exists a sequence (uj) C T>(A) such that Uj — > u in H and 
A 1/2 Uj A 1/2 u in H. We have (A 1/2 u j \A 1/2 u j ) H > ^{uj\uj) H for all j EN 
and hence, this inequality remains valid in the limit j ^ 00. □ 

Remark 4.1.1. A self-adjoint operator A satisfying (Au\u)h > i(u\u)h for 
all u E T>(A) and some 7 E 1R is called semibounded from below. 



4.1.2 Boundedness of the Spectrum 

For self-adjoint operators there exists an important connection between semi- 
boundedness and boundedness of the spectrum similar to the finite dimen- 
sional case. The following result makes it possible to deduce inequalities by 
studying spectra. 

Proposition 4.1.1. Let A : T>(A) C H — > H be a self-adjoint operator on 
a Hilbert space H and 7GI. Then, (Au\u)h > l(u\u)u for all u E T>(A) if 
and only ifmia(A) > 7. 

Proof. See [26], p. 278. □ 



4.2 Sturm— Liouville Operators 

In this section we review some aspects of Sturm-Liouville theory. We will 
also give proofs for most of the results stated below since they are very 
instructive. For more detailed expositions see e.g. [39], [54]. 

4.2.1 Basic Definitions 

Absolutely continuous functions We introduce a new function space 
which turns out to be useful in connection with ordinary differential oper- 
ators. Let u : [a, b] — > C. Then, u E AC[a, b] if there exists a function 
v E L 1 (a, b) and a c E (a, b) such that u(x)—u(c) = J x v(s)ds for all x E [a, b}. 
In particular, it follows that u is continuous and possesses a weak derivative 
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given by v. Clearly, AC [a, b] is a vector space and it is called the space 
of absolutely continuous functions. As usual, the local version AC\ oc (a, b) is 
defined as AC\ oc (a, b) := {u G AC[c, d] : [c,d] C (a, b)}. 

Sturm— Liouville operators A Sturm-Liouville operator is generated by 
a formal differential expression a of the form 

era := — (—(pu'Y + qu) 
w 

where p, q : (a, b) — > R, l/p,q £ L\ oc (a, b) and w is a weight function. 
Formal integration by parts yields the Green's formula 



(au)(x)v(x)w(x)dx = [u,v] p (d) — [u,v] p (c) + j u(x)(av)(x)w(x)dx 

° (4-1) 

for any [c,d] C (a, b) where [u, v] p (x) := u(x)(pv')(x) — (pu')(x)v(x). 

The initial value problem Consider the equation 

— {pu')' + qu = Xwu 

on (a, b) C R for a complex parameter A. The assumptions on p, g are 
sufficient to guarantee existence and uniqueness of the initial value problem. 
More precise, we have the following theorem. 

Theorem 4.2.1. Let c G (a, b) and 6 C. Then, there exists a unique 
function u(-, A) G AC\ oc (a, b) with pu'(-, A) G AC\ oc (a, b) such that —(pu')' + 
qu = Xwu and u(c, A) = £ ; (pu')(c, A) = r\. Moreover, the function u(x, •) is 
holomorphic on C for any x G (a, 6). 

The proof is similar to the classic Picard-Lindelof theorem and can be 
found e.g. in [39]. 

The maximal operator We set H := L^(a, b) and define the maximal 
operator A± : T>(Ai) C H — > H associated to a by 

V(Ai) := {u G H : G AC loc (a, 6), aw G H} 

and Ai« := aw for u G Note that functions u in satisfy the 

minimal requirements to give sense to era and to guarantee u G H as well as 
era G H. That is why Ai is called the maximal operator. Note further that 
u, v G X'(Ai) implies the existence of 

[u, v] p (a) := lim [u, v] p (x) and [w,vL(6) := lim [u, v] p (x) 



as follows from the Green's formula eq. (4.1). 
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Endpoint classification The endpoint a of the interval (a, b) is classified 
as follows. 

• The endpoint a is said to be regular if a > — oo and there exists a 
c G (a, b) such that 1/p, q,w G L 1 (a, c). 

• a is said to be in the limit-circle case if there exist u,v G ^(^i) such 
that [w, f ] p (a) 7^ 0. 

• Finally, a is said to be in the limit-point case if [w,u] p (a) = for all 
m,d6 V(Ai). 

An analogous classification is applied to b. 

There exists a close connection between endpoint classification and inte- 
grability of solutions of the equation otu = known as the Weyl alternative. 
This fact will be discussed later on. 



The minimal operator The minimal operator A : T>(Aq) C H — > H 
associated to a is defined by 

V(A ) : = {u G V(Ai) : [u,v] p (a) = [u,v] p (b) = for all v G V{A 1 )} 

and A u := au for u G V(A Q ). Obviously, we have A C A\. Furthermore, 
(A u\v)h = (u\Aiv)h for all u G T>(A ) and v G i.e. the operators 

Aq and Ai are adjoint to each other. This observation leads to the first easy 
relationship between A and A\. 

Lemma 4.2.1. We have the inclusion uxyAq C (ker Ai) 1 - . 

Proof. Let / G imAo, i.e. there exists a u G P(>lo) such that Aqu = f. 
Choose any v G ker A\. Then, since Aq and A\ are adjoint to each other, we 
have (f\v) H = (A u\v) H = {u\A x v) H = 0. □ 



4.2.2 Regular Sturm— Liouville Operators 

We prove some properties of Aq and A±. However, it is easier to consider the 
regular case first. Thus, we assume that both endpoints a and b are regular, 
i.e. l/p,q,w G L 1 (a, b). The regular case allows major simplifications. First 
of all, any u G C[a, b] belongs to H = L^(a, 6) thanks to the inequality 

J \u(x)\ 2 w(x)dx < || |M| 2 ||(7[a,()] IMU 1 ^^) < 00 anc ^' secon dly, the initial value 
problem can be uniquely solved at the endpoints, i.e. the point c in Theorem 



4.2.1 can be chosen to be a or b. Moreover, any solution u of A\u = satisfies 



u,pu' G AC[a, b]. 
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Relation between imA and ker Ai 



Lemma 4.2.2. Let f G H. Then, we can find a function u G T>(Ai) satis- 
fying u(a) = (pu')(a) = and A\U = f . In particular, A\ is surjective. 

Proof. Let / G H. Since A\ is regular we can find two linearly indepen- 
dent functions Ui,u 2 G T>(Ai) satisfying A\Uj = and Uj,pu'j G AC [a, 6] 
for j = 1,2. Moreover, [u 1 ,u 2 ] p (x) is a constant ^ as follows from the 
Green's formula eq. (4.1) and thus, by normalization we can assume that 
[ui,u 2 ] p (x) = 1. We set 

PX PX 

u(x) :— U\{x) j U2(s) f '(s)w(s)ds — u 2 (x) / ui(s)f(s)w(s)ds. 

J a J a 

The involved integrals exist thanks to the inequality 



Uj(x) f (x)w(x)dx < 



b \ 1/2 

\uj(x)\ 2 w(x)dx 



H 



< 



\Uj\ 2 \\c[a,b]\\ W \\L 1 (a,b)) 



1/2 



H < OO. 



A direct computation reveals u,pu' G AC[a,b], u(a) = (pu')(a) 

Am = f. 



and 

□ 



Now we are ready to prove the converse statement to Lemma 4.2.1 
Lemma 4.2.3. Let f G (ker A^ 1 . Then, f G imA . 



Proof. Let / G (ker Ai)- 1 . According to Lemma 4.2.2 we can find a u G T>(Ai) 
such that u(a) = (pu')(a) = and A%u = f. Let Ui,u 2 G ker Ax satisfy 
u\(b) = {pu' 2 )(b) = and u 2 {b) = {pu'^ib) = 1. Such functions exist thanks 
to Theorem |4.2.1| and the regularity of A\. Invoking Green's formula eq. 
(4.1 ) we calculate 

= (f\uj) H = (Axu\uj) B = [u,Uj) p (b) - [u,Uj] p (a) + {^Axu^h = [u,Uj) p (b) 

for j = 1, 2 which yields u(b) = (pu')(b) = 0. Thus, u G T>(Aq) and we have 
Aqu = A\u = f which shows that / G imAo- □ 



Remark 4.2.1. Taking together Lemmas |4.2.1| and |4.2.3| we have proved that 
inn4o = (kerAi" 



4.2.1 



lm- 



Note that the existence and uniqueness Theorem 
plies that ker A\ is two-dimensional and in particular, as a finite-dimensional 
normed vector space, it is closed. Thus, by taking the orthogonal comple- 
ment, the relation imv4 = (ker Ai) 1 - implies (imAo) -1 = ker A]_. 
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Density of T>(Aq) Next we prove that T>(Aq) is dense in H which shows 
that A and A\ are densely defined. 

Lemma 4.2.4. The domain T>(Aq) is dense in H . 

Proof. It suffices to show that any element which is orthogonal to T>(Aq) is 
zero. Let / G T^Ao) -1 , i.e. (u\f)n = for all u G T>(Aq). Invoking Lemma 



4.2.2 we see that there exists a v G D(A\) such that A\V = f. Thus, we 



have = {u\f)n = (u\Aiv)h = (A u\v)h for all u G T>(A ) which shows that 



v G (imAo) . However, according to Remark 4.2.1 we have (imAo) = ker A\ 



and thus, v G ker A% which yields = A\v = f. □ 

The adjoints We calculate the adjoints of Aq and A\ (which are now 
known to exist since Aq and Ax are densely defined). 

Lemma 4.2.5. The operator A\ is the adjoint of Aq, i.e. Aq = A\. 



Proof. According to the Green's formula eq. (4.1) we know that (Aqu\v)h = 
(u\Aiv)h for all u G T>(A ) and v G T>{A\) which means that A\ C Aq. Let 
v G T>(Aq), i.e. there exists an / G H such that {u\f)u = (^o^l^)^ f° r a h 



u G T>(Aq). However, since Ai is surjective by Lemma |4.2.2 there exists a 
v G T>(Ai) such that A\v = f. Hence, we have 



(Aqu\v) h = (w|/)h = {u\A x v)h = (Aqu\v] 



H 



which shows that (A u\v — v)h = 0. We conclude that v — v G (imAo) 1 - = 
ker A\ C T>(Ai) (Remark |4.2.1 ) and therefore, since v G D(Ai), we infer that 



v G T>(A\). Thus, we have shown that Aq C A\ which finishes the proof. □ 

Remark 4.2.2. In particular it follows that A\ is closed since it coincides with 
the adjoint of a densely defined operator which is always closed. 

Lemma 4.2.6. The operator Aq is the adjoint of A\, i.e. A\ = Aq. 



Proof. Applying "*" to the relation A\ = A* Q (Lemma 4.2.5) yields A\ = 
Aq* D Aq. Thus, it remains to show that A\ C Aq. The relation A% D Aq 
implies A* C A* and we infer A\ C A x . Let v G V(A\) C By 
definition of the adjoint we have (A\u\v)h = (u\A*v)h = (u\A\v)h for all 
u G P(Ai). On the other hand, by the Green's formula eq. ( |4.1 ), we have 



(A\u\v)h = [u, v] p (b) — [u,v] p (a) + (u\Aiv)h for all u G T>(A\) and thus, 
[u, v] p (b) — [u,v] p (a) = for all u G V(A 1 ). However, since A\ is regular 
we can choose u(a), (pu')(a), u(b), (pu')(b) arbitrarily and we conclude that 
v(a) = (pv')(a) = v(b) = (pv')(b) = which shows that v G T>(Aq) and 
therefore, A\ C A . □ 

Remark 4.2.3. Again, it follows that A is closed. 
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A special case The operator A is symmetric but it is not self-ajoint since 
Aq = Ai 7^ A . Thus, possible self-adjoint extensions of A lie between A 
and Aq = A\. We give a simple example. Let T>(A) := {u G T>(Ai) : u(a) = 



u(b) = 0} and Au 
A is self-adjoint. 



au for u eV(A). Then, A G A G A 1 . We claim that 



Lemma 4.2.7. The operator A is self-adjoint. 

Proof. First of all we note that (Au\v)h = {u\Av)h for all u, v G by 
the Green's formula eq. (4.1) and therefore, since A is densely defined by 
Lemma 4.2.4 A is symmetric. Hence, it remains to show that A* G A. 

The relation Aq G A G Ai together with Lemmas |4.2.5 and 4.2.6 imply 
A 1 = A* d A* D A\ = A . Let v G V(A*) C X>(Ai). By definition we 

(m|Ai«)h for all u e V(A). On the other 



have (Au\v)h = (u\A*v)h 
hand, according to the Green's formula eq. (4.1), we observe that (Au\v)h = 
[u,v] p (b) — [u,v] p (a) + {u\A\v)h for all u G T>{A) which yields [u, v] p (b) — 
[u,v] p (a) = for all u G D(A). However, since A is regular we can choose 
(pu')(a) and (pu')(b) arbitrarily which yields v(a) = v(b) = 0. Thus, v G 
T>(A) and therefore we have A* G A. □ 



4.2.3 Singular Sturm— Liouville Operators 

The results of the previous section are not very useful for concrete appli- 
cations since most of the interesting Sturm-Liouville problems are singular. 
Thus, we have to generalize the theory to the singular case, i.e. in the sequel 
we merely assume l/p,q,w G L\ oc {a, b). 

An auxiliary operator We define an auxiliary operator A : T>(Aq) G 
H -> H by 

T>(A ) := {u G T>(A\) : u has compact support} 
and A u := au for u G V(A ). 

Regularized operators The main idea is to restrict the problem to a 
fixed interval I := [c, d] G (a, b) where everything is regular and then use the 
fact that I is arbitrary. We define H 1 := L^(c, d). Consider the mapping 
i : H 1 — > H, u \— ► i(u) defined by 



{i(u)){x) :-- 



u(x) if x G [c, d] 
if x G (a,b)\[c,d] 
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Since i : H 1 — > H is injective we can identify u with i(u) and in this sense we 
have the inclusion H 1 C if. We define the regular operators A ! and /!{ on 
if 7 analogous to A and Ai where (a, 6) is substituted by (c, <f). According 
to the results of the previous section we have {Aq)* = A[ and {A[)* = Aq. 
Note further that u G £>(v4q) implies i{u) G £>(Ad) and hence, we have 
V(Aq) C T>(Aq) in the same sense as H 1 C i/. 



Density of v4o We show that £>(Ao) is dense in H. 
Lemma 4.2.8. The operator Aq is densely defined. 



Proof. Let f E H with (u\f) H = for all u G V(A Q ). Since P(A 7 ) C X>(i ) 
it follows that (u\f)n = for all u G V^Aq) and we can substitute the inner 
product on H by the inner product on H 1 which yields (u\f\i) H i = f or all 
u G V(A I ). Since A^ : X>(A^) C H 1 ->• if J is densely defined by Lemma [42^4 



we conclude that /|j = 0. However, J was arbitrary and therefore we infer 
that / = almost everywhere. □ 



Lemma 4.2.8 together with (A u\v)h = {u\A v)h = for all u,d£ ^(^4o) 
imply that Ao is symmetric. Furthermore, since Aq G Aq G A\, we observe 
that A and Ax are densely defined as well. 



The equation au = f We consider the inhomogeneous equation au = f. 

Lemma 4.2.9. Let f G H. Then, there exists a u G AC\ oc (a,b) with pu' G 
AC\ oc (a,b) satisfying au = f. 



Proof. According to Theorem 4.2.1 there exist two linearly independent func- 



tions Ui,U2 G AC\ oc (a, b) with pu'j G AC\ oc (a } b) satisfying auj = (j = 1, 2). 
Without loss of generality we can assume [ui,U2] p (x) = 1 (Green's formula 



eq. (4.1)). We define u by 



u(x) := cxUi(x) + c 2 u 2 {x) + ui(x) / u 2 {s)f(s)w(s)ds 



— U2{x) I ui(s)f(s)w(s)ds 

J XI 



where the constants Ci,c 2 G C and Xi,x 2 G (a, b) can be chosen arbitrarily. 
A direct computation shows that u,pu' G ACi oc (a, b) and au = f . □ 



4.2. STURM-LIOUVILLE OPERATORS 



47 



The adjoint operator Since Aq is densely denned by Lemma 4.2.8 the 
adjoint Aq is defined. 

Lemma 4.2.10. The operator A\ is the adjoint of Aq, i.e. Aq = A\. 



Proof. According to the Green's formula eq. (4.1) we have (A u\v)h = 
{u\A\v)h for all u G T>(Aq) and all v G T>(Ai) which means that A± is 
adjoint to A , i.e. A\ C Aq. Thus, it suffices to show that A\ D Aq. Let 
v G V(Aq), i.e. there exists an / G H such that (A u\v)h = (u\f)n for all 
u G T>(Aq). In particular it follows that {Aqu\v) h = (u\f)n for all u G T>(Aq) 
since T>(Aq) C V(A ). Invoking Lemma |4.2.9 we find a v G AC\ oc (a, b) satis- 
fying av = f. Hence, we have (A u\v)h = (u\f)n = (u\(xv)h = {u\av\i) H i = 
(u\A[v\j) h i = (Aqu\v\i) h i for all u G V(Aq). On the other hand we can write 
(Aqu\v) h = (Alu\v\i) H z for all u G V(Aq). This yields {A^u\(v - v)\j) H i = 
for all u G T>(Aq) which shows that (v — v)\j G (imAQ) 1 - = ker A[. Thus, 
a(v — v)\i = which implies av\j = Since this relation holds for all 
/ = [c,d] C (a, b) we infer av = f G H. Hence, we have v,pv' G AC\ oc (a, b) 
and v,av G H which shows that v G T>{A\). Therefore, we have shown that 
Aq <Z A\ and we are done. □ 



The special case limit— point, limit— point We consider the special case 
of a singular Sturm-Liouville operator where both endpoints are limit-point. 

Lemma 4.2.11. Suppose that both endpoints a and b are in the limit-point 
case. Then, the maximal operator A\ is self-adjoint. 



Proof. According to Green's formula eq. (4.1), the operator A\ is symmetric 
and thus, Ai G A\. However, we have Aq C A x which implies A\ C Aq = A\ 
by Lemma |4.2.10| □ 



Corollary 4.2.1. Let both endpoints a andb be in the limit-point case. Then, 
the maximal operator A\ is the closure of Aq. 



Proof. Lemma 
Lemma 4.2.11 



4.2.10 



tells us that Aq 



Ai. This implies A* * = A\ 



and we are done since Aq* is the closure of Aq. 



A 1 by 

□ 



The special case limit— circle, limit— point We study the special case 
of a singular Sturm-Liouville operator on (a, b) where a is limit-circle and 
b is limit-point. The question is what kind of boundary condition one can 
choose in order to construct a self-adjoint operator. 

Lemma 4.2.12. Suppose a is in the limit-circle case and b in the limit-point 
case. Fix a x G T>(Ai) such that there exists av G T>{A\) with [v ,x]p( a ) 
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(such a x exists since a is limit-circle) and define T>(A) := {u G T>(Ai) : 
[u, x]p( a ) = 0}; Au := au for u G T>(A). Then, the operator A : T>(A) C 
H — > H is self-adjoint. 

Proof. For ui,u 2 ,Us,U4 G D(A{) we have the so-called Pliicker identity 

[u 1 ,u 2 ] p (x)[u 3 ,u 4 ] p (x) + [u 1 ,Us)p(x)[u4,u 2 )p{x) + [uuuA p {x)[u 2 ,uz\ p {x) = 0. 

Thus, choosing u±,u 2 G T>(A), u% := \ an d ^4 sucn that [x,u^\p(a) ^ 
we obtain [ui,u 2 ] p (a) = for all Ui,u 2 G T>(A). Hence, the operator A is 
symmetric (b is limit-point) and it remains to show that A* C A. 

We have A Q G A which implies A\ = A^ D A* by Lemma 
v G V(A*) C V(Ai). By definition we have (Au\v)h = (u\A*v)h = {u\A\v)h 
for all u G Thus, the Green's formula yields [u, v] p (a) = for all 

u G T>(A). Note that x ^ ^(^) by definition and hence, the above formula 
with u = x yields [x, v] p (a) = which shows that v G T>(A). Therefore, we 
have A* C A. □ 



4.2.10 



Let 



4.2.4 The Weyl Alternative 

Finally, we come back to the already mentioned relationship between inte- 
grability of solutions of au = and endpoint classification. The following 
theorem, known as the Weyl alternative, is very important for applications 
since it simplifies the endpoint classification for concrete Sturm-Liouville 
operators a lot. 

Theorem 4.2.2 (Weyl alternative). The endpoint a is in the limit-circle 
case if and only if there exist two linearly independent functions U\,u 2 G 
AC\ oc (a,b) with pu'j G AC\ oc (a, b) which belong to H near a Q and satisfy 
(i-a) Uj = (j = l,2). 

Sketch of Proof. Given two functions Ui,u 2 satisfying the properties stated 
in the theorem one can easily construct u, v G T>(Ai) such that u = U\ 
and v = Vi near a (use appropriate cut-off functions). Since U\,u 2 are 
linearly independent it follows that u, v have the same property which implies 
[u,v] p (a) 7^ and hence, a is limit-circle. 

Conversely, let a be limit-circle, i.e. there exist XiV £ ^(^1) such that 
[x,v]p{ a ) 7^ 0- By specifying an appropriate boundary condition at b (or 
none, if b is limit-point) we can construct two self-adjoint operators A x 
and A v where u G V(A X ) implies [u, x]p( a ) = an d u G T>(A V ) satisfies 

1 Onc says that the function u belongs to L 2 w (a, b) near a if there exists a c € (a, b) such 
that u|( a , c ) S L^(a, c). 
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[u, T]] p (a) = (cf. Lemma 4.2.12). We choose an / G H with compact 
support and define u x := {% — A x )~ l f , := {% — A^ 1 f . Since /|( a , c ) = 
for some c G (a, 6), it follows that (z — a)w x |( a , c ) — (* — «)w r? |( ajC ) = 0. By 
playing around with the variation of constants formula one can show that / 
can be chosen in such a way that neither u x |( a , c ) n o r w^|( a>c ) are identically 
zero. Moreover, from [u x ,x\ p {a) = [u v ,r]] p (a) = it follows easily that u x 
and linearly independent. Hence, by solving an initial value problem 

at c we can extend w x |( ,c) an d Wjylra.c) to (a, 6). By construction the resulting 
functions belong to H near a and solve (i — a)u = 0. □ 



Remark 4.2.4. Clearly, Theorem 4.2.2 is equally valid for the endpoint b. 
Moreover, one can show that i can be substituted by any A G C. 
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Chapter 5 



Semigroups and Abstract Wave 
Equations 

In this chapter we prove a generation result for an abstract second order 
Cauchy problem on a Banach space. Then we consider the inhomogeneous 
abstract Cauchy problem and prove its well-posedness under certain assump- 
tions. Next, using a fixed point argument, we show existence and uniqueness 
of solutions of a nonlinear abstract wave equation. 

5.1 The Abstract Wave Equation 

We prove well-posedness of an abstract second order Cauchy problem on a 
Banach space where the involved operator satisfies certain conditions. 

5.1.1 Reduction to a First Order Equation 

Statement of the problem In what follows we will assume that B : 
T>(B) C H — > H is a densely defined closed linear operator on a Hilbert 
space H satisfying (Bu\Bu)h > 7 for all u G T>(B) and a 7 > 0. 

Since B is densely defined, there exists the unique adjoint B*. The fol- 
lowing theorem is due to von Neumann. 

Theorem 5.1.1. Let B : T>(B) C H — > H be a densely defined closed linear 
operator on a Hilbert space H . Then, B*B is self-adjoint and V(B*B) is a 
core of B. 



Proof. See [26], p. 275. 



□ 
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We set A := B*B and consider the abstract second order Cauchy problem 

Mt) = ~M(t) for t > 

^(0) =-00,^(0) =Vi 1 j 

for a function ^ : [0, oo) — > H. 

Function spaces Define Y := T>{B) and ||u||y := a/ (Bu\Bu)h for u G 
V{B). 

Lemma 5.1.1. The normed vector space (Y, || • ||y) is a Banach space. 

Proof. Let (uj) be a Cauchy sequence in (Y, \\ ■ ||y), i.e. {Buj) is a Cauchy 
sequence in if. if is complete and therefore (Buj) has a limit / G if . Since 
< 7~ 1 (-Bm|5m)^ = 7~ 1 ||w||y for all w G Y, (uj) is a Cauchy sequence 
in if as well and hence it has a limit u G if. Therefore, we have — > u and 
-Bwj — > / in if . Since B is closed, it follows that u G T)(B) = Y, Bu = f 
and we have \\u — Uj\\ Y = \\B(u — MjOII-H" = 11/ — Buj\\ H — > 0. □ 

Lemma 5.1.2. 27ie Banach space (Y, || • ||y) zs continuously embedded in 
(if, || • ||#), z.e. the inclusion map i : (Y, || ■ ||y) — > (if, || • ||#) defined by 
i(u) = u for u EY is bounded. 

Proof. Let u G Y and calculate ||i(u)||ff = IMIh ^ 7 _1 I|£ u IIh = 7~ 1 IMIy- 

□ 

Remark 5.1.1. We write Y "—>■ H whenever (Y, || ■ ||y) embeds continuously 
in (H, || -Wh). 

We define V(L) : = V(A) X V(B ) and X := Y x if. Introducing the 
norm ||(u,u)||x : = \/IMI y + IMI# f° r ( u > v ) e II * llx) becomes a 



Banach space thanks to Lemma 5.1.1 We define L : T> (L ) C X — > X by 



L(u,v) := (f, — Au) for (it, t> ) G i?(f) and consider the Cauchy problem 

|u(t) = Lu(t) for * > 
u(0) = u 



for u : [0, oo) — > X and u G X which is (formally) equivalent to eq. (5.1). 
Then the following holds. 



Proposition 5.1.1. The evolution problem eq. (5.2) is well-posed. 
5.1.2 Well— Posedness 



We prove Prop. 5.1.1 by showing that the operator L generates a strongly 
continuous one-parameter semigroup on X, thus we verify the assumptions 
of the Hille-Yosida Theorem. 
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Analytic properties 

Lemma 5.1.3. The operator L is densely defined. 



Proof. According to Theorem 5.1.1, T>{A) is a core of B, i.e. there exists a 
closeable operator C : V{C) C H -> H such that V(C) = V(A) and C = B. 
We define the graph G(B) of B by G(S) := {(u, Bu) G H x H : u G £>(£)}. 
We equip H x H with the norm \\{u, v)\\hxh '■= vlHl5 + IMIh an d hence, 
(if x if, || • ||ifxi?) is a Banach space. The fact that 5 is closed is equivalent 
to G(B) being a closed subset of (H x H,\\ ■ \\hxh), i-e. = 
C(C) is a subset of and G(C) = G(B). 

Now let u G y = 15(5). Then, (u, Bu) G and, since G(C) is 

dense in G(B), there exists a sequence ((uj,Cuj)) C G(C) with (uj,Cuj) — ► 
(u, 5m) m. H x H. Observe that Cttj = £?itj for all «j G T>(C) since S is an 
extension of C. Therefore we have 



\u — it 



12 < \\Bu - BujWjj + \\u-uA\ 2 H = \\( U j,Cuj) - {u,Bu)\\ 2 HxH -> 



j II y 



which shows that is dense in (y, || • ||y). 

By assumption we have T>(B) dense in if and hence X'(L) = X'(A) x T>(B) 
is dense in X = Y x H. □ 

Lemma 5.1.4. The operator L is closed. 

Proof. Let ((uj,Vj)) C T>{L) be a sequence with (uj,Vj) — > (u, u) in X and 
L(uj,Vj) — > (/, p) in X, i.e. Wj — > w in y, u,- — > u in H, Vj — > f in Y and 
—Auj — > g in if. Convergence in Y implies convergence in H since Y ^ H 



by Lemma 5.1.2 and hence we also have Vj — > / in H which implies / = v 
by uniqueness of limits. Therefore, v EY = T>(B). 

By the same argument we have Uj — > it in if and together with —Auj — > g 
in ii and the closedness of A we conclude w G X^A) and — Au = g. 

Hence, we have shown that (u,v) G T>(L) = T>(A) x T>(B) and L(u,v) = 
{f id) which proves closedness of L. □ 

Spectral properties 

Lemma 5.1.5. The spectrum a (A) of A satisfies o~(A) C [7,00). 

Proof. A is self-adjoint and satisfies (Au\u)h = (Bu\Bu) H > j(u\u)h for all 



u G T>(A). Invoking Prop. 4.1.1 finishes the proof. □ 



Lemma 5.1.6. The spectrum o~(L) of L satifies o~(L) C {A G C : —A 2 G 
a{A)}, i.e. A G a(L) -A 2 G a{A). 
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Proof. Consider the equation 



v J \ 9 

for (it, v) G T>(L) and (/, g) G X which is equivalent to the system 

v = Xu — f 
(X 2 + A)u = Xf + g ■ 

Suppose -A 2 G p(A). Then X 2 + A : V(A) C H -> # is invertible and 
we define u : = (A 2 + A) _1 (A/ + g) and t> := Xu — f for given f <EY C H and 
g E H. Then, G Z>(A) x F = V(L) and (A - = (f,g) which 

shows that A — L is surjective. Let (it, i>) G P(i^) and L(u, v) = 0. It follows 
that (A 2 + A)u = which implies u = since —A 2 G p(^4). We conclude that 
v = Ait = and hence (u, t>) = (0,0) which proves injectivity of A — L and 
therefore, A — L is bijective which implies A G p(£). 

Thus, we have shown —A 2 G p(A) =^> A G p(L) which is equivalent to 
A G cr(L) =>- —A 2 G cr(v4) and this is the claim. □ 

Corollary 5.1.1. The interval (0,oo) is contained in the resolvent set of L, 
i.e. (0, oo) c p(L). 



Proof. Suppose A G (0, oo) and A ^ p(L), i.e. A G cr(L). From Lemma 5.1.6 
it follows that — A 2 G cr(A) but this is a contradiction to Lemma 5.1.5 which 



states that a (A) C [7, 00) and 7 > 0. □ 

Lemma 5.1.7. The resolvent Rl of L satisfies \\Rl(^)\\b<x) < k f or a ^ 
A>0. 



Proof. Let A G (0, 00). From Corollary 5.1.1 we know that -Rl(A) = (A — 
L)^ 1 G B(X) exists and we set (u,v) := Rj J (X)(f, g) for (/, g) G X. Then we 
have Au + Xv = g. We take the inner product with v and obtain (Au\v)h + 
'MMIh = {q\ v )h- Substituting v = Xu — / yields X(Au\u)h + A||i>||g = 
(g\ v )H + (Au\f)n- Using the Cauchy-Schwarz inequality we estimate 

K\\Bu\\ 2 H + \\v\\ 2 H ) = (g\v) H + (Au\f) H = (g\v) H + (Bu\Bf) H 

< \\ 9 \\h\\v\\h + \\Bu\\ H \\Bf\\ H < (\\Bu\\ 2 H + \\v\\ 2 H ) ll2 {\\Bf\\ 2 H + WgWl) 1 / 2 

which yields A || (it, f) || x < || (/, g)\\x and this is equivalent to ||i2z,(A)(/, g)\\x < 
l\\(f,g)\\x- Since (f,g) G X was arbitrary this implies ||-Ri(A)||B(x) < 

□ 
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Generation of the semigroup Taking together Lemmas 5.1.3 5.1.4 5.1.7 



and Corollary 5.1.1 we have shown that L satisfies the requirements of the 
Hille-Yosida Theorem and hence, L generates a strongly continuous one- 



parameter semigroup S : [0, oo) — > B(X) on X satisfying ||5'(t)||e(x) < 1 
for all t > 0. Thus, the abstract Cauchy problem eq. (5.2) is well-posed as 
claimed in Prop. |5.1.1 



Summary We summarize the results of this section in the following theo- 
rem. 

Theorem 5.1.2. Let H be a Hilbert space and B : T>(B) C H — ► H a 

densely defined closed linear operator which satisfies (Bu\Bu) H > j(u\u)h 
for all u e V(B) and a 7 > 0. Define A := B*B, V(L) := V(A) x V(B), 
\\ U \\ Y ■= \\Bu\\ H for u E V(B), Y := V(B), X := Y x H and L : V(L) C 
X ^X by 



Then, L generates a strongly continuous one-parameter semigroup S : [0, 00) - 
B(X) on X satisfying \\S(t)\\is(x) < 1 for allt > 0. In particular, the abstract 
evolution problem 

f f t u(t) = Lu(t) fort>0 
\ u(0) = u 

for u : [0, 00) — > X and u G X is well-posed. 



5.1.3 A Simple Example 

The wave equation As a simple example we apply the generation result 
to the one-dimensional wave equation with Dirichlet boundary conditions, 
i.e. we consider the Cauchy problem 

ip tt (t, x) = ip xx (t, x) for (t, x) e (0, 00) x (0, 1) 
tp(t,0) = 0,i>(t,l) = for t > 
^(0, x) = if)o(x),i/)t{Q, x ) — i>i{ x ) f° r x G [0, 1] 



Operator formulation As a Hilbert space we take H := L 2 (0, 1) and set 
V(B) := {u e H l {Q, 1) : u(0) = u(l) = 0}. We define B : V(B) C H -»• H 
by Bu := v! for u E V(B). 

Lemma 5.1.8. The operator B : T>(B) C H ^ H is closed. 



Proof. See e.g. [34], p. 29, Example 3.10. 



□ 
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5.2. THE CASE 7 = 



We claim that -B*Bu = u" for u G C£°(0, 1) and thus, 



^ tt (t) = -B*Bij){t) for t > 
^(0) =M t (0) 



is an operator version of the one-dimensional wave equation. 

Lemma 5.1.9. For u G C^°(0, 1) we have -B*Bu = u" . 

Proof. Integration by parts immediately yields (Bu\v)h = (u\ — Bv)h for all 
11,1)6 V(B) and thus, -B is adjoint to B, i.e. -B C B*. Let u G Q°(0, 1). 
Then we have -B*Bu = -B*u' = Bu' = u". □ 

Now we show well-posedness using Theorem |5.1.2| 

Lemma 5.1.10. The operator B satisfies (Bu\Bu)h > {u\u)h f or a ^ u ^ 
V(B). 



Proof. Let u G T>(B) and observe that 



\u(x)\ 



u'(s)ds 



< 



\u'(s)\ds < 



1 \ 1/2 

\u'(s)\ 2 ds 







W Iff 



for all x G [0, 1] by Cauchy-Schwarz. Integration yields (u\u)h 



\u\ 



H 



U 



'112 
Iff 



(-Bul-Bu)u. 



< 
□ 



Thus, our previous results in this section (Theorem 5.1.2) imply that 
the first-order operator evolution problem which is associated to the one- 
dimensional wave equation is well-posed. 



5.2 The Case 7 = 



In Theorem 5.1.2 the bound 7 is assumed to be strictly positive. The obvious 
question is whether this requirement can be weakened. In this section we 
discuss the case 7 = 0, i.e. we assume that B : T>(B) C H — > H is a densely 
defined closed linear operator on a Hilbert space H satisfying (Bu\Bu)h > 0. 
The main difficulty one encounters in this case is the fact that the normed 
vector space (Y, || • \\y) defined by Y := V(B) and ||w||y := is not 

complete. Hence, the semigroup cannot act on Y x H (which is the energy 
space). To go around this problem one has to introduce a slightly different 
Banach space in order to recover the existence of the semigroup. This yields 
a well-posedness result but the growth estimate becomes worse. However, 
this can partly be compensated by energy conservation. 
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5.2.1 Generation of the Semigroup 

Without loss of generality we can assume B to be self-adjoint (B*B is self- 
adjoint and nonnegative and hence, there exists the self-adjoint nonnegative 



square root of B*B, cf. sec. 4.1). The idea is to consider the operator B + ie 



for e > instead of B and "repair" this defect by a bounded perturbation. 

Splitting of the operator We define the operator L by L := L + V 
where L Q and Ll are given by 



L :-{ B *l_ £ 2 I) andZ/:-( £ ° 2 ° 

We show that L generates a strongly continuous one-parameter semigroup 
and apply the bounded perturbation theorem. 

Theorem 5.2.1 (Bounded Perturbation Theorem). Let Lq : T>(L$) C X —>■ 
X be the generator of a strongly continuous one-parameter semigroup So : 
[0, 00) — > B(X) on a Banach space X satisfying 



\\So(t)\\ B (x) < e 

for an to E R and all t > 0. If V E B{X) then L := L + V : V(L ) C 
X —>■ X is the generator of a strongly continuous one-parameter semigroup 
S : [0, 00) — ► B(X) satisfying the estimate 

||^(t)|| B (x)<e^ + ll L 'll^))* 

for all t > 0. 

Proof. See [17], p. 158. □ 
Generation of the semigroup Since ((B+ie)u\(B+ie)u)jj = (Bu\Bu) H + 



e 2 (u\u)h > e 2 (u\u)h for all u G T>(B), Theorem 5.1.2 implies that the opera- 
tor L generates a strongly continuous one-parameter contraction semigroup 
on X e . Note that the space X e on which the semigroup acts depends on 
e since the norm || • is given by \\(u, v)\\xs = \\{B + ie)u\\ 2 H + \\v\\ 2 H for 
(u, v) G T>(B) x H\ The perturbation V satisfies 

\\L'(u,v)\\ X e = \\(0,e 2 u)\\ X c =e 2 \\u\\ H < e\\(B + ie)u\\ H < e\\(u,v)\\ X e 

for all (u,v) G X e . Thus, applying the bounded perturbation theorem we 
conclude that L = L + V generates a strongly continuous one-parameter 
semigroup S £ on X s satisfying || S' e (t) < e£t f° r a h t > 0. 
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5.2. THE CASE 7 = 



5.2.2 Removing the e— Dependence 

We show how to remove the bothersome e-dependence of the underlying 
Banach space. We define X := X 1 , i.e. X = V(B) x H and v)\\ x = 
\\(B + i)u\\ 2 H + \\v\\h for (u,v) G X. Since \\(B + ie)u\\ 2 H = \\Bu\\ 2 H + e 2 \\u\\ 2 H , 
we have 

e 2 \\{B + i)u\\ 2 H < \\(B + ie)u\\ 2 H < \\(B +i)u\\ 2 H 
for all u G T>(B) and any < e < 1. Thus, we conclude that 

e|M|x < ||u||x« < ||u|U 

for all u G X and any < e < 1. Hence, we infer the estimate 

£||S e (t)u|| x < ||S e (t)u|| xe < e £ *||u|| xe < e et ||u|| x 

for all u G X which implies HS^t) \\b(X) < -e et for any < e < 1. Fur- 
thermore, all S* £ (for different e) have the same generator L which does not 
depend on e and hence, they all coincide. Thus, we can drop the superscript 
e and, for any < e < 1, we obtain ||£(£)||b(x) < -e e * for all t > 0. Finally, 
choosing £ = | for i > 1 we arrive at 

\\S(t)\\ BiX) <tfoTa31t>l. (5.4) 

5.2.3 Energy Conservation 

We define the energy "norm" on X by || (u, v)\\ 2 E := ||.Bu||^ + ||f ||# for (m, t>) G 
X [j] Thus, we have u)||.e < \\(u, v)\\x for all (u,v) G X. Consider the 
time evolution of initial data (u , Vq) G T>(L). According to semigroup theory, 
the solution (u(t), v(t)) := S(t)(u , v ) stays in T>(L) for all t > and satisfies 
the equation 

r |«(t) = v(t) 

in the strong sense, i.e. it is a classical solution. Recall that the norm 
||(wo,«o)IU is § iven °y ll( M o,^o)|lx = \\( B + O m oIIh + IKHI and note that 
11(5 + i)uo||jf = \\Bu \\ 2 H + \\u \\ 2 H . Hence, we can write Q ||(t*o,«o)||i = 
|| «o || b + II^oIIh an d this implies that the derivative 4:u(t) exists with respect 
to the graph norm of B, i.e. lini/^o h~ l \\u(t + h) —u(i)\\B exists for all t > 0. 

^dn general this is only a seminorm since might be an eigenvalue of B and thus, there 
might exist elements (u,v) € X with u ^ but \\(u, v)\\e = 0. 

2 Recall the definition of the graph norm || • \\g given by \\u\\b := ||£?m||h + \\u\\h for 
ueV(B). 
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Therefore, since B is a closed operator, we have j^Bu(t) = Bj- t u(t). Having 
these observations in mind we readily calculate 

j t [(Bu(t)\Bu(t)) H + (v(t)\v(t)) H ] = 2Re(Bu(t)\Bu(t)) H 

+ 2Re(v(t)\v(t)) H = 2Re(Bv(t)\Bu(t)) H - 2Re(v(t)\B*Bu(t)) H = 

for alH > where' := -|. Thus, the function t \— > \\S(t)(u , v )\\e is constant. 
This shows that for all classical solutions the energy is conserved. 

Now let u G X. Then, there exists a sequence u j C T>{V) such that 
u j — > u in X. Since S{t) is bounded we conclude that S(t)u j — > S(t)uo in 
X for any t > 0. Recall that < \\-\\ x which shows that S^Uoj — > 5 , (t)u 
with respect to \\ ■ \\e as we U- I n particular we have 

H-S^UoIIe = lim ||S , (t)u ,-|U = lim \\u 0j \\ E = ||u ||e 

for all t > 0. This shows that energy conservation holds for generalized 
solutions as well. 



5.2.4 Summary 

We formulate the results of this section as a theorem. 

Theorem 5.2.2. Let H be a Hilbert space and A : V(A) C H — > H a 
self-adjoint operator satisfying (Au\u) H > for all u e f(^4)- Define X : = 
V{A 1 ' 2 ) x if and ||(u,u)|ft := ||w||^i /2 + \\v\\ 2 H for (u,v) E X. Then, the 
operator L : V(L) Cl^I, defined by V(L) := V(A) x V(A 1/2 ) and 




generates a strongly continuous one-parameter semigroup S : [0, oo) — > B(X) 
on X satisfying 

\\S(t)\\ B{x) <t 

for all t > 1 . 

Furthermore, if (uo,v ) e X and (n(t),n(t)) := S , (t)(no ? ^o); the function 
t ^ \\A 1/2 u(t)\\ H + \\v(t)\\ H is constant for all t>0. 



5.3 The Inhomogeneous Problem 

As a next step we discuss semigroup theory for inhomogeneous evolution 
problems. This approach relies on the notion of an integral of a semigroup 
whose definition requires a little background in measure theory which will be 
outlined first. 
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5.3. THE INHOMOGENEOUS PROBLEM 



5.3.1 Basic Aspects of Measure Theory 

We briefly discuss the construction of the Lebesgue integral for Banach space 
valued functions on intervals. An introduction to measure theory can be 
found in e.g. [4]. For the definition of the Bochner integral we also refer 
to [59]. 

o"— algebra, measure Let f2 be a set and A C V{Vt) where V{Vt) denotes 
the power set of Q, i.e. the set of all subsets of Q. A is said to be a a-algebra 
if 

• G A 

• A e A =>- Q\A e A 

• A n e A for n e N [j neN A n eA 

Let C C V(£l). The smallest <r-algebra which contains C is denoted by a{C). 
We remark that a(C) exists for any C C V(VL) since V(Vl) is a cr-algebra 
itself. 

A mapping (j, : A — > [0, oo] is called a measure on A if /i(0) = and 
MLLn A) = E^=i MA.) for A e A n G N and n A,- = for i ^ j. 

Borel a - — algebra, Lebesgue measure Let f2 = R and set 7i := {0} U 

{(a,b] : a, b G K, a < 6}. Then, i3 := cr(7i) is called the Borel a-algebra. One 
can show that there exists a unique measure A on B such that \((a,b}) = b — a 
for all half-open intervals (a, b] G 7i. The measure A is called the Lebesgue 
measure. A similar construction can be applied to Q = R n . 

Measurable functions, simple functions Let B be the Borel a algebra 
on 1, (1 a set and A C V(VL) a cr-algebra. A function / : R — > f2 is said to 
be measurable if / _1 (A) G £> for all Ae A 

Now we restrict ourselves to real-valued functions. We denote the char- 
acteristic function of a set A by xa, i-e. 

, x _ f 1 for i e A 
~ \ for x £ A 

A function u : M. — > R is called simple if there exist A,- G £> for j G N which 
satisfy IJjeN A = ^ anc ^ A fl A,- = for i ^ j and w = Y^j=i c jXAj where 
cj G [0, oo), nGN. 
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Lebesgue integral One defines the Lebesgue integral J udX over a simple 
function u : R — > R with u = YTj=i c jXA ] by 



/ 



where A denotes the Lebesgue measure. 

Let / : R — > R be a nonnegative measurable function. One can show that 
there exists a sequence of simple functions (u n ) such that u n < u n+ i for all 
nGN and /(x) = lim^oo u n (x) = sup n6N u n {x). Then one defines 



/fd\ := sup / u n dX. 
neN J 



For general measurable functions / : R — > R one sets / _ (x) := max{-/(i), 0} 
and f + {x) := max{/(x), 0}. One easily shows that the nonnegative functions 
/~ and / + are measurable. If f f + dX < oo and f f~dX < oo one says that 
/ is integrable and defines 



J fdX := J f + dX - J f-dX. 



The generalization to complex-valued functions is obtained by consider- 
ing the real and imaginary parts separately. 

Lebesgue's theorem on dominated convergence We have the follow- 
ing important convergence theorem. 

Theorem 5.3.1 (Dominated convergence theorem). Let (Q,A,fj) be a mea- 
sure space (i.e. A C V(£i) is a a -algebra and \i a measure on A) and 
f,f n :Q^WL measurable functions. Furthermore, assume that f n — > / 
pointwise almost everywhere and there exists a nonnegative measurable func- 
tion g with f gdji < oo and \f n \ < g almost everywhere for all n € N. Then, 
f n is integrable for all n G N and we have 

lim / f n dfi = / fdfi. 



Banach space valued functions Let X be a Banach space. Again, B 
denotes the Borel cr-algebra on R and A the Lebesgue measure. 

A function u : R — > X is said to be simple if there exist Aj e B, j G N with 
UjeN = R, Ai fl Aj = for i ^ j and Xj G X such that u = YTj=i x iX.A v 
n G N. 
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5.3. THE INHOMOGENEOUS PROBLEM 



Measurable functions A function u : R — > X is said to be weakly mea- 
surable if, for any / G X*, the complex-valued function t i— > f(u(t)) : K — > C 
is measurable, u is said to be strongly measurable if there exists a se- 
quence (u n ) of simple functions and a i?o G £> with X(Bq) = such that 
||w(£) — u n (t)\\ x — > for all t G K\-B , i.e. u n (i) — > it(£) in X for almost all 
tel. It turns out that, if X is separable, the notions "weakly measurable" 
and "strongly measurable" are equivalent. 



The Bochner integral Let u : K — > X be a simple function with u 
Sj=i ^jXAj- Then, the Bochner integral J wdA G X of u is defined by 



/ 



3=1 

Let u : R — > X be a strongly measurable function. One can show that the 
real-valued function £ i— > ||tt(£)||x : K. — * K. is measurable. The function 
u is said to be Bochner integrable if there exists a sequence (u„) of simple 
functions such that 

lim / \\u — u n \\xd\ = 0. 

n— >oo / 

It turns out that this condition implies the existence of the limit of the 
sequence (J u n dX) in X. Then, one defines the Bochner integral of u by 



/wdA := lim / u n d\. 



Theorem 5.3.2 (Bochner). Let (X, || ■ \\ x ) be a Banach space. A strongly 
measurable function u : R — > X is Bochner integrable if and only if the 
function t \— > : R — *• R integrable. In this case the estimate 



holds. 



J udX 


X <-J 







Notation To improve readability we will adapt the usual " Riemann-like" 
notation for integrals, i.e. consider a function u : I — > X where I C 1 is 
some interval and X a Banach space. Then we define a function u : R — > X 

w \ G X for t <£ I 

and set 

u(t)dt := I udX. 
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5.3.2 The Abstract Problem 

We consider the inhomogeneous abstract Cauchy problem, see [17]. We also 
refer to [44] for an extensive treatment of this subject. 

Statement of the problem Let (X, || ■ \\ x ) be a separable [^Banach space, 
L : T>(L) C X —>■ X a linear operator and / : [0, oo) — > X a function. 
Consider the abstract inhomogeneous evolution problem 

f t u(t) = Lu{t) + fit) for t > 
u(0) = u 

for a function u : [0, oo) — > X and initial data u G X. Now we extend the 
notion of well-posedness to the inhomogeneous problem. 

Definition 5.3.1. A function u : [0, oo) — > T>(L) C X is called a classical 



solution of eq. (5.5) if it is continuously differentiable and satisfies eq. (5.5). 



The abstract evolution problem eq. (5.5) is said to be well-posed if for 



every uq G T>(L) there exists a unique classical solution u : [0, oo) — > V(L) 



of eq. (5.5) such that the mapping u i— »■ tt(t) : X>(£) C X — > X is uniformly 



continuous in t on compact intervals [0, T] for any T > 0. 

Construction of solutions Suppose that L generates a strongly continu- 
ous one-parameter semigroup S : [0, oo) — > Bpf) on X satisfying HS'(t) ||s(x) < 
C(t) for a continuous positive function C and all t > 0, i.e. the associated ho- 
mogeneous problem (/ = 0) is well-posed. Assume further that / : [0, oo) — ► 
X is continuous. Then, the function s i— > 5*(t — s)f(s) : [0, t] — » X is continu- 
ous and hence, for any (7 G X*, the function s 1— > ^(S 1 ^ — s)f(s)) : [0, t] — >• C 
is continuous and thus measurable. This shows that s 1— > 5(t — s)f(s) : 
[0, t] — > X is weakly measurable. Since X is separable we conclude that 
s 1 — > S(t — s)f(s) is strongly measurable. Moreover, for s G [0, t] we have 
\\S(t-s)f(s)\\ x < C(t - s)\\f(s)\\ x < sup sem C(s)\\f(s)\\ x < 00 which 
shows that 

\\S(t-s)f(s)\\ x ds 





exists. Thus, Bochner's Theorem implies that s i— > S(t — s)f(s) is Bochner 
integrable. 
Define 

u(t) := 5(t)«o + / 5(t - s)f(s)ds. (5.6) 



Then, u is called the mzW solution of eq. (5.5). 



3 The assumption of separability is introduced here for convenience only, it is by no 
means necessary. 



64 



5.4. A NONLINEAR PROBLEM 



Classical solutions It is easy to show that every classical solution is a 
mild solution. In particular, this implies that a classical solution is unique. 
It turns out that, if / satisfies certain additional regularity conditions, a 



classical solution can be constructed by formula eq. (5.6). If, for instance 



/ has a weak t-derivative which is integrable and uq G T>(L) then u defined 



by eq. (5.6) is the classical solution of eq. (5.5). More precise, we have the 



following theorem. 

Theorem 5.3.3. Let L be the generator of a strongly continuous one-parameter 
semigroup S : [0, oo) -> B(X) on X. If u G V(L) and f G W 1 ' 1 ^, oo), X) 



then u : [0, oo) — > X defined by eq. ( 5.6) is the unique classical solution of 
eq. (5~5\ ). 



Proof. See [17], p. 439. □ 

Remark 5.3.1. The space 11 /1 ' 1 ([0, oo), X) is a Sobolev space of Banach space 
valued functions u : [0, oo) — > X analogous to iy 1 ' 1 (0, oo) for complex-valued 
functions. One can define W /1 ' 1 ([0, oo), X) as follows. Consider the vector 
space of functions u : [0, oo) — > X such that u is integrable (in the sense of 
Bochner) and there exists an integrable function v : [0, oo) — > X such that u 
can be written as 

u(t) = u(t ) + / v(s)ds 

for a t G [0, oo). Two such functions are identified if they coincide for almost 
all t G [0, oo). Then, the quotient space defined by this equivalence relation 
is denoted by W 1 ' 1 ^, oo), X). 

Well— posedness of the inhomogeneous problem We summarize the 
results of this section in a theorem. 

Theorem 5.3.4. Let X be a separable Banach space, L : T>(L) C X —>■ 
X a linear operator and f G W /1 ' 1 ([0, oo), X). If L generates a strongly 
continuous one-parameter semigroup S : [0, oo) — ► B(X) then the abstract 



solution u is given by eq. (5.6) for uq G T>(L) 



inhomogeneous Cauchy problem eq. (5.5) is well-posed. Its unique classical 



5.4 A Nonlinear Problem 



In this section we consider existence and uniqueness of an abstract nonlinear 
evolution equation for a function u : [0, oo) — * X on a separable Banach 
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space A given by 

f f t u{t) = Lu{t) + g(u(t)) (57) 
1 w(0) = u 

where uq G A and g : X — > A is Lipschitz-continuous, i.e. there exists a 
C > such that\\ g(x)—g(y)\\x < C||rr— y||x for all x, y G A. We assume that 
L generates a strongly continuous semigroup S 1 : [0, oo) — ► -B(A) satisfying 
\\S(t)\\ x < e ut for all t > and an cu G E. 



5.4.1 Existence and Uniqueness 

Let T > and denote the vector space of continuous functions u : [0, T] — > A 
by C([0,T],A). For u G C([0,T],A) the real-valued function t ^ ||w(i)|U 
is continuous and hence attains its maximum on the compact interval [0,T]. 
We define a norm || • ||c([o,t],x) on C([0,T], X) by 

||«||c([o,T],x) := sup ||w(t)||x- 
te[o,T] 

It is easy to show that (C([0,T],X), \\ ■ ||c([o,t],x)) is a Banach space. 

Now let Mo G A be fixed and define a mapping K : C([0,T], X) — > 
C([0,T],A)by 

(K(v))(t) := S(t)u + [ S(t- s)g(v(s))ds (5.8) 

Jo 

for t G [0, T]. The function s i— >■ S(t—s)g(v(s)) is continuous since g : X — > X 
is continuous. Hence, it is Bochner integrable and the above definition of K 
makes sense. 

Note that the function u := Kv is the unique mild solution of the linear 
inhomogeneous problem 

£u(t)=Lu(t)+g(v(t))foTte[0,T\ 
u(0) = u 



as discussed in sec. 5.3 Thus, a fixed point u of K (u = K(u)) is a mild 



solution of the nonlinear problem eq. (5.7). Hence, to show existence of 



solutions of eq. (5.7) it suffices to show existence of a fixed point of the 



mapping K. To do so we will invoke the Banach fixed point theorem. 

Theorem 5.4.1 (Banach fixed point theorem). Let (Z, || ■ \\z) be a Banach 
space and K : Z — » Z a mapping that satisfies 

\\K{x)-K{y)\\ z < k\\x-y\\ z 



for all x,y G Z and a k < 1. Then, there exists a unique fixed point of K. 
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With the help of this theorem we can prove existence and uniqueness of 



mild solutions of eq. 5.7 for small times. 

Lemma 5.4.1. There exists a T > such that the mapping K : C([0,T],X) 



C([0,T],X) defined by eq. (5.8) has a unique fixed point. 
Proof. Let T > 0, u, v G C([0,T],X) and observe that 



\\K{u)-K{v)\\ c{[W) < sup / \\S(t-s)[g(u(s))-g(v(s 

te[o,T] Jo 



x 



ds 



< Ce' 



where we have used ||S(£)||e(x) — e 
continuity of g. Thus, we have 



u{s) - v{s)\\ x ds < CTe" T 

* for all t > and the Lipschitz- 



sup \\u(s) - v(s)\\x 
se[o,T] 



\K(u) - K{v) 



\C([0,T],X) 



< CTe wT \ 



u - V\\c([0,T],X)- 



Hence, choosing T such that Te^ T < 1/C shows that K satisfies the con- 
traction property required by the Banach fixed point theorem. Therefore, 
for such a T, K has a unique fixed point. □ 



Thus, we have shown that eq. (5.7) has a unique solution u : [0, T] — » X 
for any T satisfying Te wT < 1/C. 

From the local result we immediately obtain a global result. 

Lemma 5.4.2. There exists a unique global (mild) solution u G C([0, oo),X) 



of eq. (5.7). 



Proof. Applying Lemma |5.4.1| we obtain a constant T > and a local-in- 
time solution u\ G C([0,T], X) of eq. (5.7) for given initial data u G X. 
Now we set w := Ui(T) and apply Lemma 5.4.1 again to obtain a solution 
u 2 G C([0,T],X) satisfying 1*2(0) = U\(T). Repeating this process yields 
a sequence (uj) of solutions of eq. (5.7) in C([0,T], X) satisfying Uj(0) = 
%-i(T) for j G N. Now we define u(t) := Uj(t-{j-l)T) for t G [(j-l)T,jT), 
j G N. By construction, w G C([0, 00), X), u satisfies = Lw(i) for each 

£ > and w(0) = w - □ 



5.4.2 Dependence on Data and Growth Estimates 



By construction, the global solution u G C([0, 00), A) of eq. (5.7) satisfies 
the equation 

u{t) = S{t)u + S(t - s)g(u(s))ds (5.9) 
Jo 
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for all t > 0. In order to obtain an estimate for u(t) in terms of the initial 
data u we apply Gronwall's inequality. 

Lemma 5.4.3 (Gronwall's inequality). Let £ : [0, T] — > R be a nonnegative 
integrable function which satisfies 

m<c 2 + c 1 [\( s )ds 

Jo 

for constants C\,C2 > and almost allt £ [0,T]. Then, 

Z(t) <C 2 (1 + dte Clt ) 

for almost all t e [0,T]. 

Proof. See e.g. [19], p. 625. □ 

Lemma 5.4.4. The unique global solution u G C([0, oo),X) satisfies the 
estimate 

\\u{t)\\ x < e<ik||x (l + Cte^e cte ^ 
for all t > and a constant C > 0. 



Proof. Using eq. (5.9), \\S(t) \\ B (x) < for all t > and the Lipschitz- 



continuity of g we readily estimate 

h(Ollx<e^|| M olU + Ce^ f \\u(s)\\ x ds 



for all < t' < t. Application of Gronwall's inequality yields 

!K*')IU < e^||n || x (l + Ct , e- t e ct ' e ") 

for all < t' < t. Setting t' — t we infer the desired result since t > is 
arbitrary. □ 



Although the estimate stated in Lemma 5.4.4 is very weak, it suffices to 



ensure continuous dependence on the initial data. 
5.4.3 Summary 

We reformulate the results on the nonlinear problem as a theorem. 
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Theorem 5.4.2. Let L : T>{L) Cl->I be a linear operator on a Banach 
space X which generates a strongly continuous one-parameter semigroup S : 
[0, oo) — > B(X) satisfying S(t) < e wt for all t > and an uo 6l, Moreover, 
let g : X — > X fre a Lipschitz-continuous function. Then, the nonlinear 
abstract evolution problem 



for uo E X has a unique mild solution u : [0, oo) — > X which depends contin- 
uously on uq. 




Chapter 6 



The Cauchy Problem for Wave 
Maps 



We formulate the Cauchy problem for the wave maps system eq. (2.1 ) 

□$ A + ^ v T A BC {<$>)(d^ B )(d v <5> c ) = 

and state some basic results. As a technical requirement we need fractional 
Sobolev spaces which are introduced first. Then, we state a general result 
concerning local well-posedness of nonlinear wave equations which can be 
applied to the wave maps system. We also mention some recent developments 
in connection with global existence for solutions with small data. 



6.1 Fractional Sobolev Spaces 

We introduce noninteger Sobolev spaces by using the Fourier transform (cf. 
e.g. [19], [59]). 

Lebesgue spaces We have already defined Lebesgue spaces for functions 
defined on intervals. The generalization to complex-valued mappings on 
open subsets of M. n is similar and goes as follows. Let U C lR n be open and 
consider the set C^°(U) of smooth functions from U to C having compact 
support. For p > 1 we define a norm || • \\lp(u) on C£°(U) by 

||w||ip(iO := (J \u(x)\ p d n x 

where integration is understood with respect to the ordinary Lebesgue mea- 
sure on W 1 . The Lebesgue space LP(U) is defined as the completion of C^°(U) 
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with respect to || • ||lj>([/)- We also define local versions L\ oc (U) by 
L\ oc {U) := {u G L P (V) : V C U, V compact}. 

The Fourier transform Now suppose u G L 2 (IR n ) n LMW 1 ) and define 

where £ ■ x := YT j= i & & for x = (x\ . . . , x n ), £ = (£\ . . . , C) G R n Since 
u G L 1 (IR ri ) and |e~^' x | = 1, u is well-defined and it is called the Fourier 
transform of u. It turns out that the mapping T : u h-* u can be extended 
to L 2 (IR n ). We list some important properties. 

• JF : L 2 (IR ra ) — > L 2 (IR n ) is an isometric isomorphism, i.e. it is linear, 
invertible and norm-preserving (Plancherel's theorem). 

• For u G L 1 (R") n L 2 (lR n ) the inverse J 7-1 is given by 

• n T maps derivatives into multiplication", i.e. TD a u{^) = (^) a jFw({) 
if D a u G L 2 (R n ) where a = (a u ...,«„) G N£, -D a- u := d" 1 • • • d^u 
and r := 1 1 '' ,(£•'')" for f = (f, ..,(")£ M n . 

Fractional Sobolev spaces For s G K we define the Sobolev space i7 s (M n ) 
by 

# s (R n ) := {« G L 2 (M n ) : f h-> (1 + |£| s )«(£) G L 2 (M")} 

where |£| is the Euclidean norm of the vector £ G M n , i.e. |£| 2 := YTj=\ \£ j \ 2 - 
Furthermore, we set 

1/2 



M| h .<r») := ( [ Ki + ieiXOI 2 ^) 



for u G H s (R n ). Note that H°(R n ) = L 2 (R n ). Equipped with this norm, 
H s (W n ) becomes a Banach space. Another commonly used notion is the 
homogeneous Sobolev space ii s (IR n ) which is defined as 

H s (R n ) := {u G L 2 (R n ) : £ i-> |£| s w(0 G L 2 (M n )} 

and ^ ^ 

Ha.(r-) := f / \z\ 2s \m\ 2 d n z) . 

For brevity we will write H s instead of H s (R n ). 
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6.2 Local Well— Posedness 

We state some known results concerning local well-posedness of nonlinear 
wave equations and wave maps. 

6.2.1 Semilinear Wave Equations 

We consider the Cauchy problem for a nonlinear wave equation of the form 

□V> = dtp), ip\ t=0 = ip , d t ip\ t =o = fa (6.1) 

for a function ip : R x R n — ► R where fa, fa : R n — > R, / is smooth and 
satisfies /(0) = 0. f{ip, dip) is a shorthand notation to indicate that / might 
depend on the function ip and its first partial derivatives. The necessity 
to prescribe ip\t=o and d t ip\t=o as initial data follows from the fact that the 
equation is second order in time: One has to know the function and its first 
time derivative on the initial surface t = to be able to formally calculate 
higher derivates. 



By a (local) solution of the Cauchy problem eq. (6.1 ) we mean a function 



ip e CQO,^,^ 1 ) with d t ip G (C[0,T},L 2 ) that satisfies 



T 







ip{t,x)U(p{t,x)d n xdt - / d t ip(t,x)(p(0,x)d n x 



T 



+ / fat, x)d t (p(0, x)d n x = / / f{ip,dip){t,x)y{t,x)d n xdt (6.2) 

JR n JO JR 71 

for all test functions ip G C C °°((-T,T) x R n ) and a constant T > (which 
might depend on t/> and ^i)- F° r simplicity we write ip(t,x) instead of 
ip(t)(x). We state the classical local well-posedness theorem. 

Theorem 6.2.1. Let f : R x R n -> R 6e smooi/j and /(0) = 0. T/jen ; /or 
an?/ s > | + 1 and (fa, fa) G -£P x H^ 1 there exists a T > suc/i i/iat i/ie 
Cauchy problem 

Dip = f(ip, dip), ip\ t=0 = ip , d t ip\ t =o = fa 

has a unique solution (ip,d t ip) G C([0, T], H s ) x C([0, T], H s ~ l ). Moreover, 
T is bounded below by a strictly positive continuous function of \\iPq\\h s + 
H^illtf 8 " 1 an d the mapping (fa, fa) i-> (ip, d t ip) : H^H 8 ^ 1 — > C([0, T], iiP) x 
C([0, T], if s_1 ) continuous. 



We say that the Cauchy problem eq. (6.1) is locally well-posed in H s 
for any s > | + 1. The proof (see e.g. [29]) of Theorem 6.2.1 relies on a 
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fixed point iteration together with energy estimates for the linear wave equa- 
tion, a Sobolev embedding, the so-called Moser inequality and the Gronwall 
inequality. We remark that Theorem 6.2.1 is equally valid for systems of 
equations, i.e. vector-valued ip and /. 



6.2.2 Wave Maps 

Observe that the wave maps equation is in fact a system of semilinear wave 
equations and the involved nonlinearity satisfies the requirements of Theo- 
rem [6T27T Thus, we immediately obtain local well-posedness of the Cauchy 



problem for wave maps in H s for any s > | + 1. At this point we should 
remark that wave maps have orginally been required to be C°° (cf. ch. |2j. 
However, from the point of view of partial differential equations this is an 
unnecessary restrictive assumption and therefore we relax it. 

It turns out that this local well-posedness result can be improved. The 
important observation in this respect is the fact that the nonlinearity in the 
wave maps system is not generic. It satisfies the so-called null condition 
(cf. [27]). Exploiting this special algebraic structure it is possible to show 
local well-posedness of the Cauchy problem for wave maps in H s for s > | 
and n > 2 (see [29] and references therein). This result is sharp in the sense 
that there exist wave map systems which are not locally well-posed (i.e. 
ill-posed) in H s for s < f , see [14]. 



6.3 Global Results 

In general, global existence of solutions for semilinear wave equations is ex- 
pected to hold only if the data satisfy certain "smallness" conditions. We 
cite some results in this direction for wave maps. On the other hand, if the 
solution does not exist for all times, the question arises how the breakdown 
occurs. For our particular wave map model this issue will be studied in the 
next chapter. 

There are many recent results concerning the global well-posedness for 
the Cauchy problem of wave maps and we are unable to mention them all 
(see e.g. [53], [32], [28] and references therein). However, of most interest for 
our purposes is Tao's work [51] which deals with wave maps from (n + 1)- 
dimensional Minkowski space to the (m— l)-sphere for m, n > 2. In this work 
it is shown that, roughly speaking, a wave map with smooth initial datum 
(ipo, ipi) which is small in the homogeneous Sobolev space H n l 2 x H n l 2 ~ x 
extends globally in time and stays smooth. Thus, for time evolutions starting 
from smooth data (^o,^i), blow up can only occur if HV'oll rW/2 + HV'illjU/a-i 
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is sufficiently large. 

Finally, we mention two earlier global existence results for wave maps 
from (3 + 1) Minkowski space to the three-sphere which have been obtained 
by Kovalyov [30] and Sideris [48]. We also refer to the work of Shatah and 
Tahvildar-Zadeh [47] for the special case of equivariant wave maps. 
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Chapter 7 

Self— Similar Solutions 



7.1 Blow Up 



We discuss solutions of the wave map problem eq. (2.3) with smooth initial 
data which become singular after a finite time. Such a behaviour is called 
blow up. We remark that this phenomenon can already be observed for 
ordinary differential equations. Consider for example the equation u' = u 2 
with initial data u(0) = 1/T. The solution is given by u{t) = (T — t)^ 1 and 
thus it ceases to exist at t = T. We have already mentioned regularity results 
which state that the solution is smooth for all times provided the initial data 
are smooth and small in some Sobolev space. Hence, blow up can only occur 
if the data are large enough. 

Technically we note that the derivations in this section have an informal 
character, i.e. we relax the mathematical rigor and restrict ourselves to a 
heuristic discussion. 



7.1.1 Scaling and Criticality Class 

We consider the conserved energy 



EM) := r (#(*,r) +#(f,r) + r 2 dr. (7.1) 



of the wave map equation 



^ + =^ = 0. (7.2) 



A mapping (t,r) (£/A, r/A) for a constant A > is called a dilation. 
Note that eq. (7.2) is invariant under dilations: Suppose ip solves eq. (7.2). 
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Then, ip\ denned by ip\{t,r) := ip{t/ A,r / A) is also a solution of eq. (7.2) 



provided that r in eq. (7.2) is substituted by r/A. This scale invariance 



can be used to classify conserved quantities of the equation. Let ip be a 



solution of eq. (7.2). The energy eq. (7.1) scales as E^ x (t) = X a E 1 p(t/X) 
for a = 1. One says that the scaling of energy is subcritical, critical or 
supercritical ifa<0,a = 0ora>0, respectively. Thus, the energy for the 



wave map equation eq. (7.2) is supercritical. An informal principle states 
that solutions of energy supercritical equations develop singularities for large 
initial data while they stay regular for small ones. Based on the criticality 
classification one gains a heuristic understanding of certain aspects of the 
dynamics. In the supercritical case it is favourable for solutions to shrink 
since this process is connected with a decrease of the local energy. Such a 
shrinking may eventually lead to singularity formation. Conversely, in the 
subcritical case shrinking is forbidden since it requires a larger and larger 



amount of energy. Hence, we expect eq. (7.2) to possess blow up solutions 



7.1.2 Characteristics and Finite Speed of Propagation 

We discuss a fundamental feature of wave equations: Finite speed of propa- 
gation of information. 

Characteristics It turns out that information propagates along certain 
curves in spacetime which are called characteristics. In what follows we will 
explain what is meant by this statement. Consider a first-order differential 
equation 

u t + Au r + f(u) = 

for a vector-valued function u = (ui, . . . , u n ) and a nxn-matrix A depending 
on t and r. We assume A to be diagonalizable with eigenvalues {Aj : % = 
1, . . . ,n} (the Aj's are functions of t and r as well). Hence, there exists an 
invertible matrix V (depending on t and r) such that V~ 1 AV is diagonal. 
Defining v := V~ l u we can write the system in component form 

d t Vi + Xid r Vi + g(v)i = (7.3) 

where g{y) := V^i^fiVv ) + AV r v + V t v). Now let r$ be a function of t such 
that fi(t) = Xi(t,ri(t)) for all t where ' := 4. Then, the spacetime curve 
t i — > (t,ri(t)) is called a characteristic. 



Finite speed of propagation Let v be a solution of eq. ( |7.3 ). We calcu- 



late the directional derivative of v along the characteristic t (t,ri(t)). 

^-Vi(t,ri(t)) = d t Vi(t,ri(t)) +d r v i (t,r i (t))f i (t) 
at 
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= -Ai(t,r,(t))9 r fi(t,ri(t)) +r i (t)d r v i (t,r i (t)) - g(v(t,ri(t)))i 

= -g(v(t,ri(t)))i. 

Integrating this equation yields 

Vi{t,ri(t)) = ^(0,^(0)) - / g(v(s,ri(s)))ids. 

Jo 

Thus, the solution at a spacetime point (t, r) depends solely on the value of 
the solution along the characteristics through that point. If fj(t) < oo for 
all % and t it follows that information encoded in the initial data propagates 
with finite speed. Hence, initial data given on a compact subset of the initial 
surface can only influence a compact spacetime domain in the future, the 
domain of influence. On the other hand, the value of the solution at a fixed 
spacetime point depends solely on the value of the field in a compact region 
of spacetime in the past, the domain of dependence. 

Note that the notion of domain of dependence is of fundamental im- 
portance for the numerical treatment of wave equations. Discretizing the 
equation in a way which is compatible with the characteristic structure is 
absolutely necessary in order to obtain a stable scheme. 



Characteristics of the Wave Map Equation We set U\ :- 

«3 := Vv and write eq. (7.2) in the form 



ip, u 2 := ipt 



where 



u t + Au r + f(u) = 



A 




and 



f(u)(t,r) 



-u 2 {t,r) 
u 3 (t,r) + sin( y r)) 




The eigenvalues of A are —1,1,0 and hence, there are three characteristics, 
r\{t) = ri(0) —t, r 2 {t) = r 2 (0) +t and r 3 (t) = r 3 (0). In a spacetime diagram 
the characteristics r\ and r 2 are lines with slope 1 or —1 and hence, the speed 
of propagation of information is limited by 1. 

Finally, we remark that the notion of characteristics can be defined in 
the more general context of a fully nonlinear first-order partial differential 
equation, see e.g. [19]. 
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7.1.3 Blow Up Solutions 



We intend to construct an explicit example of a solution of eq. (7.2) with 
smooth initial data which develops a singularity in finite time. The most 
promising strategy is to look for self-similar solutions. We have already 



discussed the dilation invariance of the wave map equation (7.2). Hence, it 
is natural to look for a solution which shares this invariance, i.e. we seek 
solutions ip with the property ip(Xt,Xr) = ip(t,r) for any A > 0. Such 
solutions are called self -similar. To this end we plug the ansatz ip(t, r) = 



f(fr-f) m ec l- (7-2) and obtain 



/" + V-^Vo (7.4) 
p p z {± — p z ' 



where p := ' '= j- and T > is an arbitrary constant. Shatah [45] 



showed that eq. (7.4) has a smooth solution. This solution has been found 



in closed form by Turok and Spergel [56] and is given by 

fo(p) '■= 2 arctan p. 

Set ipo(t, r ) := /o(t^*)- Then, ip is perfectly smooth for t < T but d r ip (t, 0) = 
(T — t)^ 1 and hence the spatial derivative at the center r = blows up for 



t — > T— . Thus, ipo is an explicit example of a solution of eq. (7.2) with 
smooth initial data Vo(0, •) and d t ipo(0,-) which develops a singularity in 
finite time. 

One might argue that this example is of no physical relevance since the 
solution ip is not a finite energy solution, i.e. E^ o (0) = oo. However, as we 
have seen in the previous section, the speed of propagation of information 
is limited by 1. This fact can be used to construct a blow up solution with 
finite energy. Consider smooth initial data which equal ipo(0, •) and dtipo(0, •) 
for < r < T and are identically zero for r > IT . These data have finite 
energy and due to finite speed of propagation the singularity at r = will 
form before any information of the region r > T reaches the center. Hence, 
the existence of self-similar solutions together with finite propagation speed 
implies the existence of a physically relevant blow up solution. 



7.1.4 Some Numerics 

A natural question is whether blow up occurs for generic initial data or the 
example given above is an "exceptional case" in a certain sense. In order to 
answer this question we employ some very simple numerics. 
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Coordinate transformation Since interesting things are expected to hap- 
pen around r = it is useful to make a coordinate change r i— > x := log(a+r), 
a > 0, in order to gain a better resolution near the center. The value of the 
constant a determines the quality of this primitive mesh refinement. The 



wave map equation (7.2) in these new coordinates reads 



Vfc - e 2x ip xx - e 2x - -i) x + l> = 0. (7.5) 



a \e d - — a \ 



Characteristics We calculate the characteristics of this equation. Setting 



U\ := tp, u 2 := ipt an d «3 := iftx we write eq. (7.5) in first-order form 

u t + Au r + f(u) = 

where 

A = 

The eigenvalues of A are — e~ x , e~ x , and hence the characteristic curves X\ 
and x 2 satisfy the differential equations xi(t) = —e~ Xl ^ and x 2 (t) = e~ X2 ^\ 
Thus, the characteristic speeds Xi and x 2 at the center x = log a are given 
by ±i and therefore they diverge as a approaches zero. 

Discretization We choose one of the simplest second-order accurate schemes 




to discretize eq. (7.5), the centered space and centered time method. We 
label discrete spacetime points by (n, k) where n, k are natural numbers in- 
cluding zero and use the following approximations for derivatives of if). 

4j tt {nAt,kAx) « — 



i/j x (nAt, kAx) 



Ax 2 

r k+1 - n-x 



2Ax 

where we have applied the us ual abbreviation ip^ := ip(nAt, kAx). Plugging 



these approximations in eq. (7.5) and solving for ip? we obtain 



r k +1 = a{kAx)At 2 rk+l ^|±jLl + b (kAx)At 2rk+1 Vk ~ x 



Ax 2 v ' 2Ax (7.6) 

n- 
k 



+c(kAx)At 2 sin(2^) + 2^ - ^ 



80 



7.1. BLOW UP 



where a(x) := -e 2x , b(x) := -e 2x $^, c(x) := (e ^ a)2 . 

We use a finite spatial grid which covers the interval [0, X] for some large 
X > and impose Dirichlet conditions ip(t,0) = and if)(t,X) = for all 



t at the endpoints. Since the discretization scheme (7.6) requires two time 
steps in the past we have to calculate the first time step by Taylor expansion. 
This yields the initialization formula 

At 2 ( , . .fk+i — 2fk + fk-i 



At 2 ( 

i>l = fk + 9kAt + — I a(kAx 



Ax 2 



+ b{kAx) h+l 2 J x k ~ l - c(kAx) sin(2/ fc ; 



where fk := ip(0, kAx) and := ^t(0, kAx) are the initial data. 



Domain of dependence As we have already mentioned, the character- 
istic speeds at the center are given by In order to obtain a stable 
discretization scheme it is necessary to make sure that the physical domain 
of dependence is included in the numerical domain of dependence. This 
requirement is known as the Courant-Friedrichs-Lewy condition (cf. [25]). 
Hence, we have to choose Ax, At and a in such a way that the condition 
^ < a is satisfied. We observe that in order to improve the quality of the 
mesh refinement (i.e. decrease a) we have to decrease At as well if Ax is kept 
fixed. Thus, one cannot increase the resolution around the center without 
increasing the computational effort. 



Numerical results We choose a Gaussian with amplitude A as initial 
data. This pulse splits into an ingoing and an outgoing wave packet. We 
make the grid large enough and focus on the ingoing pulse since nothing 



interesting happens to the outgoing one. Fig. |7.1| shows the time evolu- 
tion of a Gaufi pulse with small amplitude in the original (t, r)-coordinates. 
Then we successively increase the amplitude A and monitor the r-derivative 



of the solution ip at the center, given by a 1 ip x (t, log a). Fig. 7.2 shows 
maxj ip x (t, log a) of the solution ip plotted against A. 

We observe that the spatial derivative of the solution at the center di- 
verges when A approaches a critical value from below. A similar result can 
be produced when using other types of initial data. From this observation we 
conclude that the blow up is a generic phenomenon which occurs whenever 
the initial data are large enough. 
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Figure 7.1: Time evolution of a GauB pulse 



Figure 7.2: Blow up 



7.1.5 Universality of Blow Up 



Bizori et. al. [10] have studied eq. (7.2) numerically and based on their 



observations they have formulated some conjectures concerning the blow up. 
They claim that there exists a large open set of initial data which lead to 
blow up and, in addition, the asymptotic shape of the blow up solution 
approaches the Turok Spergel solution tp as t —>■ T— locally near the center 
r = 0. Hence, the self-similar blow up behaviour defined by the solution i/jq is 
conjectured to be universal in this sense. Furthermore, families of initial data 
depending on a parameter p which interpolate between dispersion and blow 
up have been studied. There exists a critical value p* of the parameter p such 
that initial data with p < p* lead to dispersion while data with p > p* blow 
up. Considering initial data which lie exactly at the boundary, i.e. p = p*, 
another self-similar solution ipi which plays the role of an "intermediate 
attractor" has been identified. This means that the solution approaches ipi 
locally around the center for a certain time and eventually disperses or blows 
up via ipo since p = p* exactly is numerically impossible. With the help of 
the code developed in the previous section it is possible to reproduce these 



results. Fig. 7.3 for example shows the last stages of the self-similar blow 



82 



7.2. PROPERTIES OF SELF-SIMILAR SOLUTIONS 



up of a solution if) with initial data of the form 

A sin r for < r < 



ifj(0,r) 



A exp( 



)2\ £ Qr r > 



and ipt(0, •) = 0. The dashed line is a plot of the Turok Spergel solution ip 
with appropriately chosen T. One sees that for small r the two solutions if) 
and if>o coincide. 




Figure 7.3: Self-similar blow up via ipo 



Hence, the simple equation (7.2) shows very interesting behaviour and a 



better (mathematically rigorous) 
sireable. 



understanding of these phenomena is de- 



7.2 Properties of Self— Similar Solutions 



Numerical studies of eq. (7.2) suggest that self-similar solutions play an 



important role in the dynamics of time evolution. Hence, it is necessary to 
take a closer look at the equation 



f"+-f- 



sin(2/) 
p 2 (l - p 2 



0. 



(7.7) 



7.2.1 Existence of Self— Similar Solutions 



Eq. (7.7) can be solved numerically using a shooting and matching technique. 



This has been done first by Aminneborg and Bergstrom [3] and later in [10] 
and [33]. The singular behaviour of eq. ( |7.7 ) at p = and p = 1 yields 
the regularity requirements f(0) = and f(l) = | for smooth solutions /. 
One imposes these boundary conditions and integrates the equation away 
from the singularities towards p = \ with a standard ODE integrator. By 
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this, one obtains two solutions f\ and f r on [0, |] and [|,1], respectively. 
Varying the free parameters //(0) and /^(l) one tries to smoothly match the 
two solutions at p = |. It turns out that there exists a countable family 
{/n : n = 0, 1, 2, . . . } of different self-similar solutions whose existence has 
been proved by Bizoh [6]. 



Theorem 7.2.1. 

eq. (7.7) satisfying the boundary conditions f n (0) 



There exists a countable family of smooth solutions f n of 

and f n {l) = f . The 
2, . . . denotes the number of intersections of f n with the line 



index n = 0, 1 
/ = § (the equator of S 3 ) on p G [0, 1). 



Additionally it is shown in [6] that, for n — > oo, the solutions f n converge 
to the limiting solution f^ = | pointwise for all p G (0, 1]. Furthermore, we 
note that the self-similar solutions f n can be interpreted as harmonic maps 
from the hyperbolic space H 3 to S 3 (cf. [13]). 



7.2.2 Stability Properties 

Stability of self-similar solutions is an important issue. Highly unstable 
solutions are not expected to play a role in the time evolution of generic 
initial data since they cannot be approached while stable solutions may act as 
attractors. In order to make these ideas more precise it is useful to transform 



the wave map equation (7.2) to the new space coordinate p = jF-: which 
yields 



2n_, 2(1 -p 2 )^ , sin(2^) 



(T - tyiPtt + 2p(T - t)Vt p - (1 - p')ip PP - v > p + = 0- ( 7 - 8 ) 

Hence, all self-similar solutions with blow up time T are static solutions of eq. 



(7.8). Based on the numerical observations of sec. 7.1 we expect the Turok 
Spergel solution ip to act as a static attractor for solutions of eq. (7.8). 
However, the situation is more delicate. Since the Turok Spergel solution is 
not a single solution but a family of solutions (depending on the parameter 



T), we expect a solution of eq. (7.8) with blow up initial data to converge 
to a certain i^q for an appropriate T. Hence, if one fixes T, not all blow 
up solutions will converge to ipo but only the ones with the "right" blow 
up time. Numerically this can be tested by considering a family of initial 
data depending on a parameter p. Then, it should be possible to adjust 
the parameter p such that the fine-tuned solution converges to ipo with a 
prescribed T. 
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7.2.3 Hyperbolic Coordinates 



Eq. (7.8) is not well-suited for a rigorous mathematical analysis since the 
coefficients depend on t. Furthermore, the mixed derivative is bothersome. 
Hence, we intend to introduce a new time coordinate a in order to simplify 
the structure of the equation. To this end we interpret the new coordinates 
a, p as functions of t, r and calculate the new coordinate vector fields d a , d p 
with the help of the equations d t = o- t d a + p t d p and d r = a r d a +p r d p . In order 
to avoid bothersome off-diagonal terms we require the new coordinates to be 
orthogonal, i.e. t](d a , d p ) = where n is the (coordinate representation of the) 
Minkowski metric. Setting p(t, r) = this yields the partial differential 
equation 

T-t 

a t o r = 

r 

which is a transport equation for the function a. The general solution is 
given by 

a(t,r) = g(t 2 - 2Tt - r 2 + c) 
where g is a free function and c an arbitrary constant. We transform eq. 



(7.8) to the new time coordinate a. By construction, the mixed derivative 
ip ap drops out. We use the remaining freedom to make the coefficients of 
the equation cr-independent. This yields c = T 2 and g = a log where a is a 
constant. Therefore, we obtain a(t, r) = a log((T — t) 2 — r 2 ). Since we want 
a to increase with t we choose a negative and it turns out that a = — | is 
convenient. Hence, we arrive at a = — log v/(T — t) 2 — r 2 . The wave map 



equation (7.8) transforms to 



- 2,v - (1 - m„ - + (1 -' 2) f ( ^ = 0. (7.9) 

p p 2 

We will refer to (a, p) as hyperbolic self-similar coordinates or simply hy- 
perbolic coordinates since the lines a = const are hyperbolae in a spacetime 
diagram. Note that a is only defined for r < T — t and hence, the hyperbolic 
coordinates cover the interior of the backward lightcone r = T — t of the blow 
up point {t,r) = (T, 0). The inverse transformation is given by 

1 p 

t = T ; and r = -j^=^= . 

- p 2 e a y / l-p 2 



Fig. TA shows lines of constant a in a spacetime diagram where T = 1. The 
dashed line is the backward lightcone of the blow up point. 
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Figure 7.4: Hyperbolic coordinates 



7.3 Well-Posedness 



We prove well-posedness of the Cauchy problem for the linearization of eq. 



(7.9) which governs the linear flow around a self-similar solution /„. 



7.3.1 Linearization 

Linearized flow around f n We substitute the ansatz ip(<j,p) = f n (p) + 



0(<r,p)ineq. (7.9) where f n is a self-similar solution of eq. (7.2). Expanding 
the nonlinear term in powers of and neglecting all terms of order higher 
than 1 yields the linear evolution equation 

A <M h ^ 2(1 -p 2 ) 2 ^ 2(l-p 2 )cos(2/ n ) ^ 

ffCT - 20 CT - (1 - p ) pp p H = (7.10) 

for the perturbation 0. The wave map /„ is said to be linearly stable if so- 
lutions of eq. (7.10) do not grow (with respect to some suitable norm) as a 
increases. However, we emphasize that our discussion is still on a heuristic 
level since we have not specified what we mean by a solution so far. Further- 
more, we note that it is by no means clear whether the nonlinear flow around 
f n can be approximated by this linear equation. There are explicit examples 
of nonlinear evolution equations where certain phenomena cannot be treated 
by linear perturbation theory, although the perturbations are small in a cer- 
tain sense (e.g. [9]). However, it is generally believed that instabilities in the 
linearized problem lead to instabilities in the nonlinear case and therefore it 
is useful to study the linearized equation. 



Symmetries We informally discuss the role of the time translation symme- 
try of eq. ( 7.2 ) . Let ip be a solution of the original wave map equation ( 7.2 ) in 
(t, r)-coordinates. Define a function ip £ by ip £ (t,r) := if)(t — e,r). Then, ip e is 
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also a solution of eq. (7.2). Hence, the mapping $ e : ip 

$ £ o$,5 and $" 



ip e maps solutions 
- Let ^ n be a 



to solutions and satisfies $ = id, 
self-similar solution with blow up time T given by ip n (t, r) = f n (jrb) ■ Then, 

&e(i>n) = i>n an d V ; n(^ r ) = fn { x+e-t ) • Thus, $ £ maps self-similar solutions 
with blow up time T to self-similar solutions with blow up time T' := T + e. 
The generator of the orbit {& £ (ip) : e G IR} is given by ^\ £= o& e (ip) and 
we readily calculate ^| e= o$ £ (V ; ) = — F° r a self-similar solution ^„ in 
(cr, p)-coordinates we obtain 



a,p) :-- 



d 

de 



$ e (^.)(a,p) = e> v /l-pV;(p). 



By direct calculation one easily verifies that 0^ solves eq. (7.10) and hence, 



the time translation symmetry of eq. (7.2) is reflected by an exponentially 



growing solution of the linearized equation (7.10). This exponential instabil- 
ity is referred to as the gauge instability. 

Although 0^ is not different iable at p = 1 (the backward lightcone of the 
singularity) we expect this solution to play a role in dynamical time evolution 
since p = 1 is not in the domain covered by the hyperbolic coordinates (only 
< p < 1 is valid, cf. fig. 7.4). This fact spoils the linear stability analysis 



to a certain degree since the best we can expect is a growth estimate like 
< Ce a \\(f)(0, -)|| for solutions of eq. (7.10) which would only rule out 



a. 



the existence of solutions that grow faster than the gauge instability. 



7.3.2 The Operator A 

Our aim is to give a rigorous operator formulation of the evolution problem 



eq. (7.10) and apply the semigroup theory developed in sec. 5.1 to show 



well-posedness. 



Simplifications First of all we make the simple transformation i— ► 
where 4>(cr,p) := e~°"0(cr, p) to get rid of the first order term — 2(p a in eq. 



(7.10). The transformed equation reads 



(1-P 2 ) 2 < 



pp 



2(WV 2 (1 -p3)cos(2/,,)-^ = a 

P p l 



We split the "potential" and write eq. (7.11) as 

2(1 -P 2 ) 2 



:i-p 2 )Vpp 



p 



2(1 - o 2 ) 2 ~ 

• P + ~ + 9n<P = (7. 12) 

p 2 
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where 

2(l-p 2 )cos(2/ n (p))-p 2 -2(l-p 2 ) 2 



9n(p) ■-- 



P 1 



Note that g n is regular at p — 0, i.e. g n G C[0, 1] which can easily be checked 
using de rHospital 's ru le. The idea now is to give a well-posed operator 



formulation of eq. (7.12) without g n (j) and to apply a perturbation argument. 



Thus, we first consider the formal differential expression a given by 



p 2 \ dp dp 

We set w(p) := nz^gp and define the Hilbert space H by H : = L^(0, 1). 

The method of Frobenius We give a brief description of a well-known 
method which provides us with asymptotic estimates for solutions of the 
equation au = 0. The method of Frobenius is a standard approach for ob- 
taining series expansions for solutions of linear second order ordinary differ- 
ential equations with meromorphic coefficients around singular points. Con- 
sider a differential equation u" + pu' + qu = in C where p and q are 
holomorphic functions except for a set of isolated points in the complex 
plane. Suppose that z = z$ is a regular singular point which means that 
p := \im z _> zo (z — z )p(z) and q := lim 2 ^ 20 (^ — z ) 2 q(z) exist. Then, there 
exist two linearly independent solutions wo, u± whose asymptotic behaviour 
for z — > zo can be stated explicitly: We denote the solutions of the so-called 
indicial equation s(s — 1) + pos + q$ = by s_ and s + (the indices) where 
|s_| < \s + \. 

• If the difference s + — s_ is not an integer then uq(z) = (z — z ) s+ u (z) 
and U\(z) = (z — z ) s -Ui(z) where Uq an d U\ are holomorphic around 
z = Zo with Uq(z ) 7^ and U\{zq) ^ 0. 

• If the difference is an integer then we have Uq(z) = (z — zq) s+ uo(z) 
and U\(z) = cuq(z) \og(z — z ) + (z — z$) s ~ui(z) where again u and u\ 
are holomorphic around z = z with Uo(z ) ^ and U\(zq) ^ 0. The 
constant c may also be zero and thus it is possible that the logarithmic 
term does not appear. 



In either case these representations are valid in the largest open circle around 
Zq which contains no other singularity of p or q. 
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Endpoint classification Consider the Sturm-Liouville problem au = 
on (0, 1). The endpoints p = and p = 1 are regular singular points and 
therefore, we can apply the method of Frobenius to obtain asymptotic esti- 
mates for solutions. Around p = the indices are —2 and 1. Thus, there is 
only one solution which belongs to H near p = and Sturm-Liouville theory 
(sec. 4.2) tells us that p = is in the limit-point case. Around p = 1 the 
indices are both equal to \ which shows that there does not exist a solution 
which belongs to H near p — 1. Thus, p = 1 is in the limit-point case as 
well. 

Definition of the operator A We set p(p) := p 2 and q(p) = 2. Then, 
au = reads 

— (-(pit')' + qu) = 0. 
We define V(A) := {u G H : G AC\ oc (0, 1), aw G #} and Au := au for 



w G X'(A). According to Lemma 4.2.11, the operator A : T>(A) C if — > H is 
self-adjoint. 

7.3.3 Properties of yl 

We claim that A : T>(v4) (Z H —> H satisfies the estimate (Au\u)h > 7(m|m)// 
for some 7 > and all u G T>(A). In order to show this we apply Hardy's 
inequality. 

Lemma 7.3.1 (Hardy's inequality). Let u G C x [a, b] and u(a) = 0. Then, 

dx < 4 / u x ) 2 dx. 



L {x-a) 2 
Similarly, if u(b) = then 

\u(x)\ 2 



dx < 4 / Iw'faOPofe. 



/a (x-b) 2 

Proof. Let m G C 1 [a, 6] with tt(a) = 0. Integration by parts yields 



\u(x)\ 2 ^_ \u{x)\ 2 '■ 



. u'(x)u(x) + u{x)u'(x) 
x — a 



a J a 



l a (x — a) 2 (x — a) 

Using de l'Hospital's rule and u(a) = we conclude that 

hm = 0. 

x^a+ (x — a) 
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Observe that — < and hence we have 

(6-o) — 

' v ; ' dx < 2 / ^^k/xcte 
{x — ay J a x — a 

by Cauchy-Schwarz. Dividing by f f Q 6 |"^| 2 dx ~\ and squaring the resulting 

inequality yields the claim. The same calculation can be applied for the case 
u{b) = 0. □ 

Lemma 7.3.2. There exists a 7 > such that the operator A : T>(A) C H — > 

H satisfies the estimate 

(Au\u)h > i(u\u)h 

for all u G V(A). 

Proof. Suppose u G 2? (A) has compact support. Then, u G C 1 [0, 1], u(0) = 
tt(l) = and integration by parts shows 



(Au\u) H = / p 2 \u(p)\ 2 dp + 2 \u(p)\ 2 dp 



Observe that 



1/2 i.1/2 

p 2 \u\p)\ 2 dp + 2 I \u(p)\ 2 dp>2 I \u{p)\ 2 dp> 



2 


/■1/2 




'0 




p 2 


(1 


-P 2 ) 



>C r -r J —^-Mp)\ 2 dP- 



We use the letter C for a generic real constant greater than zero which is not 
assumed to have the same value every time it appears. Furthermore, using 
Hardy's inequality we estimate 

p 2 \u'(p)\ 2 dp + 2 f \u{p)\ 2 dp>C f \u'{p)\ 2 dp 

1/2 Jl/2 A/2 



Adding up these two estimates we arrive at 



f p 2 \u\p)\ 2 dp + 2 f\u(p)\ 2 dp>C f \ P \^ \u{p)\ 2 dp = O »[»)//■ 
Jo Jo Jo l i — P 



Since {u G P(-A) : u has compact support} is a core for A (Corollary 4.2.1) 
this inequality is valid for all u G V(A) and we conclude that (Au\u)h > 
■j(u\u)h for a 7 > and all u G T>(A). □ 
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Well-posedness We define Y := V(A 1/2 ), X :=Y x H, V(L) := V(A) x 



5.1.2 



we 



V(A 1 ^ 2 ), L(u,v) := (v, — Au) for (u,v) G T>{L). Applying Theorem 
conclude that L generates a strongly continuous one-parameter semigroup 
S : [0, oo) — > B(X) on X satisfying ||S(0")||b(x) < 1. In particular, the 
Cauchy problem 

^u(cr) = Lu(ff) for a > 
u(0) = u 

for u : [0, oo) — > X with initial data u G X which is an operator formulation 
of the unperturbed equation 



2(1 -p 2 f~ 2(1 -p 2 f~ 



^2 



P P 

is well-posed. 



7.3.4 Bounded Perturbations 

Now we turn to the full linear equation 

(p aa - (1 - p ) (p pp <p p H - + g n <p = 0. 

p p A 

An operator formulation can be given as follows. We adopt the notation of 
the last section, i.e. 

-£u(a) = Lu(cr) for o > 
u(0) = u 

for u : [0, oo) — > X with initial data u G X is an operator formulation of the 
unperturbed equation (g n = 0). Now we introduce a perturbation operator 
L' n : X — > X defined by L' n (u,v) := (0, —g n u) for (u, u) G A. This definition 
makes sense since g>„ G C[0, 1] and hence, — g n w G if u G K C ff. We 
readily estimate ||L^u||x < ||Pti]|c[o,i] ll u llx for all u G A. Thus, L' n is a 
bounded operator on A and ||L^||b(x) < ||<7n||c[o,i] f° r an Y n — 0, 1,2, ... . 
Therefore, 

£u(a) = Lu(ff) + L>(o-) for a > , ? ^ 



u(0) = u 

for u : [0, oo) — ► A with initial data u G A is an operator formulation of 



the full linearized problem eq. (7.12). 



We apply the bounded perturbation theorem 5.2.1 to the operator L + L' n 
defined above which shows that L + L' n generates a strongly continuous one- 
parameter semigroup S n : [0, oo) — > B(X) on A satisfying 

\\S n (a)\\ B{x) < (7.14) 
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for all a > 0. It follows that the Cauchy problem eq. (7.13) describing the 



linearized flow around a self-similar solution f n in hyperbolic coordinates is 
well-posed. 



We note that the above estimate eq. (7.14) is not very satisfactory because 

C[o,i] = 15 which yields 



for the Turok Spergel solution / we have \\go\ 



\So(<r) 



\B(X) 



< e 



15a 



while the original intention was to rule out solutions that grow faster than the 
gauge instability. Translated to the semigroup approach this would require 
the estimate HSo^llspr) — 1 (remember the transformation <p i— > <p) for the 
semigroup S generated by L + L' . In order to achieve this we have to study 
the spectrum of the operator L n := L + L' n in more detail which is the topic 
of the next chapter. 
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Chapter 8 

The Spectrum of Lq 



We study in detail the spectrum of the operator L which is the generator 
of the semigroup describing the linearized flow around the Turok Spergel 
solution fo(p) = 2arctan(p). This will lead to a significant refinement of the 
growth estimate 7.14 and eventually we will be able to derive a result which 
rules out the existence of solutions of the linearized equation (7.10) that grow 
faster than the gauge instability. 



8.1 The Operators A n and L n 



We adopt the notation of sec. 7.3, i.e. H := L^(0, 1) with w(p) := nf p 2\2 , 
p(p) := p 2 , q(p) = 2, au := ^-(—{pu')' + qu) and T>(A) := {u G H : u,pu' G 
^4Cioc(0, l),aw G H}, Au := au. First we discuss properties of the operator 
A n defined by V{A n ) := V(A) and 

A n u := Au + g n u 

for u G T>(A n ) and 

2(l-p 2 )cos(2/ n (p))-p 2 -2(l-p 2 ) 2 



9n(p) ■-- 



P 2 



where f n is the n— th self-similar wave map (cf. Theorem 7.2.1). Hence, 
L n : T>(L) C X — > X is given by L n (u,v) = (v, — A n u) for (it, v) G T>(L n ) = 
V{A)xV{A 1 / 2 ) whereX := V(A 1/2 ) x with ||(«,v)||^ := || A 1 l 2 u\\ 2 H + \\v\\ 2 H 
for (n, f ) G X. Thus, A n is the essential nontrivial part of L n . 

First of all we note that the operator A n is self-adjoint. This can either 
be shown directly via Sturm-Liouville theory or it follows from the fact that 
A is self-adjoint, the boundedness of the symmetric multiplication operator 
u I— > g n u on H and the following theorem. 
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Theorem 8.1.1. Let A : T>(A) C H — * H be a self-adjoint operator on a 
Hilbert space H . If B : H — > H is bounded and symmetric then A + B is 
self-adjoint. 

Proof. See [26], p. 287. □ 

Hence, we see that the spectrum of A n is real. The following lemma shows 
how the spectra of A n and L n are related. 

Lemma 8.1.1. Let A G C. Then, A G cr(L n ) if and only if —A 2 G cr(A n ). 



Proof. Invoking Lemma 5.1.6 we conclude that a(L n ) C {A G C : —A 2 G 
a(A n )}. Hence, it remains to show that —A 2 G o-(A n ) implies A G o~(L n ) 
which is equivalent to A G p(L n ) =>- —A 2 G p(A n ). 

Suppose A G p(L n ), i.e. (A — L n )~ l : X — > X exists as a bounded operator 
on X. Let f E H and define (u,v) : = (A - L n ) _1 (0, -/). Then, (0, -/) = 
(A — L n ) (u, v) = (Aw — v, A n u + Xv) and hence, (—A 2 — A n )u = f which shows 
that (-A 2 - An) : V(A) C H -> if is surjective. Now let u G with 
(-A 2 - A n )u = 0. We have (A - L n )(u, Xu) = (Xu - Xu, X 2 u + A n u) = (0,0). 
However, since A — L n is injective, we conclude that u = which shows that 
— A 2 — A n is injective as well. Thus, we infer —A 2 G p(A n ). □ 

The above result shows that the spectrum of L n can be calculated from 
the spectrum of A n and therefore, it suffices to study the operator A n . Fur- 
thermore, since —A 2 G R if A G a(L n ), we observe that a(L n ) is a subset of 
the union of the real and imaginary axis. 



8.2 The Spectrum of A 

We explicitly calculate the spectrum of A . 
8.2.1 Prerequisites 

A technical lemma In what follows we will frequently encounter two types 
of singular behaviour. Terms of the form f(x) J* g(s)ds where / becomes 
unbounded for x — > and expressions like f(x) g(s)ds where the integral 
is divergent for x —>■ while / goes to zero. In either case the singular 
behaviour of one factor might be compensated by the other one. Note that 
we cannot treat these problems with de P Hospital's rule since normally the 
function g belongs to some Lebesgue space only and hence it cannot be 
evaluated at single points. The following technical lemma shows how to deal 
with such problems. 



8.2. THE SPECTRUM OF A n 95 



Lemma 8.2.1. Let f E L 2 (0, 1). 

1. If a > — \ then u, defined by u{x) := x~ a ~ x f£ s a f(s)ds, belongs to 
L 2 (0,1). 

2. If a > \ then u, defined by u(x) := x a ~ l J 1 s~ a f(s)ds, belongs to 
L 2 (0,1). 

Proof. Let / E L 2 (0, 1). 

1. We define v(x) := J* s a f(s)ds where a > — |. Using integration by 
parts and Cauchy's inequality with e we readily calculate 



i ri 

2j / ™-2a-2 



dx = J x a \v(x)\ dx 

1 







-2a - 1 



-2a-l|../'„M2| 1 / „-2a-l! 



u(x)| 2 '- / x- 2a - 1 Re[u / (x)u(a:)]da; 



o 

■l 

< Cx- 2q -V(x)| 2 |J + C / ^""-^(a;)! |a;-V(a;)|d2; 

Jo 

< C x- 2a - 1 \v(x)\ 2 \ 1 + 1 / x~ 2a ~ 2 \v(x)\ 2 dx + C [ x 2a \v'(x)\ 2 dx 

= Cx- 2a -Xa;)| 2 |J + ^ / \u(x)\ 2 dx + C [ \f(x)\ 2 dx 

2 Jo Jo 

and hence, we arrive at the inequality 

ll«lll w) <c'||/||! 2( o il) + c'x- 2 «- 1 Kx)| 2 |;. 

However, using the Cauchy-Schwarz inequality we estimate 

k*)i 2 < ri/(^)i 2 ^<^ 2a+i ii/iii 2( o,i) 

JO Jo 

which implies 

lim x' 2a ' 1 \v(x)\ 2 < oo 

x—*0 

and the claim follows. 

2. Define v(x) := J s~ a f(s)ds. For a > | the very same calculation as 
above can be applied. For a = \ one encounters logarithmic terms but 
nevertheless the same reasoning goes through. 

□ 



Remark 8.2.1. Of course, the choice of the interval (0, 1) in Lemma 8.2.1 



is 



completely arbitrary and therefore analogous results are true for a general 
finite interval (a, b) at either endpoint a or b. 
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u(x) 
v(x) 



The ~ notation In connection with asymptotic estimates it is useful to 

\a,b) — > C and c G [a, b]. We write 

< oo j^J Note that due to the usage 

is not required to converge for x — > c. For 

instance we have sin - ~ 1 for x — > 0. Simpler examples are e x ~ 1 and 
sinx ~ x for x — > 0. But note carefully that according to our definition we 
also have sinx ~ 1 for x — > 0. 



introduce a common notation. Let u : 
u(x) ~ v(x) for x — > c if limsup 2 

of lim sup in the definition, ^£4 



8.2.2 The Point Spectrum of A 

Singular points We consider perturbations around the Turok Spergel so- 
lution fo(p) = 2arctan(p). The function g Q is given by 

2(l-p 2 )cos(2/ (p))-p 2 -2(l-p 2 ) 2 
9o{P) = 5 • 

p 2 

However, the term cos(4 arctan(p)) can be written as a rational function 

1 - 4p 2 + p 4 



cos(4 arctan(p)) 



;i + p 2 )^ 



which reveals two singularities p = ±z in the complex plane. Thus, go is 
meromorphic and the equation (A — Aq)u = has six regular singular points 
p = 0, ±l,±z,oo. With the help of the transformation p i— > z := p 2 this 
number can be reduced by two. Hence, the general (formal) n solution of 
(A — Ao)u = can be given in terms of Heun's functions which are the 
solutions of the general Fuchsian equation with four regular singular points 
(cf. [43]). However, since the study of Heun's functions relies heavily on 
numerical techniques, this observation is not very useful for us at the present 
stage. 



Asymptotic estimates Nevertheless we can apply Frobenius' method to 
obtain asymptotic estimates of (formal) solutions of (A — A )u = around 
p = and p = 1. The equation (A — A )u = reads 

1 + A _ 2cos(2/q) \ 
(1-p 2 ) 2 p\\-p>)) U 

1 u(x) ~ v(x) is sometimes denoted as u — 0(v) where O is a so-called Landau symbol. 

2 We remark that the term "formal" in this context refers to the fact that an actual 
solution must belong to T>(Aq). We do not consider so-called formal power series solutions, 
i.e. series solutions whose convergence radius is zero. 



9 

u H — u 
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and thus, the indices at p = and p = 1 are {—2, 1} and {— ij—, 1+ ^ }, 
respectively. It is clear that the solution which behaves as p~ 2 for p — > 
does not belong to ZZ and hence, eigenfunctions u are holomorphic around 
p = and satisfy w(0) = 0. Furthermore, the requirement u £ H yields 
lim^i u(p) = 0. To conclude, we have the result that every eigenfunction u 
of A is in C 1 (0, 1) and satisfies u(0) = u(l) = 0. 

From the asymptotic estimates around p = 1 we immediately infer that 
there are no eigenvalues A > since neither of the two solutions is in H 
for A > because Re 1 " 1 ^ 3 ^ = Re 1+ ^~^ = |. Hence, we have shown that 
(T p {Ao) C (-oo,0). 



Nonexistence of eigenvalues smaller than We apply a standard os- 
cillation argument to show that there are no negative eigenvalues of Aq. The 
proof is based on the observation that the gauge mode 9, defined by 

1 + P Z 

(formally) satisfies Aq6 = as already discussed in connection with the gauge 
instability (cf. sec. 7.3). However, note that 6 ^ H and hence, it is only a 
formal solution. We remark that 9 has no zeros in (0, 1) which will be crucial 
now. 

Lemma 8.2.2. The equation (A — A )u = has no nontrivial solutions for 
A < 0. 

Proof. Suppose u £ V(A ) is a nontrivial solution of (A — A )u = for A < 0. 
It follows that u(p) ~ p for p — ► and u(p) ~ (1 — p) a for p — ► 1 where 
a := > | by Frobenius' method. The function u satisfies the equation 

„ 2(1 -p 2 ) 2 , 2(l-p 2 )cos(2/o)-p 2 

Am = —(1 — p ) u u H ^ u 

p p 2 

2 

for p £ (0, 1). Let [a, 6] C (0, 1). We multiply by nzzffi ®{p) an d integrate by 
parts to obtain 

Ja (1-P 2 ) 2 

where 0) := — u'9 denotes the Wronskian of the functions u and 9. 
For a^Owe obtain 

A t n ^ 2 , 2 u(p)9(p)dp = b 2 W(u,9)(b). 
Jo (l-P 2 ) 2 
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Note that the limit b — > 1 of the left-hand side exists thanks to the asymptotic 
behaviour u(p)9(p) ~ (1 — p) 1 / 2+Q and | + a > 1. Without loss of generality 
we assume u'(0) > 0. Suppose u has its first zero at b £ (0, 1]. It follows that 
u'(b) < and 

A f h ^ 2, 2 <P) 9 (p) d P = ~b 2 u'(b)9(b) 
Jo I 1 - P ) 

which is a contradiction since the left-hand side is strictly negative (A < 0!) 
while the right-hand side is either zero or positive. Hence, u does not have 
a zero in the interval (0, 1] which is a contradiction to u(p) ~ (1 — p) a for 
p — > 1 and a > |. □ 

Remark 8.2.2. We note that this observation has first been made by Bizoh 
[10]. 

Combining this result with a p (A ) C (— oo, 0) we conclude that there are 
no eigenf unctions at all and hence, the point spectrum of Aq is empty! We 
formulate this result proposition. 

Proposition 8.2.1. The point spectrum a p (A ) of A is empty. 



8.2.3 Invert ibility of A - A for A < 

Next, we study the continuous spectrum of Aq, i.e. we consider the inhomo- 
geneous equation (A — A )u = f for a given / £ H. We are interested in 
solutions defined on the open interval (0, 1) and fix A < 0. Let Uq and u\ de- 
note nontrivial (formal) solutions of the homogeneous equation (A — A )u = 
with the asymptotic behaviour Uo(p) ~ p for p — > and Ui(p) ~ (1 — p) a 
for p — > 1 where a = > \- The solutions u , U\ exist by Frobenius' 

method and they are both defined on (0, 1) since the minimal distance from 
p = (resp. p = 1) to the next singularity of the equation in the com- 
plex plane is 1. The two solutions uq and u\ are linearly independent as the 
following lemma shows. 

Lemma 8.2.3. The functions Uq and U\ are linearly independent. 



Proof. Suppose Uq and U\ are linearly dependent. Then, the asymptotic 
behaviour of U\ is «i(p) ~ p for p — > which means that U\ £ V(A ). Thus, 



u\ is an eigenfunction of Aq which is a contradiction to Prop. 8.2.1 which 



states that a p (A ) = 0. □ 
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By normalization we can always have [uq, U\\ p (x) = 1 which will be as- 
sumed from now on. 

We apply the variation of constants formula to obtain a formal solution 
of (A — Aq)u = f for / G H . The general solution can be written as 

u{p) = c u (p) + c 1 u x {p) - u (p) I ui(g) f(£)w(£)d£ 

+ u 1 (p) [ P Uo(£)f(Z)w(£)d£ (8.1) 
J pi 

for p G (0, 1) where Co, C\ G C and po, Pi £ (0, 1) are free constants. The 
so-defined u satisfies u,pu' G v4Ci oc (0, 1) and by differentiation one easily 
verifies that it (formally) satisfies (A — A )u = f. Thus, the question is 
whether the constants c ,Ci, p , pi can be adjusted in such a way that the 
resulting u belongs to H. The answer is yes as the following lemma shows. 

Lemma 8.2.4. For a given f G H define 

u{ P ) :=u ( P ) I «i(£)/(£M0# + ui(p) r«o(£)/(£M£K 

J p Jo 
for p G (0, 1). Then, u G V(A ) and (A - A )u = f. 

Proof. With the help of the asymptotic estimates wo(p) ~ P, u i(p) ~ P~ 2 f° r 
p — > and Uq{p) ~ (1 — p)' 3 , U\{p) ~ (1— p) a for p — > 1 where a := > - 



2 ^2' 

^ := < | together with Lemma 8.2.1 one easily observes that y/wu 

belongs to L 2 (0,1). □ 

Hence, for any A < and / G H we can explicitly construct a function 
u G ^(Aq) such that (A — A )u = f which shows that the operator A — A is 
surjective. Combining this result with o- p (A ) = we arrive at the following 
proposition. 

Proposition 8.2.2. The set (— oo, 0) is contained in the resolvent set p{A ) 
of the operator A , i.e. (A — A ) _1 G B(H) exists for A < 0. 



8.2.4 The Operator A - A for A > 

Finally, we study invertibility of the operator A — A for A > 0. We fix A > 
and consider the homogeneous equation (A — A )u = 0. The method of 
Frobenius implies the existence of two linearly independent solutions around 
p = 1 which both behave as (1 — p) 1 ^ 2 for p — > 1. Furthermore, around p = 
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we have a solution which behaves as p for p — > 0. Thus, we can find two 
linearly independent solutions u and U\ on (0, 1) which satisfy u (p) ~ p for 
p — > and ~ (1 — p) 1 / 2 for p — > 1. Again, by normalization we can 

assume that [uq, ui] p (x) = 1. According to the variation of constants formula 



eq. (8.1 ), solutions u of (A — A )u = f for given f G H have the form 



rp 

U(p) = C M (P) + CiMi(p) - M (p) / ttl(0/(0^(0^ 



Po 

+«i(p) f uo(o/(o^(Ode 

for p G (0,1). Now we try to adjust the constants c ,Ci,p ,pi in such a 
way that the resulting u belongs to V(A ). Since «i(p) ~ p~ 2 for p — > 
we are forced to choose C\ = — ^ Uo(0f(0 w (0d£ to compensate this bad 
behaviour. Hence, we have 

u( P ) = c Q u (p) - u {p) r«i(o/(o^(e)^+«i(p) f P uoiomMOdc 

J po JO 

(8.2) 

For any choice of po G (0, 1) we can certainly find an / G H satisfying 
UiitmZMSdt = c and [ Uo(OmMm 0. 

po Jo 

For such an / we have u(p) ~ U\{p) for p — > 1 since the first two terms in eq. 



(8.2) cancel in the limit p — > 1. Obviously, u (fc H and hence there exists an 
/ G H which is not in the range of A — A . Thus, the operator A — A is not 
surjective for A > 0. This shows that every A > belongs to the continuous 
spectrum of A since self-adjoint operators do not have residual spectra and 
hence, we have cr c (A ) = [0, oo). 

This result concludes our discussion of the spectrum of Aq and eventually 
we arrive at the following theorem. 

Theorem 8.2.1. The spectrum of the operator A is given by ct(Aq) = 
a c (A ) = [0,oo). 



8.3 Improved Growth Estimate 

With the results of the previous sections at hand we are able to improve the 
insufficient growth estimate 



\\S (a)\\ Bi x) < e 
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derived in sec. 7.3 for the semigroup So governing the linearized flow around 
the Turok Spergel solution in hyperbolic coordinates. 

Since inf a(A ) = we conclude that (A u\u)h > for all u G V(A ) 



(Lemma 4.1.1). Thus, Theorem 5.2.2 tells us that Lq generates a strongly 



continuous one-parameter semigroup So : [0, oo) — > B(X) satisfying So(a) < 
a for all ex > 1. 

Consider the conserved energy given by ||(u(cr), t>(cr)||# := ||Aj!/ 2 w(cr) \\ 2 H + 
\\v(a)Wfj for a classical solution (u(a), v (a)) := So(a)(uo,v ) with initial data 
(uo,vo) G V(L ). Since ^ ct p (A ) we see that there does not exist a u G 

1 /2 

V(A ), u ^ with A u = 0. Hence, the energy is a norm on T>(L) and we 
have || (u(a), v(a)) \\e = \\(uo,Vq)\\e for a classical solution with initial data 
(uo,vo) G V(L ). Actually we are interested in the function e a (u(a),v(a)) 
(recall the rescaling <j) \— > <fi) and hence, we conclude that there does not exist 
a classical solution which grows faster than the gauge instability. According 



to the discussion in sec. 7.3 this is the best result we could have expected. 
It shows that the Turok Spergel solution is as stable as it can possibly be 
in these coordinates. Hence, this result strongly supports the conjecture of 
linear stability of /o- 



8.4 Discussion 

We give a physical interpretation of the spectral behaviour of the operator 
/-,»• 



8.4.1 The Spectrum of L c 



According to Lemma 8.1.1, the spectrum a(L ) of L is given by a(L ) = 
{A G C : —A 2 G <j(Aq)}. Thus, since ct(Aq) = [0,oo), it follows that a(L ) = 
iM., i.e. the spectrum of L fills the imaginary axis. However, we emphasize 
that the point spectrum a p (L ) is empty and hence, there are no eigenvalues. 
Nevertheless, this spectrum is expected to lead to an exponential growth of 
solutions of the equation 



:i-^- 2 ^0, + 2(1 -^ COS(2/ °V = O (8.3) 



,^2, 

'•>•• 

p p z 



which describes the time evolution of linear perturbations of the Turok 
Spergel solution (remember the rescaling 0(a, p) = e _fT 0(cr, p) where L drives 
the evolution of 0). We have already discussed that the gauge instability is 
supposed to lead to solutions that grow like e a for o — > oo. However, this 
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instability cannot be responsible for the whole spectrum. Hence, the spectral 
behaviour of L calls for a physical (or more intuitive) explanation. 

8.4.2 Physical Explanation 

It turns out that the structure of the spectrum of L is closely related to 
the nature of the coordinate system (cr, p). A heuristic explanation of the 
instabilities caused by the spectrum of Lq can be given by the following 
argument. Consider a perturbation which has the form of a Gaussian. The 
time evolution will lead to outgoing wave packets which travel from smaller 
p to larger p. However, they cannot leave the backward lightcone of the 
singularity since the hyperbolic coordinates break down at p = 1. Hence, 
the wave will eventually cumulate near p = 1 and one observes exponential 
growth of the solution. Thus, this seemingly unstable behaviour is due to a 
defect of the coordinate system (a, p). 




Figure 8.1: Outgoing wave packet and hyperbolic coordinates 



Fig. 8.1 illustrates this phenomenon. The hyperbolic lines are sections 



a = const. The shaded region represents an outgoing wave which eventually 
leaves the backward lightcone (the dashed line). In (a, p) coordinates the 
wave comes closer and closer to p = 1 but never reaches it. 



8.4.3 Numerical Verification 

We employ some simple numerics to integrate the equation 



h « ^-P 2 ) 2 ! , 2(l-p 2 )cos(2/o)-p 2 ; - 

- (1 - p ) (j) pp (j) p + ^ <P = 



J.\2, 

'o„ 

P p' 



in order to illustrate the behaviour of solutions of the perturbation equation 



(8.3). We discretize the equation using the same scheme as in sec. 7.1 
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The characteristic speeds are given by ±(1 — p 2 ) and hence, they attain 
their maximal absolute value 1 at p = which implies that Act = 0.9Ap is 
sufficient to satisfy the CFL-condition. 




Figure 8.2: Time evolution of a GauB pulse 




Figure 8.3: Plot of o 



a, -)\\h against a 



Fig. 8.2 shows the time evolution of a Gaufi pulse. As expected, the wave 
packets slow down a nd ev entually freeze as they approach the backward light- 
1. Fig. " 



cone at p 
same time evolution 



8.3 



shows a plot of the function a \— > \\<p(a, -)\\h of the 
For late times the norm behaves like a constant and 
hence, the corresponding <j)(a, p) = e a <f)(<j, p) grows exponentially. This insta- 
bility is caused by the continuous spectrum of the operator L as discussed 
in the previous section. 
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Chapter 9 

The Functional Calculus 



For the sake of completeness we discuss another, completely different method 
for the treatment of ordinary differential equations on Banach spaces. This 
method is the standard operator theoretic approach in quantum mechanics 
for studying dynamics of the Schrodinger equation. It relies on the concept 
of functions of self-adjoint operators, i.e. the functional calculus provides 
a method to define an operator f(A) for a complex-valued function / and 
a self-adjoint operator A. Thus, solutions of ordinary differential equations 
on Banach spaces can be given in "explicit" form as functions of self-adjoint 
operators. However, this approach is not an equivalent substitute to semi- 
group theory since it is by construction restricted to self-adjoint operators 
while the notion of self-adjointness does not play a role in the formulation 
of semigroup theory. However, the examples we have studied so far have 
been formulated in a self-adjoint manner and hence, the functional calculus 
provides an alternative approach to these problems. 

9.1 The Spectral Theorem 

The material outlined in this section can be found in much more detail in 
standard textbooks, e.g. [59]. 

9.1.1 Spectral Families, Measures 

Spectral Families Let if be a Hilbert space and consider a family {E(\) G 
13(H) : A G K} of orthogonal projections E(X) depending on a real parameter 
A. Hence, for each A G R, E(\) : H — > H is a linear bounded self-adjoint 
operator on H which satisfies E(X) 2 = E(X). The family {E(X) G B(H) : 
A G R} is called a spectral family or resolution of the identity if it satisfies 
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• E(X 1 )E(X 2 ) = E(X 2 )E(X 1 ) = £(min{Ai, A 2 }) for all Ai, A 2 G R 

• lirriA^.oo E(X)u = and lim^oo E(X)u = u for all u G H 

• lim e ^ 0+ E(X + e)u = E(X)u for all u G H 



Measures Once there is given a spectral family it is possible to construct 
certain measures. Fix u G H and define F U (X) := (E(X)u\u)h- Then, F u : 
R — > R is monotonically increasing and right continuous, i.e. F u (Ai) < 
F n (A 2 ) if Ai < A 2 and lim £ ^ 0+ F U (X + e) = F U (X). Hence, F u is a distribution 
function in the sense of measure theory. Given a distribution function F u one 
can show that there exists a unique measure \l u : B — > [0, oo] which satisfies 
/4t((Ai, A 2 ]) = i ? u (A 2 ) — -^(Ai) for Ai < A 2 where B denotes the Borel cr- 
algebra on R. Note that /x u (R) = ||w||^ < oo and hence, the measure is 
finite. Moreover, observe that due to the special properties of E(X) we have 
A 2 ]) = || [E(X 2 ) - EiX^uWl for X, < A 2 . 
For fixed u, v G H we define another, complex-valued finite measure 
fi UjV : B -> C by 

^u,v{B) := | (n u+v (B) - /i u - v (B) + iji u _ iv (B) - ifi u+iv (B)) 

for all B E B. Invoking the polarization identity it follows that /i Uj „((Ai, A 2 ]) = 
(E(A 2 )«|«)h - (^(Ai)«|«)h for Ai < A 2 . 



9.1.2 Operators Defined via Measures 

Projections Fix u £ H, B £ B and consider the mapping F u> b '■ H — > C 
defined by F u>B (v) := ix UtV {B). First, we claim that F UiB (o;f) = 
for any a G C. Since the measure fi UjV is uniquely determined by its values 
on half-open intervals it suffices to show properties for B := (Ai,A 2 ] where 
Ai < A 2 and they automatically remain true for arbitrary B G B. Hence, the 
claim follows immediately from the formula // u ,u((Ai, A 2 ]) = (E(\2)u\v)g — 
(E(Xi)u\v)h- Furthermore, we have F U) b(v + w) = F U) b(v) + F Ui b{w) for 
v,w G H and \F UtB {v)\ < \\u\\h\\v\\h- Therefore, F uB '■ H — > C is a bounded 
linear functional on H. Invoking the Riesz representation theorem 3.1.1 we 
conclude that there exists a unique w G H such that F u> b(v) = (v\w)h for 
all v G H. Hence, for any B G B there exists a well-defined unique mapping 
P{B) : H -> H such that ^b(v) = (v\P(B)u) H for all u,v £ H. It is easy 
to see that P(B) is linear and bounded. Note that by construction we have 
A 2 ]) = E(X 2 ) — E(Xi) and hence, we have extended these projections 
to arbitrary Borel sets. Another way to write P(B) is (P(B)u\v)h = fJ>u,v(B) 
for all u,v G H. 
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Properties of P Note that (u\P(B)v) H = (P(B)v\u) H = fJ> v ,u(B) = 
H UtV (B) = (P(B)u\v)h by definition of /i UjV and hence, P(B) is self-adjoint. 
Furthermore, {P{B)P{C)u\v) H = ^ P{C ) u ,v{ B ) = Vu,v( B n C ) = ( P ( B n 
C)u\v)h for all v G H and B,C half-open intervals. By uniqueness of the 
measure we have P(B)P(C) = P(B n C) for arbitrary B G B and half-open 
intervals C. However, using the self-adjointness of P(B) we can interchange 
the role of B and C and hence, P{B)P{C) = P(BnC) holds for all B,C E B. 
In particular it follows that P(B) 2 = P(B) and thus, P(B) is an orthogonal 
projection. 



The domain of Af With the help of the measures fi u and fj, u>v we can 
define an operator Af : T>{Af) C H — > if for a given measurable function 
/ : K. — > C in the following way. Set 



V(A f ) :=^ueH: J \f\ 2 dfi u < ooj 



Then, u G f(^4/) implies au G P(A^) for any a G C since by definition we 
have // au ((Ai, A 2 ]) = \a\ 2 fi u ((\i, A 2 ]) for Ai < A 2 and, by uniqueness of the 
measure fi u , it follows that fi au = \a\ 2 /j, u . Similarly, for u,v G H we have 
< 2(fi u + fi v ) by the triangle inequality and hence, u, v G X> (A/) implies 
■u + v G X>(A/) which shows that T>(Aj) is a subspace of H. 

Furthermore, T>(Af) is dense in H which can be seen as follows. Define 
B n := oo, n]) for n G N. Since / is measurable, B n G B for all n G N. 

Let u & H and set u n := P(B n )u. Then we have 



2|| l|2 ^ 

"u < oo 



since // U „(C) = /ip(B„) u (C) = n C) for all C G Therefore, it n G 

T)(Af) for all n G N and one easily sees that w n — >■ -u in iJ. 

The operator A/ Fix u G f(^4/) and consider the mapping F : H — > C 
defined by 

One readily observes that F is a bounded linear functional on H and hence, 
by the Riesz representation theorem there exists a well-defined operator Af : 
T>(Af) C H — > H such that F(t>) = (Afu\v) H for all v £ H. It is easily seen 
that A f is linear and symbolically one usually writes 



A f = I /(A)dE(A). 
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9.1.3 Spectral Families and Self— Adjoint Operators 

The spectral theorem It turns out that to every self-adjoint operator 
there is associated a unique spectral family. This result is known as the 
spectral theorem. More precise we have the following statement. 

Theorem 9.1.1. Let A : T>(A) C H — > H be a self-adjoint operator on a 
Hilbert space H . Then, there exists a unique spectral family {E(X) : A G IR} 
such that 

A = [ XdE(X). 
Jr 

Representations All the information concerning the spectrum of a self- 
adjoint operator A : T>(A) C H — > H is encoded in its spectral family 
{E(X) :XeR}. 

To see this we note that for any spectral family {E(X) : A G IR} we have 
the representation 

I=\ dE(X) 



for the identity operator / : H —>■ H on H. This can be immediately seen by 



inserting the definitions. Consider the operator A\ (cf. sec. 9.1.2). Then we 
have T>(Ai) = H and (A\u\v)n = J dfi UjV = fi u>v (M) = (u\v)g for all u,v G H 
and hence, A\ = I. Thus, the operator Ao — A is given by 

A - A = [ (A - X)dE(X) 
Jr 

where {E(X) : X G IR} is the unique spectral family associated to A. Fur- 
thermore, one can show that 

||(A -A)u||^= / (A - X) 2 dfi u (X) 



where \x u is the measure associated to {E(X)} defined in sec. 9.1.1 



The spectrum of A We have Ao G o~ p (A) if and only if there exists a 
u G T>(A), u^O such that 

= \\(X - A)u\\ 2 H = [ (A -A) 2 dMA). 

Jr 

We conclude that the distribution function F u given by F U (X) = (E(X)u\u)h 
has to be constant except for a possible discontinuity at A = Ao- The re- 
quirements limA-»-oo F U (X) = and limA-«x) F U (X) = \\u\ljj together with 
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right continuity of F u imply F U (X) = for A < A and F U (X) = \\u\\ 2 H for 
A > A . Hence, A is an eigenvalue of A with eigenvector u G T>(A) if and 
only if E{X)u = for A < A and E(X)u = u for A > A . 

Similarly, one can show that Ao G a c {A) is equivalent to lim e ^o- E{Xq + 
e)u = E(Xo)u for all u G H and E(X\) ^ E(X2) for Ai < A < A2. As already 
proved in sec. 3.3.3, the residual spectrum a r (A) is empty. 

Thus, we have Ao G p(A) if and only if the associated spectral family E 
is constant on a neighbourhood of Ao- This, together with the requirements 
liniA^-oo E(X)u = and liniA^oo E(X)u = u implies that the spectrum of a 
self-adjoint operator is nonempty! 



9.2 Functions of Self— Adjoint Operators 

Let A : T>(A) C H — > H be a self-adjoint operator on a Hilbert space H. 
With the help of the spectral theorem we are now able to define functions of 
A. Again, all the material presented in this section can be found in standard 
textbooks, e.g. [59]. We also mention the freely available book [54]. 



9.2.1 Definitions and Properties 

Definition of f(A) By the spectral theorem there exists a unique spectral 
family {E(X) : A G R} such that A has the representation A = J R XdE(X) 
or, less symbolically, (Au\v)h = J idudfi UtV for all u G T>(A) and v G H 
where fi UjV is the complex-valued measure associated to {E(X)} defined in 
sec. 9.1.1[ Let / : R — > C be a measurable function. Then, we define an 
operator f(A) : V(f(A)) C H - H by V(f(A)) := V(A f ) and f(A) := A f 



where Af is the operator defined in sec. 9.1.2[ i.e. we have 



V(f(A)) = u 



EH: J |/| 2 ^„<oo| 



(f(A)u\v) H = / fdfjL u 



and 



for all u G T>(f(A)) and v G H where, as before, \i u := \i u<a . As already 
mentioned, f(A) is densely defined and linear. 



Calculus for functions of self— adjoint operators One can show that 
E(X)f(A) C f(A)E(X) for any A G R, i.e. f(A) commutes with E(X). 
Furthermore, we have the following functional calculus for f(A) (see e.g. [59] 
for a proof). 
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• V(f(A)) = V(f(A)) and (f(A)u\v) H = {u\f{A)v) H for all u,v G 
V{f{A)). 

• Let g : R — > C be a measurable function. Then, we have the represen- 
tation 

(f(A)u\g(A)v) H = J fgd/j, UtV 
for all u G V(f(A)) and v G V(g(A)). 

• (af)(A)u = af(A)u for all u G ©(/(A)), a G C and (/ + g){A)u = 
f(A)u + g{A)u for all u G V(f(A)) n X>(#(A)). 

• G X>(/(A)) for u G £>(p04)) is equivalent to u G V(fg(A)) and 
we have fg(A)u = f(A)g(A)u for all it G V(fg(A)). 

• If / is finite everywhere on R then the adjoint f(A)* is given by f(A)* = 
J(A) and hence, f(A)f(A)* = f(A)*f(A) which shows that f(A) is 
a normal operator. In particular, /(A) is self-adjoint if / is finite 
everywhere and real-valued. 



9.2.2 Unitary Groups 

Let if be a Hilbert space. A mapping U : R — > £>(i?) is called a strongly con- 
tinuous one-parameter group of linear operators on H if U satisfies U (0) = 
id#, [/ (t+s) = U(t)U(s) for all t, s G R and the mapping 1 1— > U(t)u is contin- 
uous for all u E H. U is called unitary if ||?7(t)-u||// = for all w G and 
any t G R. In particular, every strongly continuous one-parameter group is a 
strongly continuous one-parameter semigroup of linear operators and hence, 
the notion of a generator of a semigroup carries over to groups. The fol- 
lowing proposition shows how to explicitly construct a strongly continuous 
unitary one-parameter group with the help of the functional calculus for a 
given self-adjoint generator. 

Proposition 9.2.1. Let A : T>(A) C H — > H be a self-adjoint operator on a 
Hilbert space H . Define U + (t) := exp(itA) and U-{t) := exp(— it A) via the 
functional calculus. Then, U± have the following properties. 

1. For any t G R, U±(t) is a bounded linear operator on H which satisfies 
\\U±{t)u\\H = \\u\\h for all u G H. 

2. [4(0) = i&h and U±(t + s) = U±(t)U±(s) for all t, s G R. 

3. For any u G H , the function t \— > U±{t)u : R — > H is continuous. 
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4- The limit lim^o h{U±(t + h)u — U±(t)) exists if and only if u G T>(A) 
and in this case we have 

U' ± {t)u := lim U^ + ^ -U^u = ±aJ±{t)Au = ±iAU±{t)u _ 
h— »o h 

Thus, ±iA is the generator of U±. 
Proof. 1. Let t G R. Application of the functional calculus yields 



{U±{t)u\U±(t)u) H = J dfi u = \\u\\ 2 H 



for any u G H where \i u is the measure defined in sec. 9.1.1 for the 
spectral family {E(\) : A G M} associated to the self-adjoint operator 
A. Hence, U(t) : H —>■ H is a unitary linear operator on H. 

2. This follows immediately from the functional calculus. 

3. Fix u G H and let to £ R- Applying the functional calculus we have 

lim ||t/±(t)w— U±(t )u\\ 2 H = lim / | exp(±itA) — exp(±ii A)| 2 <i/i M (A) = 
t-»t t-*t J R 

by Lebesgue's theorem on dominated convergence. 

4. Let u G f Applying the functional calculus we obtain 

\im\\t~ 1 (U±(t)u-u)TiAu\\ 2 H = lim / |t _1 (exp(±itA)-l)=F«A| 2 d// u (A) 



d 
It 



exp(±itA) 



T «A 



i=0 



rf/i u (A) = 



by Lebesgue's theorem on dominated convergence. Hence, we have 
shown U' ± {Q)u = ±iAu for all u G T>(A). However, using the group 
property and strong continuity of U± we infer 

lim h- l {U±{t+h)u-U±{t)u) = U±{t) lim hr\U±(h)u-u) = ±iU±{t)Au 

h—>0 h^O 

and similarly we obtain U±(t)u = ±iAU(t)u. 
Now let A be the generator of z fiU±, i.e. 



P(, I) = { u G H : lim j(U±(t)u - u) exists 
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and A := ^zlim^ \(U±(t)u — u). Then we have 

(Au\v)h = =R limt -1 (u\U±(t)*v — v)h 

= (u\ ±i\imr 1 (U±(-t)v - vj) = (u\Av) H 

V t^O J h 

for all u, v G T>(A) since U±(t)* = U±(—t) by the functional calculus. 
This shows that the operator A is a symmetric extension of A. However, 
since A is self-adjoint there do not exist proper symmetric extensions 
and therefore we have T>(A) = T>(A). 

□ 

Remark 9.2.1. We have shown that iA generates a strongly continuous uni- 
tary one-parameter group if A is self-adjoint. It turns out that the converse 
is also true, i.e. that any strongly continuous unitary one-parameter group 
U on a Hilbert space H is of the form U(t) = exp(itA) for a self-adjoint 
operator A : T>(A) C H — > H. This result is known as Stone's theorem 
(cf. [59]). 



9.3 Well— Posedness of Wave Equations 

Applying the theory outlined in the previous sections we are now able to con- 
sider abstract wave equations on Hilbert spaces. This yields a well-posedness 
result for a certain class of equations and moreover, it is possible to write 
down the solution in a rather explicit form with the help of the functional 
calculus. For more information we refer to [5]. 



9.3.1 Well— Posedness for Strictly Positive Generators 

The first result is the analogue of Theorem |5.1.2| in semigroup theory. Let 
A : T>(A) C H — ► H be a self-adjoint operator on a Hilbert space H satisfying 
(Au\u)h > for all u G T>(A). Then there exists the square root A 1 / 2 of A 



which is again self-adjoint (Theorem 4.1.1 or directly via functional calculus) 



We consider the unitary groups 1 1— > exp(±itA 1 / 2 

Lemma 9.3.1. Let A : V(A) C H — >• H be a self-adjoint operator satisfying 
(Au\u) H > for all u G V(A). For u G V(A) define u(t) := exp(itA 1/2 )u 
via the functional calculus. Then, the function u : R —>■ H satisfies the wave 
equation u" '(t) = —Au(t) for allt el where u'(t) := lim^o j^(u(t+h)—u(t)) . 
Furthermore, the same holds true for u(t) := exp(— itA 1 / 2 )uQ. 
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we 



Proof. Define u(t) := exp{itA l / 2 )uQ for Uq G T>(A). From Prop. 9.2.1 
already know that u'(t) = iA l l 2 u(t) for all tel. Furthermore, expiitA 1 
commutes with A 1 ! 2 . Thus, since Uq G V(A), we have A x ' 2 u G G V(A l/2 ) and 
hence, applying Prop. |9.2.1 again, we obtain 

u"(t) = j^A 1 ' 2 exp(itA^ 2 )u = ^iexpiitA^A^Uo = -Au(t). 

The same reasoning goes through for u{t) = exp^itA 1 ^ 2 )^. □ 

Now we assume the slightly stronger condition (Au\u)h > i(u\u)h for 
all u G T}{A) and some 7 > and consider the operators sin(tA 1/ ' 2 ) and 
cos(ty4 1 / 2 ). According to the functional calculus they are given by 

sin(L4 1/2 ) = -. (exp(itA 1/2 ) - exp(-itA 1/2 )) 

and 

cos(tA 1/2 ) = X - (exp(^ 1/2 ) + exp(-zL4 1/2 )) . 

Furthermore, since (Au\u)h > 7(m|m)_h", it follows that cr(A) C [7,00) (cf. 
Lemma 5.1.5) and in particular, A is bounded invertible. Hence, the operator 
4-1/2 : = A 1/2 A~ 1 : H -> V(A 1/2 ) C H exists and is inverse to A 1/2 . Now we 
have collected all the necessary tools to prove the analogue of the generation 
result Theorem 15. 1.21 

Theorem 9.3.1. Let A : T>(A) C H — > H be a self-ajoint operator on a 
Hilbert space H satisfying (Au\u)h > j{u\u)h for all u G T>(A) and a 7 > 0. 
Then, for given Uq,U\ G T>(A), there exists a unique function u : R — > H 
such that u"{t) = —Au(t) for all t G R and u(0) = u , u'(0) = U\. The 
unique solution u can be given explicitly as 

u(t) = cositA 1 / 2 )^ + smitA^A- 1 / 2 ^. 

Moreover, the real-valued function t 1— > (A l l 2 u(t)\A l / 2 u(t))H + (w '(t)\u' '(t)) h 
("the energy ") is constant for all tel. 



Proof. Let uo, u\ G T>{A). Applying Lemma 9.3.1 it follows immediately that 
u given by 

u{t) = cos(tA 1/2 )u + sm(tA 1/2 )A- 1/2 Ul 

satisfies u"(t) = —Au(i) for all t G R and u(0) = uo, u'(0) = u\. Moreover, 
inserting for u and using the functional calculus we directly compute 

{A l ' 2 u{t)\A l l 2 u(t)) H + (u'(t)\u'(t)) H = WA^uoWl + \\ Ul \\ 2 H 
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which shows conservation of energy. 

Let u : R — > H be another solution with initial data u(0) = Uq and 
u'(0) = U\. Then, u — u is again a solution with zero initial data. However, 
conservation of energy implies that \\A l l 2 {u(t) — u(t))\\ H + \\u'(t)— u'(t)\\ H = 0. 
Note that by definition we have \\u'(t) —u'(t) \\h = ^\\ u (t) — w(£)||# and hence, 

= / "rll M ( s ) _ u(s)\\ H ds = \\u(t) - u - u(t) + uq\\ h 
Jo ds 

which shows u(t) = u(t) for all t G R. □ 

Remark 9.3.1. Note that the solution operators cos(£A 1//2 ) and sm.{tA 1 / 2 )A~ 1 / 2 
are bounded and hence, they can be applied to general functions wo,wi £ H 
and not only to Mo,Mi € £>(^4)- Analogous to semigroup theory this leads to 
the notion of generalized solutions of abstract wave equations. 



9.3.2 Well— Posedness for Nonnegative Generators 

The condition (Au\u)h > 7(m|m)h for a 7 > assumed in the previous section 
was necessary to assure existence of the operator A~ 1 ^ 2 . However, the opera- 
tor A~ x l 2 appears only as sin{tA 1 l 2 )A~ 1 / 2 in the solution formula. Thus, the 
question is whether sui(tA 1 l 2 )A~ l l 2 can be defined reasonably even if A^ 1 ^ 2 
does not exist. This problem is analogous to considering the real-valued 
function / : (0, 00) — > R defined by f(x) := S1 °^ . Since \im x ^ f(x) = 1 
by de l'Hospital, / can be continuously extended to [0, 00). Thanks to the 
functional calculus it is exactly this construction which can be used for the 
operator sin(tA 1//2 )A~ 1 / 2 as well. 

Let A : T>(A) C H — ► H be a self-adjoint operator on a Hilbert space H 
which satisfies (Au\u)h > for all u G T>(A). Then, the self-adjoint square 
root A x l 2 exists and it satisfies (A 1 ^ 2 u\u)h > for all u G T>(A l l 2 ) (Theorem 



4.1.1| ). For t G R we define f t : R -> R by 

, - J sin(tA)A _1 for A > 
/i(AJ : " \ t for A < 

Then, f t is continuous and hence, ft(A 1 / 2 ) is well-defined via the functional 
calculus. Note that the inequali ty (A 1 I 2 u\u)h > for all u G V^A 1 ^ 2 ) implies 



a(A ' ) C [0,oo) (cf. Lemma 5.1.5) and hence, the distribution function 



A 1— > (E(\)u\u)h for w G if is constant on (— 00, 0) where {E(X}} is the 
spectral family associated to A 1 / 2 . Thus, the values of f t on (—00, 0) do not 
contribute. 
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Lemma 9.3.2. The function u : R — > H defined by u(t) := f t (A 1 ^ 2 )uo for a 
given uq G T>(A) satisfies u"(t) = —Au(t) for all t Gl. 

Proof. Let {E(X) : A G M} be the spectral family of A 1 / 2 (spectral theorem) 



and djji u the spectral measure associated to {E(X)} defined in sec. |9.1.1 
Observe that V(f t (A 1/2 )) = H since f t is bounded and therefore, the operator 
ftiyA 1 / 2 ) : H — > H is bounded for any tel. Applying the functional calculus 
we obtain 

lim Wh-^ft+hiA^Uo - ft{A 1/2 )u ) - cos(tA^ 2 )uo\\ 2 H 

h— >0 



= lim / h-\f t+h (X) - f t (X)) - cos(tA) ^(i^ (A) = 
h-*o J 

for any Uq G H by Lebesgue's theorem on dominated convergence and there- 
fore, u'(t) = cos(ty4 1/,2 )w . Analogously, we have 

\\u"(t) + Au(t)\\ 2 H 

POO 

= lim / |/T 1 (cos((t + h)X) - cos(iA)) + X 2 f t {X)\ 2 d^i Uo = 

for any uq G £>(v4) which is the claim. □ 

Remark 9.3.2. In particular it follows that f t (A 1/2 )u G V(A) for all t G R if 

u G 

Now we are able to prove the well-posedness result for nonnegative gen- 
erators. 

Theorem 9.3.2. Let A : T>(A) C H —>■ H be a self-adjoint operator on a 
Hilbert space H satisfying (Au\u)h > for all u G T>(A). Then, for given 
uq, u\ G T>(A) there exists a unique function u : R — > T>(A) C H such that 
u"(t) = —Au{t) for all t G R and w(0) = u , u'(0) = U\. This unique u can 
be given explicitly as 

u(t) = cositA 1 / 2 )^ + ftiA 112 )^ 

where f t ■ R — > R is defined by 

. m J sin(tA)A~ 1 for X> 
/t(A) :_ { t for X <0 

Furthermore, the energy 

t h- (A 1 / 2 ^)!^ 2 ^))* + ( u '(t)| u '(*)) fl 

zs constant. 
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Proof. Let Uq,Ui G T>(A). Applying Lemmas 9.3.2 and 9.3.1 we conclude 
that u, given by u(t) := cos(tA l / 2 )u + f t {A x / 2 )ui, satisfies u"{t) = —Au(t) 
for all t G R and u(0) = Uq, u'(0) = U\. By inserting for u(t) and direct 
computation (functional calculus) we obtain 

{A^u^A^uit^H + {u'{t)\u\t)) H = (A^uolA^uoU + ( Ul \ Ul ) H . 



Uniqueness can be obtained analogously to Theorem 9.3.1 



□ 



Remark 9.3.3. Again, we also have generalized solutions since the solution 
operators cos(L4 1//2 ) and f t (A 1 / 2 ) are bounded and hence, they can be applied 
to any element in H. 



9.3.3 Application to the Linearized Wave Map Prob- 
lem 

Now we return to the equation 

I n « 2 a-P 2 ) 2 i , 2(l-p 2 )cos(2/ )-p 2 I 

0w - (1 - P ) <Pp P P H — 2 = (9.1) 

which describes the (rescaled) linearized flow around the Turok Spergel so- 



lution in hyperbolic coordinates (cf. sec. 7.3). 



Operator formulation As before we define H := L^(0,1) for w(p) = 

2 

(!_f p 2)2 an d Aq : T>(Aq) C H — > H denotes the operator constructed in sec. 
8.1. Hence, the equation u"(t) = —A u(t) for a function u : R — > H is an 



operator formulation of eq. (9.1). 



Well— posedness According to sec. |8.1[ the operator Aq is self-adjoint 
and from Theorem 8.2.1 we know that ct(Aq) = [0,oo). Thus, Aq satisfies 
(A u\u)h > for all u G T>(Aq) (Lemma 4.1.1). Applying Theorem 9.3.2 we 
infer that the Cauchy problem 



u"(cr 



—A u(a) for a > 

w(0) = Mq, u'(0) = U\ 



(9.2) 



for given uq, U\ G T>{A) and a function u : [0, oo) — > H is well-posed. 



(9.2) the energy a t— > ||Aqu((7 



Moreover, Theorem 9.3.2 tells us that for any classical solution of eq. 



\ 2 H + 



U cr 



\ 2 H is conserved. 
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Nonexistence of growing solutions As discussed previously (sec. 8.3), 

1/2 

the mapping u t— > \\A u\\h is a norm on V(A ) and hence, conservation of 
energy implies the nonexistence of growing solutions with respect to u h 
\\Aq u\\h- Thus, there are no solutions of the equation 

A u n 2{ ^-P 2 )\ , 2(1-p 2 )cos(2/q) ^ 
<p ff<r - 2(p a - (1 - p ) (p pp <p p H <p = 

P P 2 

that grow faster than the gauge instability which behaves as e a for a —>■ oo. 
This shows that the Turok Spergel solution is as stable as it can be in the 
hyperbolic coordinates. 



Discussion We conclude that the functional calculus yields essentially the 
same result as the semigroup approach but it is more explicit. In either case 
the main effort lies in determining the spectrum of the operator A . 

Thus, we have exploited the self-adjoint approach (hyperbolic coordi- 
nates) to the linear stability problem and have obtained the best possible 
result. To gain further insight one has to change coordinates. However, 
in a different coordinate system the involved operators become much less 
convenient since they are not self-adjoint anymore. 
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Chapter 10 

The Spectra of A 



This chapter is devoted to the study of the spectra of the operators A n for 
n > 1. We show that the operator A n has exactly n negative eigenvalues 
and give a rough lower bound for the smallest eigenvalue. Furthermore, we 
investigate the spectral behaviour of A n for n — > oo. 

Througout this chapter we adopt the previously used notation, i.e. H : = 
L2,(0,l) with w(p) := (izjw, -4 : V{A) C F -> H is the self-adjoint 
operator generated by the Sturm-Liouville differential expression a given by 
aw := M— (jiu'y+qu) where p(p) := p 2 and g(p) = 2. Moreover, the operator 
A n : V(A n ) C H ^ H is defined by V(A n ) := V{A) and 

A„m := Au + g> n u 

for it e T>(A n ) and 

2(l-p 2 )cos(2/ n (p))-p 2 -2(l-p 2 ) 2 

J4W : = 2 

p 2 

where /„ denotes the n-th self-similar wave map. 



10.1 The Spectrum 



10.1.1 The Continuous Spectrum 



As already mentioned (Theorem 7.2.1), Bizoh has shown existence of the 
smooth self-similar wave maps f n . However, in fact the solutions /„ are not 
only smooth but even analytic on [0, 1], i.e. they can be expanded in a con- 
vergent power series around any point p G [0, 1] (cf. [11]). Furthermore, they 
satisfy / n (0) = and f n (X) = f • Hence, all the results of sec. 8.2 which do 



not depend on the explicit form of fo but solely on the asymptotic behaviour 
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for p — * and p — > 1 carry over to A n without change. In particular, the 
whole Frobenius analysis is equally valid for A n . Thus, we have the following 
proposition. 

Proposition 10.1.1. The continuous spectrum a c (A n ) of A n for any n = 
0,1,2,... is given by a c (A n ) = [0, oo) . 



10.1.2 The Point Spectrum 



As in the case n — 0, the gauge instability (cf. sec. 7.3) is present. Thus, 
the function 6 n defined by 6 n (p) '■= P\A — p 2 f' n (p) is a formal solution of 
A n u = 0. More explicitly, (A — A n )u = reads 



u + 



-u 



(l + A)p 2 -2(l-p 2 )cos(2/„ 



P 



P 2 (l " P 



2\2 



-U 



(10.1) 



According to theorem 7.2.1 the self-similar wave map /„ has exactly n inter- 
sections with the line | on [0, 1). Thus, since / n (l) = §, it follows that 9 n has 
exactly n zeros on (0, 1). Hence, an oscillation argument similar to Lemma 
8.2.2| implies that there are exactly n + 1 numbers = A > Ai > • • • > A n 

A 





such that eq. (10.1) with A 

Uj (0) 



U ; 



sis m sec. 



8.2 



j has a nontrivial solution Uj satisfying 
= (j = 0, 1, 2 . . . , n). According to the Frobenius analy- 
Uj has the asymptotic behaviour Uj(p) ~ p for p — > and 

u j(p) ~ (1 — p) aj for p 1 where aj := 1+ ^~ Aj . Thus, we immediately 
observe that Uj G V(A n ) for j > and hence, Uj is an eigenfunction of A n if 
j > 0. This shows that the operator A n has exactly n negative eigenvalues 
Ai, A2, • • • , A n . Analogous to sec. ^2 it follows that A — A is invertible for 
negative A which are not in the point spectrum and hence, we arrive at the 
following theorem. 

Theorem 10.1.1. The spectrum of the operator A n is given by 
a p (A n ) = {Ai, . . . , A„ G R : > Ai > A 2 > • • • > A n }, 
o~ c (A n ) = [0, 00) and cr r (A n ) = 

for n = 0, 1, 2, . . . . 



10.1.3 A Simple Estimate 

We apply an elementary argument to obtain a rough lower bound for the 
smallest eigenvalue of A n . 



10.2. NUMERICS 
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Lemma 10.1.1. The smallest eigenvalue X n of A n satisfies the estimate 

A„> inf ( 2 ( 1 -^T< 2/ "M)-l)>-oc 
pe(o,i) [ p 2 J 

for any nGN. 

Proof. First of all we show that the infimum exists. We abbreviate 

2(l-p 2 )cos(2/ n (p)) 

K{p) := 1. 

P 2 

Since / n (0) = for all n G N it follows that lim p _, 0+ h n (p) = oo and 
\\m p _>i_ h n (p) = —1. However, h n G C(0, 1) and therefore it must have a 
minimum on [0, 1]. 

Now suppose u G T>(A n ) is an eigenfunction of A n with eigenvalue A < 
inf pg ( 01 ) h n (p). By Frobenius it follows that -u(O) = u(l) = and u'(0) ^ 0. 
Without loss of generality we assume w'(0) > 0. In order to satisfy the 
boundary condition u(l) = 0, u has to have a maximum. Let the first 
maximum be located at p G (0, 1). Thus, we have u(p ) > 0, 'u'(Po) = and 
u "(Po) < 0. Inserting in eq. (10.1) we obtain 

a/ \ ^ — h n (p ) 
u (P°) = — (l - p §)2 u ( p °> > 

which is a contradiction. □ 



10.2 Numerics 

We intend to numerically calculate the point spectrum of the operator A n . To 
this end it is necessary to construct the self-similar wave map f n . As already 
mentioned, this has first been done in [3] with a shooting and matching 
procedure and we will reproduce these results. The point spectrum of A n 
can be obtained by the same technique and this has been done in [10]. 



10.2.1 Construction of Self— Similar Wave Maps 

We construct the self-similar solutions f n with a standard shooting and 
matching technique. It turns out that the derivative f' n {fS) increases very 
quickly with n becoming larger and hence, it is advantageous to use a loga- 
rithmic coordinate. Thus, we define x := log(a + p) for a small a > 0. Eq. 
(7.7) transforms into 

r+^/'- n f f Si tfl 5 = (10.2) 
e x — a (1 — (e x — a) 2 ){e x — a) z 
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where x G [log a, log (a + 1)]. We have the regularity conditions /(log a) 
and /(log(a + 1)) = 





For integrating eq. (|10.2|) we use the ODE solver 



provided by the GNU Scientific Library [24]. Fig. 10.1 shows the first five 
solutions calculated with a = 10~ 4 . 




Figure 10.1: The first five self-similar profiles calculated with a = 10 



10.2.2 Calculation of the Point Spectrum 

We numerically calculate solutions of the eigenvalue equation (A — A n )u = 
which is given explicitly by 



, _ ,. (l + A)p 2 -2(l-p 2 )cos(2/ w , 
u H — u H — — u 



2 

—\ 
P 



p 2 (l -p 2 ) 2 







(10.3) 



for p G [0, 1]. We have a p (A n ) C (— oo,0) by Theorem 10.1.1 and hence, we 



restrict ourselves to A < 0. The boundary conditions at the singular points 
p = and p = 1 are dictated by the requirement u G T>(A n ). According to 
the Frobenius analysis of sec. 8.2, u G T>(A n ) implies u(p) ~ p for p — > and 
u(p) ~ (1 — p) 13 for p — > 1 where f3 = However, since f3 increases as A 

decreases we encounter a technical diffculty: In order to "shoot away" from 
the singular point p = 1 we have to calculate more and more derivatives as 
A decreases. To go around this problem we define a new unknown v(p) := 
(1 — p) _/3 -u(p). It follows that v(p) ~lasp^lif«G T>(A n ) and hence, v is 
analytic if u G T>(A n ). Moreover, the "bad" solution around p = 1 behaves 
as (1 — p) _v ^ and thus, it can easily be distinguished numerically from the 
"good" one since it is singular at p = 1. This is essential in order to obtain 
a well-behaved numerical approximation. Eq. (10.3) transforms into the 
generalized eigenvalue problem 



v" + 



'2 2(1 + //) 



J> 



P 1 



v — 



2cos(2/ n ) 
p 2 (l-p 2 ) 



+ 



;i + /x)(2 + /x) 

1 -p 2 



(10.4) 
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Pn 


fl 


/ 2 


/a 


u 


Hi 


5.333625 


5.304 


5.30 


5.3 






58.0701 


57.68 


57.6 


^3 






625 


620 



Table 10.1: Approximate eigenvalues 

where we have set p := \J — A. 

We solve the eigenvalue problem with the shooting and matching method 
similar to the construction of self-similar wave maps in the previous section: 



Let Vi and v r denote the solutions obtained by integrating eq. (10.4) from 



p = to p = | and p — 1 to p — |, respectively. As a matching condition we 
require the Wronskian W(vi, tv)(§) to vanish. This has the advantage that 
we can fix t>[(0) = 1 as well as v r (l) = 1 and the only free parameter is the 
eigenvalue p. Hence, eigenvalues can be found by a simple bisection and one 
is not forced to use a multidimensional root finder which can be tricky. The 
numerical results are summarized in Table 110.11 

We avoid a detailed error analysis since we do not need these numbers 
in the sequel but merely note the interesting observation that the horizon- 
tal rows seem to converge. This numerical convergence raises the question 
whether there exists a relation like "lim^oo a(A n ) = a{A oo y i where is 



the perturbation operator around the limiting solution = | (cf. sec. 7.2) 
We will discuss this issue in the following section. 



10.3 The Operator 

We study the operator A^ which describes linear perturbations around the 
limiting solution (cf. sec. 7.2). 

The (rescaled) flow of linear perturbations around the limiting solution 
/oo = f is governed by eq. (7.11) with f n substituted by the constant f, i.e. 



pp 



2(1 -P 2 ) 2 7 -2 + P 2 t 
4> P H s — (P 



0. 



(10.5) 

P " P 

Thus, the operator A^ is generated by the formal differential expression a m 
defined by 

1 / d d \ 

(10.6) 



1 

w 



d d 
dp dp 



where w(p) 



r , p(p) := p 2 and q(p) 



-2+p 2 



Thus, we define 



(l_ p 2)2, ■— H ■— (l-p2)2- 

the underlying Hilbert space H as H := L^(0, 1) and apply the method of 



124 



10.3. THE OPERATOR A, 



Frobenius in order to obtain asymptotic estimates for solutions of a^u = 0. 
Around p — 1 the indices of the Frobenius analysis are both equal to | and 
thus, the situation is completely analogous to A n and p = 1 is in the limit- 
point case. However, around p = the indices are — 1±l ^ which means that 
all nontrivial solutions of a^u = belong to H near p = 0. Thus, p = is 
in the limit-circle case and in order to define a self-adjoint operator we have 
to specify a boundary condition at p = 0. But which boundary condition is 
the "correct" one? Bizoh [7] has studied the analogous problem for Yang- 
Mills equations and he proposed to choose the boundary condition in such a 
way that there exists a nontrivial solution of a^u = similar to the gauge 
instability for A n . 



10.3.1 The Boundary Condition 

For A < let u(-, A) be the solution of (A — a 0O ) , w(-, A) = with asymptotic 
behaviour u(p, A) ~ (1 — p) a for p — > 1 where a = (which exists by 

Frobenius' method). We are interested in the asymptotic behaviour of «(•, A) 
for p — > 0. According to Frobenius' method we have 

u(p, A) oc p^ 2 (p^(p, A) + m(A)p-^/ 2 ^(p, A)) 

where <p(-,A) as well as if)(-,\) are holomorphic around p = and satisfy 
<p(0, A) = ^(0, A) = 1. The complex number m(A) is called the connection 
coefficient. The asymptotic behaviour of u(-,X) can be written in a more 
convenient way. 

Lemma 10.3.1. For A < let u(-,X) be a real solution of (A — a^u = 
with asymptotic behaviour u(p, A) ~ (1 — p) a for p — > 1 where a = . 
Then, u(-,X) can be written as 

u(p, A) = cp" 1 ' 2 ^sin logp + 6(X) j + %, A) 

where c is a real constant, |5(A)| < | 7 h(-,\) e C 1 (0, 1) and 

lim /i(p, A) = lim ph'(p, A) = 0. 
p^o p— >o 

Proof. As already mentioned above, we have 



A) oc p^ 2 [ P ^(p, A) + m(A)p-^/ 2 ^(p, A)) 
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Since u(-, A) is real by assumption, we have u(-, A) = Rew(-, A) and hence, we 
infer 

p 1/2 u(p, A) = a(A) cos lo S PJ + b W sin lo g PJ + h (P, A) 

where a(A) and 6(A) are real-valued and defined by m(A). Moreover, h(-, A) G 
C 1 (0,1), linip_^o h(p, A) = and lim p ^o ph'{p, A) = which follows immedi- 
ately from the properties of <p(-, A) and A). Thus, the identity acosx + 
fesinx = \Ja 2 + 6 2 sin(x + <i) where <i = arctan - finishes the proof. □ 

According to Bizoh's proposal, eigenfunctions should satisfy 5(A) = 5(0). 
The following lemma shows how this requirement can be translated into a 
boundary condition. 

Lemma 10.3.2. For A < letu(-,X) be a real solution of (X — a^u = with 
asymptotic behaviour u(p, A) ~ (1 — p) a where a = Assume further 

that |5(0)| < | where 5 is defined in Lemma 10.3.1 Then, 5(A) = 5(0) is 
equivalent to 



lim p 2 (u(p, X)u'(p, 0) - u\p, X)u{p, 0)) = 

p^0 



where ' :— 4-. 

dp 



Proof. Applying Lemma |10.3.1 and using the addition theorem for trigono- 



metric functions we readily obtain 

lim p 2 (u(p, X)u'(p, 0) - u'(p, X)u(p, 0)) = sin(5(A) - 5(0)) 

p^0 

which yields the claim since |5(A) — 5(0) | < n. □ 



10.3.2 Construction of the Operator 

We construct the self-adjoint operator Aoo by imposing the boundary con- 
dition discussed in the previous section. 

To this end recall the definition of the maximal operator A max : T>(A max ) C 
H — > H generated by which is given by 

^(Amax) := {li 6 : u,pu' G ACi oc (0, 1) eH} 

and A max ti := a^u for u G V(A max ). 

We denote by \ a real-valued nontrivial function satisfying a^x = with 
asymptotic behaviour x(p) ~ (1 — p) 1 ^ 2 for p — > 1. Such a function exists by 



126 



10.3. THE OPERATOR A 



Frobenius' method and it is unique up to constant multiples since the other 
linearly independent solution of a^u = around p = 1 contains a logarithmic 
term. Note, however, that \ H due to its asymptotic behaviour for p — > 1. 
Using a smooth cut-off function we construct a x G V(A max ) satisfying 
x{p) = x{p) f° r P e (0) §)■ Obviously, there exists a i> G £>(A max ) with 
[x? t ']p(0) 7^ 0. We define the operator by 

P(A 00 ) := {u G D(A max ) : [u,x] P (0) = 0} 



and AqqU := a^u for m G V(A oq ). Invoking Lemma 4.2.12 we immediately 
infer that is self-adjoint. 



10.3.3 Calculation of the Point Spectrum 

Consider the formal eigenvalue equation (A — a^u = 0. We make a coordi- 
nate transformation p i— »■ z := p 2 and define a new unknown 



where a := — 1(— 1-MV<) and /3 := — 1(1 + V— A). The equation (A — a^u = 
transforms into 

z{l - z)v" + [c - (a + 6 + l)z]v' -abv = (10.7) 
where a := \(l + 2 V /Z A + i>/7), & := i( 3 + 2 v /Z ^ + and c := 1 + 



Eq. (10.7) is the hypergeometric differential equation (cf. [1]). Around z — 1 
there exist two linearly independent solutions t>i(-,A) and t>i(-,A) given by 
Wl (z, A) := 2 Fi(a, 6; a + 6 + 1 - c; 1 - z) and £i(z, A) := (1 - z) c - a - b 2 F l {c - 
a, c — b; c + 1 — a — b; 1 — z) where 2 -Pi denotes the hypergeometric function 
(cf. [1]). We are interested in eigenfunctions and hence, the condition u G H 
rules out the solution V\(-, A) since 2-Pi(c — a, c — b; c + 1 — a — b; 1 — z) ~ 1 
for 2; — > 1. Around z = the two linearly independent solutions t>o(-, A) and 
v (-, A) are given by v (z, A) = 2-^1(05 6; c; z) and {?o(^, A) = 2; 1_C 2 -Fi(a + 1 — 
c, b + 1 — c; 2 — c; z) . By the general theory of linear second order ordinary 
differential equations we infer that for z G (0, 1) we can write 

ui(z, A) = u (z, A) + m(A)i5o(z, A). 

The connection coefficient m(A) can be given explicitly in terms of the V- 
function [18] and reads 

/x . r(a + 1 - c)r(6 + 1 - c)r(c - 1) 

m(AJ = 



r(o)r(6)r(i - c) 
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fl 


h 


fa 


h 


foo 


/il 


5.333625 


5.304 


5.30 


5.3 


5.3009 


A*2 




58.0701 


57.68 


57.6 


57.637 


V>3 






625 


620 


619.61 



Table 10.2: Approximate eigenvalues 



The boundary condition specified in the definition of T>(A 0O ) translates into 
the requirement m(A) = m(0). This transcendental equation can be solved 
numerically using bisection and by this we ob tain the point spectrum cr^Aoo) 



10.2 



where we have set \i :- 



of Aoq. The results are given in Table 
Furthermore, we have duplicated Table 10.1 for comparison. 

Hence, the numerical results suggest that the point spectrum of A r , 
verges to dp^A^) for n — > oo. 



con- 
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Chapter 11 

Further Results and Outlook 



So far we have studied linear stability of the Turok Spergel solution f in hy- 
perbolic coordinates. We have identified two features which spoil the analysis 
to a certain degree. First, the time translation symmetry of the original prob- 
lem leads to an exponentially growing solution of the perturbation equation. 
This fact corresponds to an isolated point in the continuous spectrum of the 
operator A . However, one could easily go around this problem by defining 
a projection operator which removes this point from the spectrum of A . 
By this, one would immediately gain linear stability of / as conjectured. 
The reason why this does not work is the nature of the coordinate system 
which induces an unbounded part in the continuous spectrum of A . Thus, 
the only possibility is to introduce a new time coordinate r different from 
a. However, this destroys the self-adjoint character of the problem which 
makes it much more difficult. To illustrate the problems one has to deal with 
we consider the (possibly simplest) choice r := — log(T — t). The resulting 
linarized equation around the self-similar solution f n reads 

TT - (1 - p 2 )0 pp + 2 P rp + r - 2 -^jA<P P + 2C ° S p ( 2 2/w) = 0. (11.1) 

Thus, one obtains a mixed derivative since the coordinate lines are not or- 
thogonal anymore. Such a term turns out to be very inconvenient. We 
formally rewrite this evolution equation as a first order system of the form 

^u(t) = L„u(t) 

for a matrix differential operator L n . Now we consider the eigenvalue equa- 
tion (A — L )u = which reduces to 
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Since /o(p) = 1 \\+ p ty2 this equation has six regular singular points p = 
0, ±1, ±2, oo. Using the transformation p i— > z := p 2 one is left with the four 
singularities z = — 1,0, 1, oo. Thus, solutions are given in terms of Heun's 
functions (cf. [43]). However, most of the knowledge concerning Heun's func- 
tions is based on numerical techniques and thus, it is very difficult to obtain 
the "spectrum" of this generalized eigenvalue problem. Furthermore, it is by 
no means clear what boundary conditions one should specify at the singular 
points p = and p — 1. In the self-adjoint formulation the choice of the 
function space is dictated by the operator itself but in this non-self-adjoint 
setting we do not have such an information. Thus, as a first step one would 
have to find a function space such that Lq is properly defined and the linear 
evolution problem eq. (11.1) is well-posed which is probably not so easy. 
Secondly, one has to analyse the spectrum of L , i.e. study eq. (11.2) with 
appropriate boundary conditions. There exist partial results addressing this 
issue by the author [16] which state that eq. (11.2) has no analytic solutions 
for real A unless A = 1. Furthermore, this problem has been investigated 
numerically (cf. [8], [15]) and these studies strongly suggest that A = 1 is 
indeed the only "eigenvalue" with positive real part. We note that the origin 
of this unstable mode is well understood since it stems from the time transla- 
tion symmetry of the wave map equation similar to the gauge instability we 
have encountered in the self-adjoint formulation. Thus, if one could make 
these ideas rigorous it would be possible to use a spectral projection which 
removes the eigenvalue 1 from the spectrum of L and prove linear stability 
of/ . ' 



Of course, the ultimate goal would be proving nonlinear stability of this 
solution. We do not try to make this precise but merely note that quite re- 
cently [31] rigorous results have been obtained for a different (easier) problem 
where some of the difficulties arising there seem to be similar to our wave 
map model. Concerning self-similar blow up for semilinear wave equations 
we also mention [22] and [36], [35] although the latter deal with the energy 
(sub)critical case. There is also a notion of orbital stability pioneered by 
Weinstein (cf. e.g. [57]) which might be able to deal with the difficulties 
arising from the fact that fo is a one-parameter family of functions rather 
than a single solution. However, this approach has mainly been worked out 
for the Schrodinger equation which is quite a different problem. To conclude 
we have to admit that a rigorous proof of nonlinear stability of fo seems to 
be beyond the scope of present techniques. 

Finally, we emphasize that the results obtained in this thesis are not 
confined to the particular wave map model which has been considered. In 
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principle, they are equally applicable to any evolution equation of the form 

if it shows similar behaviour (existence of self-similar solutions). A prominent 
example is the Yang-Mills field in 5 + 1 dimensions which is in the same 
criticality class as our wave maps model and analogous phenomena have 
been observed for this system. Bizoh has studied this model [7] and he has 
proved existence of a countable family of self-similar solutions. The analysis 
presented in this thesis carries over to the Yang-Mills model with minor 
changes. 
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Appendix A 



Symbols 



N The natural numbers 1,2,... 

N The natural numbers including zero 0, 1, 2, . . . 

Z The integer numbers 0, —1, 1, —2, 2, —3, 3, . . . 

R The real numbers 

R k {(x!,...,x k ) :xj GRV1 < j < k} 

(a, b) {x G R : a < x < b} 

[a,b) {x G R : a < x < b} 

(a,b\ {x G R : a < x < b} 

[a, b] {x G R : a < x < b} 

C The complex numbers 

C k {(x 1 ,...,x k ):x j eCVl<j<h} 

x The complex conjugate of x G C 

Rex The real part ofxeC 

Imx The imaginary part of x G C 

Mj., <9 X M, t^u Partial derivative of u with respect to x 

idM The identity on a set M, i.e. idM(^) := £ for x G M 

ker(A) The null space or kernel of a linear mapping A : V ^ W 

between two vector spaces V and W, 

i.e. ker(A) := {n 6 V : Au = 0} 

im(A) The image or range of a mapping A : X — > F between 

two sets X,F, i.e. im(A) := {Ax : x G X} 
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